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1.  INTRODUCTION/SCOPE  OF  THIS  RESEARCH 
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With  the  increase  in  use  of  digital  communications  has  come  a 
renewed  interest  in  nonlinear  filtering  techniques.  In  some  signal 
processing  applications  the  suppression  of  noise  on  a  discrete-time 
signal  sequence  can  not  be  adequately  accomplished  using  linear 
filtering.  In  order  to  preserve  sharp  edges  in  a  signal  while  still 
smoothing  noise  and  removing  impulsive  components  requires  some 
type  of  nonlinear  or  adaptive  filtering. 

The  most  common  nonlinear  smoother  to  be  investigated  is  the 
median  filter.  In  this  filter  the  median  is  taken  from  a  finite 
length  window  of  values  surrounding  the  input  value.  This  scheme 
has  been  shown  to  preserve  edges  in  signals  with  a  minimum  of 
distortion  and  is  resistant  to  outliers  in  the  data.  However,  it  often 
does  not  provide  adequate  smoothing  of  noise  that  is  not  by  nature 
impulsive.  Yet  the  desirable  properties  of  median  filtering  make  it 
worth  investigating. 

Several  methods  have  been  investigated  to  overcome  the 
limitations  of  median  filters  while  retaining  some  of  its 
advantageous  properties.  These  methods  have  sought  to  retain  the 
impulse  rejection  and  edge  preservation  properties  of  the  median 
filter  while  improving  the  Gaussian  noise  suppression.  Several  of 
the  more  prominent  filters  are  discussed  in  Section  II.  This  section 
shows  the  progression  of  hybrid  filter  designs  combining  the 
desirable  properties  of  linear  and  nonlinear  filters. 

Yet  this  progression,  while  impressive  and  important,  does  not 


NOMENCLATURE 


Throughout  the  course  of  the  simulation  several  different 
methods  were  used  to  describe  the  various  filters  and  filter 
operations  under  consideration.  First,  median  and  mean  refer  to 
the  two  filtering  operations  denoted  by  those  words:  median  and 
mean  filters.  However,  the  word  median  is  also  used  in  several 
titles  to  denote  transfer  functions  that  were  calculated  using  the 
median  method  described  in  Appendix  A.  Transfer  functions  not 
calculated  using  the  median  method  were  calculated  by  averaging. 

The  a- trimmed  mean  and  a-trimmed  linear  filters  are  denoted 
by  either  a-TM  and  a-TL  or  a-TM  and  a-TL.  The  size  and 
number  of  elements  trimmed  is  also  often  specified  in  most  titles. 
An  a-trimmed  mean  filter  of  length  9  with  4  elements  trimmed  is 
specified  by  a-TM  9T4.  Note  that  a-TM  9T0  is  the  same  thing  as 
mean  9  or  mean  L9  and  that  a-TM  9T8  is  the  same  thing  as 
median  9  or  median  L9. 

The  a-trimmed  linear  filters  are  specified  in  the  same  manner 
with  an  additional  designation  to  show  which  linear  filter  is  being 
used.  The  linear  filter  designation  standing  alone  refers  to  the 
untrimmed  linear  filter.  For  example,  LPF1500  L31  refers  to  a  low 
pass  filter  of  31  elements  with  a  cut  off  frequency  of  1500Hz 
(based  on  a  10kHz  sampling  rate).  Equivalent  designations  for  this 
same  filter  are  LPF.15  L9  and  LPF1500  9.  The  numbers  following 
the  three  letter  filter  designation— LPF,  DPF,  or  HPF— indicate  a 
frequency  associated  with  that  filter.  For  the  low  pass  and  high 


pass  filters  this  number  represents  the  frequency  in  the  center  of 
the  transition  band,  the  cut  off  frequency.  For  the  band  pass  filter 
the  number  represents  the  frequency  in  the  center  of  the  pass 
band.  The  frequency  may  be  specified  as  a  normalized  frequency 
from  0  to  0.5  or  as  a  frequency  in  Hertz  based  on  a  10kHz 
sampling  rate.  Designations  for  an  a- trimmed  linear  filter  look  like 
a-TL  BPF2500  31T2.  This  designation  is  for  a  31  element  band  pass 
filter  with  a  pass  band  centered  at  2500Hz  having  2  elements 
a- trimmed. 

There  are  some  special  one  letter  designations.  BPF2500s  is  used 
to  identify  the  "skinny”  BPF  centered  at  2500Hz  with  a  sharper 
roll  off  than  BPF2500.  The  following  one  letter  designations  refer  to 
the  filters  opposite  them: 

X(tra  low  pass)  LPF300 

L(ow  pass)  LPF1500 

B(and  pass)  BPF2500s 

F(at  band  pass)  BPF2500 

H(igh  pass)  HPF3500. 

These  designations  were  used  when  a  condensed  notation  was 
necessary.  They  are  often  seen  as  F31T2  to  designate  a-TL 
BPF2500  31T2. 
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address  a  particularly  important  characteristic  of  linear  filters: 
frequency  selectivity.  There  are  many  applications  in  which  the 
signal  of  interest  is  not  obtained  using  a  simple  low  pass  filter*— the 
median  filter  being  a  low  pass  operation— but  rather  a  filter 
selecting  some  other  frequency  band.  It  is  addressing  these 
applications— matched  filter  banks  and  IF  and  RF  filters,  for 
example— to  which  this  research  effort  is  directed. 

It  is  a  relatively  simple  matter  to  design  a  linear  filter  to  pass 
a  certain  band  of  frequencies.  However,  if  the  environment 
through  which  the  signal  is  propagated  happens  to  be  impulsive 
rather  than  Gaussian,  say  under  the  polar  ice  cap  as  an  extreme 
example,  then  a  linear  filter  may  not  provide  adequate  smoothing 
of  the  noise.  In  this  case  it  would  be  desirable  to  have  some  of  the 
outlier  resistance  properties  of  the  median  filter  in  order  to  deal 
with  the  impulsive  noise. 

Section  HI  describes  the  generalized  order  statistic  filter  (GOSF) 
model.  An  a-trimmed  linear  (a-TL)  filter  is  defined  using  this 
model  which  is  an  attempt  to  meet  the  goal  of  a  frequency 
selective  filter  having  the  outlier  resistance  properties  of  the 
nonlinear  median  filter.  In  this  section  the  a-TL  filter  model  is 
defined  based  on  the  logical  progression  of  hybrid  filters  described  in 
Section  II.  Several  important  considerations  of  this  design  are  also 
discussed  at  the  end  of  Section  III. 

In  Section  IV  some  performance  measures  are  defined  with 
which  the  a-TL  filter  will  be  evaluated.  For  this  filter  model  the 
performance  measures  of  interest  are,  of  course,  outlier  resistance 


and  frequency  selectivity.  Special  attention  is  given  to  how  the 
frequency  performance  of  the  nonlinear  a-TL  filter  is  assessed.  A 
brief  discussion  of  some  of  the  important  simulation  parameters  is 
also  included. 

The  results  of  the  simulation  are  presented  in  Section  V. 
Several  effects  of  the  trimming  process  on  the  linear  filter 
characteristics  are  discussed  to  provide  insight  into  the  filtering 
mechanism.  More  importantly,  several  specific  effects  of  filter 
coefficient  parameters  on  the  performance  of  the  a-TL  filter  are 
discussed  in  detail.  This  discussion  centers  on  how  varying  certain 
filter  design  parameters  such  as  the  number  of  elements  in  the 
filter  and  the  steepness  of  the  designed  filter  rolloff  seem  to  affect 
the  performance  of  the  a- trimmed  filters  in  general. 

The  analysis  of  these  effects  is  summarized  in  Section  VI.  The 
overall  performance  of  the  a-TL  filter  is  assessed  and  possible 
avenues  for  improvement/future  research  on  the  trimmed  linear 
filter  are  discussed. 


II.  THE  ORDER  STATISTIC  FILTER 


A.  THE  L  FILTER 

Bovik,  Huang,  and  Munson  introduced  a  generalization  of  the 
median  filter  which  they  called  an  order  statistic  filter  (OSF) .  This 
filter  uses  some  linear  combination  of  the  ordered  data  in  a 
window  around  an  input  point  to  produce  the  output.  With  the 
proper  choice  of  coefficients  one  may  create  a  running  mean  filter 
(all  the  coefficients  equal  to  l/N),  a  median  filter,  or  something  in 
between.  In  addition,  a  minimum  or  maximum  filter  is  also  a 
special  case  of  the  OSF.1  Lee  and  Kassam  noted  that  an  OSF  is  the 
same  thing  as  an  L  filter  since  an  L  filter  is,  by  definition,  a 
linear  combination  of  order  statistics.2 

Lee  and  Kassam  went  on  to  develop  a  simple  representation  of 
the  L  filter  called  an  a- trimmed  mean  (a-TM)  filter.  In  this 
representation  the  filter  has  a  requirement  that  all  the  weights  in 
some  central  portion  of  the  window  be  equal  to  a  constant  with 
the  rest  of  the  weights  taken  to  be  zero.  The  number  of  samples 
trimmed  off  each  end  of  the  window,  T,  is  parameterized  by  a, 

T  =  |_a(2N+l)J,  where  0£a£0.5  and  LXJ  is  the  largest  integer  less 
than  or  equal  to  X.  This  filter  also  has  the  two  special  cases  of  the 
running  mean  filter,  a=0,  and  the  median  filter,  a=0. 5.  The 
parameter  a  is  chosen  based  on  a  priori  knowledge  of  the  noise 
distribution  or  by  an  adaptive  scheme.3 


The  output  yk  for  this  filter  is  given  by 

2N+1-T 

yk=  2  x*0)/[2(N-T)+l],  (: 

j-T+1 

where  is  the  jth  order  statistic  of  the  window  centered 
around  the  IIth  element  in  the  signal.  Lee  and  Kassam  observed 
that  the  o-TM  filter  possesses  a  clear  trade-off  between  the 
advantages  of  a  mean  and  a  median  filter.  As  a  approaches  0  or 
0.5  the  characteristics  of  the  filter  approaches  a  mean  or  median 
filter,  respectively.4 

B.  M  FILTER 

Lee  and  Kassam  also  proposed  an  M  filter  with  a  more 
favorable  combination  of  the  two  individual  filters’  characteristics. 
The  output  yk  of  an  M  filter  is  defined  as  the  solution  to  the 
equation 

k+N 

Z  0(xryk)  =  o,  (: 

i-k-N 

where  0  is  some  odd,  continuous,  and  sign- preserving  function. 
When  0  is  the  linear  function  0(x)  =  ax,  for  a= constant,  the  M 


6 


filter  reduces  to  the  running  mean  filter.  On  the  other  hand,  the 
M  filter’s  characteristics  approach  those  of  the  median  filter  as 
0(x)-+sgn(x)  under  certain  conditions.5  They  also  defined  a 
standard  type  M  filter  (STM  filter)  as  a  limiter  type  filter  for 
which 


r 

0(x)  =  * 


g(p),  x  >  P 

g(x),  M^P  (3) 

-g(p).  x  <  -p 


and  g(x)  =  ax  as  shown  in  Figure  1. 

•60 


Figure  1:  0(x)  for  a  STM  filter6 


The  STM  filter  was  found  to  have  very  favorable 
characteristics.  Its  window  is  somewhat  data  dependent  so  that  it 
tends  to  treat  as  outliers  those  values  which  are  very  large  or 
very  small  as  compared  to  the  median.  This  leads  to  the 
interpretation  of  an  STM  as  a  data  dependent  type  of  L  filter  that 
may  have  a  nonsymmetric  window.  Lee  and  Kassam  concluded 
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that  an  STM  tended  to  act  as  a  running  mean  filter  when  neither 
edges  nor  impulsive  noise  were  present  and  as  a  median  filter  over 
areas  with  edges.  It  combined  the  advantages  of  the  two  types  of 
filters  more  favorably  than  an  L  filter.  The  major  disadvantage  of 
the  STM  is  that  it  is  difficult  to  evaluate  the  output  of  the  filter. 

In  their  simulation,  Lee  and  Kassam  used  the  iterative  Newton's 
method.7 

Finally,  Lee  and  Kassam  noted  that  the  main  difference 
between  an  a-TM  filter  and  an  STM  filter  was  that  the  number  of 
samples  trimmed  from  each  end  of  the  window  of  an  STM  was 
data  dependent  and  not  necessarily  symmetric.  This  was  the 
reason  STM  filters  could  outperform  the  simple  a-TM  filter. 

Applying  this  observation  they  came  up  with  a  modified  trimmed 
mean  (MTM)  filter  which  is  simple  to  implement  and,  in  many 
cases,  produces  results  at  least  as  good  as  those  obtained  with  an 
STM  filter.* 

The  MTM  filter  first  determines  the  sample  median  mk  inside 
its  window  and  then  chooses  an  interval  [mk  -  q,  mk  +  q]  using 
some  preselected  constant  q.  All  data  samples  whose  value  lie 
outside  the  range  are  discarded  and  the  average  of  the  remaining 
values  is  taken  as  the  output  at  sample  k.  This  is  a  similar 
operation  to  an  a-TM  filter  only  for  the  MTM  filter  the  range 
around  the  median  is  the  determining  factor  in  how  many  samples 
are  discarded  instead  of  the  constant  number  of  discarded  samples 
determining  the  range  of  values  to  average  in  the  a-TM  filter.  The 
MTM  filter  also  differs  from  the  STM  filter  in  that  the  MTM  filter 
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may  discard  the  values  that  lie  outside  the  preselected  range 
whereas  the  STM  filter  only  limits  these  values  to  the  range  itself.5 6 7 8 9 
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111.  THE  GENERALIZED  ORDER  STATISTIC  FILTER 

A.  DEFINITION 

Seeing  these  developments  and  the  improvements  that  they 
allow  makes  one  wonder  at  the  possibilities  of  merging  the 
capabilities  of  linear  and  nonlinear  filtering  techniques.  Using  the 
notation  that  has  been  discussed  up  to  this  point,  we  will  look  at 
the  possibilities  of  a  filtering  scheme  proposed  by  Ritcey  that  may 
be  called  a  generalized  order  statistic  filter  (GOSF) .  The  scheme 
works  like  this:  a  signal  is  passed  through  a  set  of  time  weighting 
coefficients  designed  for  some  particular  frequency  response.  But 
before  the  results  are  summed,  the  values  are  ordered  and 
trimmed  according  to  some  scheme.  This  is  easier  to  visualize 
graphically  as  shown  in  Figure  2. 


filter 

output 


Figure  2:  Generalized  Order  Statistic  Filter 


A  notable  feature  of  this  system  is  that  many  other  filters  are 
special  cases  of  this  generalized  model.  As  with  the  L  filter  in 
general  and  the  a-TM  filter  specifically,  the  running  mean  and 
median  filters  are  easily  and  obviously  attained  by  choosing  the 
proper  coefficients  for  this  model.  In  addition,  one  notes  that  the  L 
filter  is  also  a  subset  of  this  general  model  if  we  choose  all  of  the 
an’s  to  be  unity.  Finally,  one  can  also  obtain  an  MTM  filter  by 
providing  a  mechanism  to  check  the  signal  values  as  they  come 
out  of  the  sorter  and  choose  the  bn  coefficients  based  on  the 
median  and  a  preselected  range  parameter . 

The  output  zk  is  given  by 

2N+1-T 

=  (1/g)  2  bjyk(j)  (4) 

J«T+1 

where  ykQ)  =  ordered  [yj  centered  around  xk  and 

yt  *  a^  ;  i  =  1,2, . .  .N  (5) 

and  g  is  the  gain  factor  to  be  discussed  later. 


B.  THE  a- TRIMMED  LINEAR  FILTER 


The  a-trimmed  linear  (a-TL)  filter  is  obtained  from  the 
generalized  order  statistic  filter  by  using  a  set  of  linear  filter 
coefficients  as  the  time  weighting  coefficients.  The  signal  enters  a 
window  and  is  weighted  by  the  an’s  just  as  if  it  were  going  to  be 
linearly  filtered.  Only  before  the  weighted  values  are  summed  to 
produce  the  filter  output  as  in  a  linear  filter  operation,  they  are 
sorted  and  trimmed  (the  bn’s  being  restricted  to  zero  or  one) . 

Thus  we  have  a  sort  of  trimmed  “linear”  filter  similar  to  the 
a-trimmed  mean  filter.  The  main  difference  between  the  two 
filters  is  that  the  a-TM  filter  has  a„’s  that  are  a  constant  equal  to 
1/(N-T)  while  the  a-TL  filter  has  a  set  of  FIR  coefficients  for  the 
a^s  which  are  modified  by  a  constant,  the  gain  factor  g. 

C.  GAIN 

A  serious  consideration  in  this  model  is  the  gain  of  the  system. 
Without  some  modification,  if  a  set  of  linear  FIR  filter  coefficients 
that  sum  to  a  particular  value  (say  unity  for  a  low  pass  filter) 
are  used  and  then  some  of  the  values  are  trimmed  off,  there 
would  be  some  gain  in  the  system  that  may  make  it  difficult  to 
see  what  is  actually  happening.  To  take  into  account  the  gain  of 
the  system  we  look  at  what  happens  to  a  constant  signal  passing 
through  it.  To  overcome  the  attenuation  the  linear  set  of 
coefficients  are  ordered  and  those  coefficients  that  would  not  be 


trimmed  by  the  filter  are  summed.  As  an  example,  suppose  a 
constant  signal  of  unity  was  passed  through  a  hypothetical  seven 
element  LPF  whose  coefficients  were  as  follows: 


ans  =  16  16  16  16  16  16  16 
where  E  aB  =  1. 


Ordering  these  we  would  get  -1,  -1,  2,  2,  4,  4,  and  6  sixteenths. 
If,  in  our  filtering  scheme,  we  happened  to  trim  off  the  largest 
and  smallest  values  and  summed  the  result  we  would  get  11/16  as 
the  filtered  output  of  the  constant  unity  input.  Thus  our  gain 
needs  to  be  16/11  which  corresponds  to  summing  the  a^s  that  are 
not  trimmed  and  dividing  the  filter  output  by  this  amount .  The 
gain  factor  is  given  by 


T/2 

g  =  2  an  -  E  (a^)  +  a(j|+i-|))  (6) 

1*1 


where  the  a^  =  ordered[an]  and  T  is  the  number  of  elements 
being  trimmed  (always  even).  In  general,  however,  a  filter 


IV.  PERFORMANCE  EVALUATION 


A.  PERFORMANCE  MEASURES 

In  evaluating  this  proposed  filtering  scheme  concentration  will 
be  focused  on  two  major  filter  characteristics.  Our  primary 
interest  will  be  trying  to  maintain  the  frequency  selectivity  of  a 
linear  filter  (i.e.  being  able  to  create  a  filter  to  pass  a  certain 
frequency  band)  while  adding  in  the  inherently  nonlinear 
characteristic  of  outlier  resistance.  These  two  capabilities  do  not 
exist  simultaneously  in  any  of  the  filters  discussed  previously.  In 
all  of  the  filter  models  discussed  up  to  this  point  the  application 
was  assumed  to  be  that  of  a  low  pass  filter.  What  if  the 
application  calls  for  a  selection  of  some  other  frequency  band  as 
the  signal  of  interest?  In  this  case  a  linear  filter  is  needed,  but  we 
already  know  that  a  linear  filter  is  highly  susceptible  to  the  effects 
of  outliers  on  the  input  signal.  Is  there  a  filter  design  that  will 
allow  a  frequency  selection  of  other  than  a  low  pass  filter  while 
providing  improved  performance  in  the  presence  of  impulsive  noise? 
We  are  going  to  evaluate  the  proposed  a-trimmed  linear  filter  in 
regards  to  these  two  performance  measures:  outlier  resistance  and 
frequency  selectivity. 

FREQUENCY  SELECTIVITY 

The  evaluation  of  the  frequency  selectivity  of  a  filter  is  going  to 
be  based  on  how  well  the  modified  filter's  frequency  characteristics 
match  that  of  the  unmodified  linear  filter.  This  criteria  was  chosen 


because  it  is  relatively  simple  to  evaluate.  It  does  not  require  any 
modification  of  the  filter  coefficients.  We  can  just  compare  how 
well  the  frequency  response  of  the  modified  filter  compares  to  that 
of  the  unmodified  filter.  This  allows  an  evaluation  of  how  the 
ordering  and  trimming  of  the  a- trimmed  linear  filter  affects  the 
frequency  characteristics  of  the  original  linear  filter. 

We  could  have  chosen  to  set  a  particular  frequency 
characteristic  as  the  desired  response  and,  through  some  iterative 
procedure,  changed  the  coefficients  of  the  filter  adaptively  to  best 
imitate  this  desired  response.  However,  due  to  the  nonlinear 
characteristics  of  the  a-TL  filter  this  could  prove  difficult  to  do. 
Therefore,  this  possibility  was  not  explored.  This  could  be  an  area 
for  further  study. 

In  evaluating  the  frequency  selectivity  of  the  a-TL  filter  we 
will  first  compare  visually  the  frequency  response  characteristics  of 
the  a-TL  filter  to  those  of  the  original  linear  filter.  This  may  allow 
a  quick  insight  into  the  actual  physical  process  of  the  filter  and 
may  suggest  ways  to  improve  filter  performance.  Second,  we  will 
do  an  error  analysis  of  the  modified  filter’s  frequency 
characteristics  compared  to  those  of  the  original  linear  filter.  This 
will  allow  a  quantitative  evaluation  of  how  much  degradation  in 
frequency  performance  is  introduced  by  the  a- trimming  process. 
TRANSFER  FUNCTIONS 

A  problem  encountered  here  is  how  to  measure  the  frequency 
characteristics  of  a  particular  a-TL  filter.  The  frequency  response 
of  an  FIR  filter  is  found  simply  by  using  an  impulse  as  the  input 
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to  the  filter  and  taking  the  DFT  of  the  output.  This  is  the  transfer 
function  of  the  filter  because  the  DFT  of  the  input  (an  impulse)  is 
unity  at  all  frequencies.  However,  we  cannot  find  the  transfer 
function  of  the  a-trimmed  linear  filter  in  the  same  manner 
because  it  is  designed  to  reject  impulses  in  the  input  signal.  By 
using  an  impulse  as  the  input  into  an  a-TL  filter,  with  even  only 
two  elements  trimmed  off,  the  output  of  the  filter  is  going  to  be 
zero  at  all  times. 

To  overcome  this  limitation  in  measuring  the  transfer  function 
of  the  nonlinear  filter,  white  Gaussian  noise  was  used  as  the  input 
to  the  filter.  The  transfer  functions  were  then  computed  as 
follows.  An  input  x(n)  of  white  Gaussian  noise  was  filtered  to  get 
the  output  y(n).  The  input  and  output  sequences  were  both 
Fourier  transformed  to  get  X(eJw)  and  Y (eJw) ,  respectively.  The 
transfer  function  H(eJ“)  was  given  by 


«  !Y(e>)l  /  |X(eJ«)|.  (7) 


This  calculation  gives  the  proper  magnitude  of  the  transfer  function 
which  is  more  often  calculated  by 

Tyx  =  Pyx  /  Pxx 

where  Py*  is  the  cross  power  spectral  density  of  the  output  with 
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the  input  (  =  Y(eJ")X*(eJ“)  )  and  Pxx  is  the  auto  power  spectral 
density  of  the  input  (  =  X(eJw)X*(eJw)  ).  This  calculation  preserves 
the  phase  information  of  the  transfer  function  as  well.  However, 
we  are  only  interested  in  the  magnitude  of  the  transfer  function. 

Experimentally  it  was  found  that  by  choosing  a  window 
function  for  the  input  signal  that  minimized  leakage,  the  transfer 
function  for  the  linear  FIR  filter  calculated  using  white  noise  as  the 
input  was  equal  to  the  transfer  function  using  an  impulse  to 
within  the  resolution  of  the  computer.  That  is  to  say  the  error 
between  the  transfer  functions  calculated  using  the  two  different 
inputs,  if  any,  could  easily  be  attributed  to  round  off  error  within 
the  machine.  The  window  function  is  discussed  in  more  detail  in 
Appendix  A. 

OUTLIER  RESISTANCE 

As  was  noted  previously,  FIR  filters  inherently  perform  poorly 
in  the  presence  of  impulsive  noise  on  a  signal.  This  is  due  to  the 
effect  a  large  impulse  has  on  a  linear  filter.  Being  summed  in  each 
element  as  the  entire  window  passes  over  it  creates  a  bias  in  the 
output  signal  in  the  direction  of  the  impulse.  We  would  like  to 
minimize  the  effect  of  an  outlier  by  discarding  it  as  in  a  median 


type  filter.  If  this  were  possible  then  the  effect  of  an  impulse  on 
the  input  signal  would  become  very  small.  The  only  performance 
sacrifice  would  be,  in  effect,  a  shortening  of  the  filter  length  by 
one  element  thus  decreasing  the  Gaussian  noise  suppression 
capability  of  the  original  linear  filter.  While  this  is  a  simplification 
of  the  actual  filtering  process,  it  serves  to  illustrate  the  differences 


in  effect  an  impulse  has  on  a  linear  filter  versus  a  nonlinear  filter. 

We  can  evaluate  a  filters  performance  in  rejecting  outliers  in 
two  ways.  The  simplest  and  most  straightforward  method  is  a 
visual  measure  of  comparing  graphs  of  test  signals  to  see  which 
filter  produces  the  smoothest  output  with  impulses  present  on  the 
input  signal.  While  this  method  may  seem  imprecise  it  does  allow 
an  immediate  evaluation  of  a  filter's  performance  upon  a 
particular  input  signal  and  a  quick  comparison  of  different  filters’ 
performance.  By  comparing  the  input  to  the  output  at  any  given 
instant  it  also  gives  an  insight  into  the  filter's  physical  process. 

s 

This  gives  an  understanding  of  why  a  filter  performs  well  in  some 
situations  and  poorly  in  others  and  may  suggest  ways  to  modify  a 
particular  filter  design  to  perform  better  in  some  respect. 

The  second  method  for  evaluating  performance  would  be  to 
compute  the  error  between  the  filter  output  and  the  uncorrupted 
input  signal  to  try  to  numerically  determine  which  filter  best 
smooths  noise.  This  would  allow  a  direct  comparison  of  different 
filters’  performance  on  a  particular  set  of  input  signals.  With  a 
proper  choice  of  test  signals  we  should  be  able  to  derive  some 
insight  into  how  well  a  particular  filter  design  will  perform  under  a 
given  set  of  input  conditions.  In  order  to  find  out  what  trade-offs 
are  being  made  we  will  evaluate  the  performance  of  several  filter 
designs  on  several  different  input  signals  under  three  test 
conditions:  l)  the  signal  corrupted  by  Gaussian  noise  only;  2)  the 
signal  corrupted  by  impulsive  noise  only;  and  3)  the  signal 
corrupted  by  both  Gaussian  and  impulsive  noise.  This  will  allow  us 


to  see  if,  by  adding  a  capability  to  deal  with  outliers  on  the  input 
signal,  we  are  giving  up  too  much  Gaussian  noise  suppression. 

B.  SIMULATION  APPROACH 

There  are  several  specific  aspects  of  the  simulation  that  are 
worth  mentioning.  First,  the  equations  for  the  calculation  of  the 
transfer  functions  for  the  frequency  selectivity  tests,  given  in  part 
A  of  this  section,  do  not  allow  the  transfer  functions  to  be 
calculated  deterministically.  Since  the  inputs  for  these  calculations 
were  Gaussian  noise,  and  therefore  random,  the  transfer  functions 
were  found  probabilistically.  Five  hundred  trials  of  random  noise 
were  generated  to  calculate  each  transfer  function.  Each  of  these 
were  sent  through  the  filter  separately,  Fourier  transformed 
separately  and  then  the  output  transform  magnitude  was  divided 
by  the  input  transform  magnitude  to  generate  a  single  trial 
transfer  function  and  stored  point  by  point.  When  the  transfer 
functions  of  all  500  trials  had  been  calculated  and  stored,  the 
results  for  each  point  were  then  ordered  from  smallest  to  largest. 
The  median  at  each  point  was  chosen  to  represent  the  transfer 
function  of  the  a-TL  at  that  frequency. 

The  median  was  used  instead  of  the  mean  in  determining  the 
transfer  function  at  each  point  for  two  reasons.  First,  the 
histogram  of  values  at  each  point  usually  was  quite  skewed.  Most 
of  the  values  were  lumped  at  a  certain  distance  from  zero  (there 
could  be  no  negative  values)  with  an  exponential  distribution  fading 


away  from  the  main  lobe.  Second,  and  more  importantly,  due  to 
the  nonlinearity  of  the  filtering  there  were  individual  values  which 
were  hundreds  of  times  larger  than  most  of  the  others.  These 
individual  values  could  affect  the  output  at  a  single  point  by  as 
much  as  20%  of  the  average  value.  For  a  complete  discussion  of 
the  transfer  function  histograms  see  Appendix  A. 

Second,  in  evaluating  the  outlier  resistance  performance  of  the 
a-TL  filter  and  comparing  it  to  some  of  the  other  filters  described 
previously  it  was  important  to  develop  a  set  of  test  signals  that 
would  allow  observations  on  particular  signal  characteristics.  This 
was  basically  narrowed  to  a  comparison  of  signals  with  only 
smooth,  “slow"  variations— slow  being  relative  to  the  frequency 
band  of  the  filter— to  signals  with  sharp  edges.  A  suitable  set  of 
input  signals  was  developed  and  included  several  signals  of  very 
smooth  slow  variations,  a  sine  wave  or  slow  ramp,  for  example, 
several  signals  of  step  functions  and  a  couple  of  combination 
signals.  The  combination  signals  generally  showed  the  most 
interesting  results  and  are  used  for  most  of  the  examples  shown  in 
Chapter  V. 

Once  this  set  was  developed  each  signal  was  corrupted  with 
noise.  To  each  signal  was  added  1)  Gaussian  white  noise  only;  2) 
impulsive  (Lapladan  or  double-exponential)  white  noise  only;  and 
3)  both  Gaussian  and  impulsive  additive  noise.  The  purpose  of  these 
three  testing  schemes  was  to  try  to  show  that  the  a- trimmed 
mean  and  a- trimmed  linear  filters  performed  better  than  some  of 
the  other  models  in  impulsive  noise.  The  comparisons  would  allow 
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us  to  see  how  much  Gaussian  noise  suppression  was  being  given  up 
for  this  outlier  rejection. 


| 
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V.  RESULTS 

A.  FREQUENCY  SELECTIVITY 

As  one  would  probably  expect,  the  frequency  characteristics  of 
the  a- trimmed  linear  filter  are  similar  to  those  of  the 
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corresponding  linear  filter.  However,  as  one  might  also  expect, 
the  frequency  performance  is  not  as  good.  Naturally,  trimming 
elements  off  of  a  set  of  coefficients  will  produce  some  side  effects,  a 
degradation  of  performance.  This  is  especially  pronounced  in  the 
a-TL  filter  for  a  couple  of  reasons. 

First,  we  have  a  randomness  associated  with  the  nonlinearity 
of  this  particular  trimming  process.  The  element  chosen  to  be 
trimmed,  and  thus  the  effect  a  particular  coefficient  of  the  linear 
filter  has  on  the  output  at  a  given  time,  varies  from  point  to 
point  depending  on  the  input  data  surrounding  the  point.  The 


trimming  takes  on  a  data  dependence.  This  inherent  nonlinearity 
of  the  filter  introduces  a  randomness  that  can  only  be  expected  to 
upset  the  balance  of  a  set  of  linear  coefficients. 

Second,  although  the  trimming  is  data  dependent  and  therefore 
somewhat  random,  in  this  filter  it  is  the  largest  and  smallest  (or 
rather  the  largest  positive  and  largest  negative)  values  that  are 
trimmed.  In  the  absence  of  large  outliers  on  the  input  signal,  we 
can  expect  that  these  values  will  correspond  to  some  of  the  largest 
(in  absolute  value)  filter  coefficients.  This  means  that  some  of  the 
coefficients  that  would  normally  contribute  the  most  to  the  output 
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at  a  given  point  in  the  linear  filter  are  the  elements  that  are  being 
trimmed. 

These  two  observations  help  us  understand  some  of  the  process 
that  causes  a  degradation  of  performance.  Experimentally, 
however,  several  interesting  and  perhaps  not  so  intuitive  things 
were  discovered  about  what  effects  the  a-trimming  process  has  on 
a  set  of  linear  filter  coefficients  and  the  effects  of  certain  filter 
coefficient  parameters  on  the  trimming  process  itself.  Specifically 
there  were  three  notable  items:  1)  the  effects  of  the  filter’s 
length;  2)  the  effects  of  the  designed  linear  filter’s  rolloff 
characteristics;  3)  the  symmetry  properties  maintained  through 
the  trimming;  and  4)  the  DC  gain  problem  the  a-trimming  process 
introduces.  Each  of  these  topics  will  be  discussed  in  some  detail. 
EFFECTS  OF  FILTER  LENGTH 

As  we  know  from  linear  filter  design  class,  the  more  elements 
we  have  in  a  filter  the  better  the  frequency  response  we  may 
obtain.  More  elements  in  the  filter  mean  we  can  design  a  better 
attenuation  in  the  stop  band,  a  steeper  rolloff  in  the  transition 
band,  a  flatter  response  in  the  pass  band  or  more  often  some 
combination  of  these  design  parameters.  Thus,  when  we  start 
trimming  away  some  of  the  coefficients  of  a  filter  we  can  only 
expect  that  some  or  all  of  these  basic  parameters  may  suffer.  In 
addition,  if  the  trimming  is  done  in  a  random  or  nonlinear  fashion 
we  can  expect  a  certain  amount  of  randomness  or  nonlinearity  to 
be  introduced  into  these  filter  parameters  as  well. 

During  a  normal  filter  design,  one  of  the  major  trade-offs 


encountered  is  the  cost  of  increased  complexity  caused  by  a  longer 
filter  with  the  desired  stop  band  attenuation  and  transition  band 
roUoff.  The  importance  of  the  flatness  of  the  pass  band  and  stop 
band  are  normally  accounted  for  by  the  choice  of  the  filter 
design— Butterworth,  Chebyshev,  Elliptical,  etc.— and  to  a  smaller 
extent  by  the  length  of  the  filter.  Usually  the  driving  factor  in  a 
filter  design  will  be  that  the  attenuation  in  the  stop  band  meet 
some  minimum  requirement.  This  is  balanced  against  the  cost 
(complexity/length)  of  the  filter. 

An  interesting  effect  occurs,  however,  when  we  start  talking 
about  an  a-trimmed  linear  filter  with  a  set  of  FIR  coefficients.  An 
attenuation  limit  is  reached  quite  rapidly  even  though  the  filter  is 
made  longer!  By  comparing  the  graphs  of  BPF2500s  for  the  lengths 
of  31  and  63  elements.  Figures  3  and  4,  and  the  graphs  of 
HPF3500  for  these  same  lengths,  Figures  5  and  6,  one  can  see  that 
the  attenuation  does  not  get  any  better  for  the  longer  filter.  In 
fact,  if  the  graphs  were  superimposed  it  could  be  seen  that  the 
responses  are  nearly  identical.  This  is  especially  true  for  the  T=2 
graphs.  The  attenuation  is  the  same  for  all  frequencies  with  only 
some  minor  discrepancies  in  the  rolloff  region.  As  T  increases  this 
variation  in  the  transition  band  increases  somewhat.  The  rolloff  is 
slightly  steeper  for  the  longer  filters  as  would  be  expected.  In 
addition,  the  first  hump  is  slightly  narrower,  a  little  higher,  and 
rises  more  sharply  for  the  longer  filter.  These  facts  all  seem  to  be 
in  line  with  what  we  would  expect  for  a  larger  linear  filter. 


In  order  to  complete  the  comparisons,  note  the  maximum 
attenuations  for  each  of  the  filters.  For  the  BPF  the  maximum 
attenuation  of  the  trimmed  filters  near  0  and  .5  normalized 
frequency  is  about  -32  to  -33dB.  This  is  for  both  the  31  element 
filter  and  the  63  element  filter.  Now  compare  this  to  the 
attenuation  obtained  by  the  actual  linear  filters.  Figure  7  shows 
that  the  31  element  linear  BPF  has  an  attenuation  of  -47dB  while 
Figure  8  shows  that  the  63  element  filter  has  an  attenuation  of 
-84dB,  an  additional  37dB  of  attenuation!  Yet  the  trimmed  filters 
have  nearly  identical  attenuation.  A  similar  case  exists  for  the  HPF 
example.  The  linear  filters  of  length  31  and  63  have  stop  band 
attenuation  figures  of  -56dB  and  -102dB,  respectively.  A  difference 
of  46dB!  Yet  the  maximum  attenuation  of  either  of  these  filters 
when  trimmed  is  only  about  -28dB.  This  limiting  effect  was  seen 
in  all  of  the  examples  in  the  simulation. 

This  effect  is  probably  most  closely  related  to  one  of  the  effects 
discussed  at  the  beginning  of  this  section.  In  the  a-trimming 
process,  the  data  elements  that  most  often  get  trimmed  are  those 
that  are  modified  by  the  largest  linear  coefficients.  These 
coefficients  are  responsible  for  a  large  portion  of  the  filter’s  overall 
shape.  When  the  length  of  the  filter  is  increased,  the  additional 
elements  are  typically  very  small  and  have  a  very  fine  effect. 

They  only  become  significant  when  combined  with  the  larger 
coefficients  of  the  filter.  Taking  out  the  larger  coefficients  probably 
nullifies  most  of  the  subtle  effect  the  smaller  coefficients  have  in 
lowering  the  stop  band. 
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EFFECTS  OF  DESIGNED  ROLLOFF 

If  you  were  watching  closely  in  the  previous  section  you  may 
have  noticed  something  strange  that  was  not  emphasized.  In  the 
filter  examples  given  the  linear  band  pass  filters  had  attenuations 
of  -47dB  and  -84dB  for  31  and  63  elements,  respectively.  Likewise 
the  high  pass  filters  had  corresponding  attenuations  of  -56dB  and 
-102dB  in  the  stop  band.  Yet,  when  2  elements  were  trimmed 
from  the  BPF  the  attenuation  only  dropped  to  *-33dB,  while  the 
same  amount  of  trimming  off  the  HPF  caused  the  attenuation  to 
rise  all  the  way  to  «-28dB.  The  BPF  suffered  less  degradation  than 
the  HPF!  This  seems  counter-intuitive. 

Taking  a  closer  look  by  comparing  the  two  31  element  filters 
we  see  that  the  attenuation  of  the  BPF  dropped  from  -47dB  to 
roughly  -33dB  when  two  elements  were  trimmed.  While  the  HPF, 
whose  better  original  attenuation  of  -56dB,  rose  to  about  -28dB 
when  the  same  trimming  was  done.  Obviously  there  is  some  effect 
other  than  filter  length  involved  here. 

It  turns  out  that  the  answer  is  in  the  coefficients  themselves, 
or  at  least  in  their  design.  The  main  difference  between  the  high 
pass  and  the  band  pass  filter  of  the  previous  section  is  the  width  of 
the  transition  band.  The  BPF  was  designed  with  a  steeper  rolloff 
than  the  HPF.  It  may  seem  improbable  that  this  is  indeed  the 
reason  for  the  discrepancy  in  the  stop  band  attenuation  since  the 
filters  pass  different  frequency  bands,  so  the  following  example  is 
offered. 

A  second  band  pass  filter  was  designed  with  the  same,  more 
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gentle  rolloff  as  the  HPF  recently  discussed.  The  frequency  response 
of  the  31  element  filter  is  shown  in  Figure  9.  Note  that  the  linear 
filter  has  a  -56dB  attenuation  in  the  stop  band,  the  same  as  the 
HPF  of  the  previous  section.  In  fact,  in  other  filter  designs  it  did 
not  matter  to  where  the  actual  pass  band  was  moved,  if  the 
rolloff  was  kept  at  a  width  of  0.1  (normalized  frequency)  and  the 
pass  band  was  kept  reasonably  wide,  then  the  attenuation  in  the 
stop  band  was  very  nearly  -56dB  for  all  examples.  Now  note  in 
Figure  10  that  the  attenuation  in  the  stop  band  does  not  go  below 
about  -22dB  when  2  elements  are  trimmed.  Yet  this  filter  passes 
nearly  the  same  frequency  band  as  the  previously  discussed  BPF 
whose  attenuation  went  from  -47dB  to  -33dB  when  2  elements 
were  trimmed. 

It  appears  then  that  the  coefficients  for  an  FIR  linear  filter  are 
more  resistant  to  a- trimming  when  the  filter  is  designed  with  a 
steeper  rolloff,  or  at  least  the  effect  of  a-trimming  is  dependent  on 
the  characteristics  of  the  coefficients  in  some  way.  Somehow  in 
these  two  examples  the  coefficients  of  the  filter  with  the  steep 
rolloff  were  more  resistant  to  a-trimming  than  the  other  two 
shallower  filters.  (We  can  say  this  for  a-trimming  in  general 
because,  although  it  is  hard  to  determine  what  is  a  good  measure 
of  the  attenuation  in  the  stop  band  for  trimming  greater  than  two 
elements,  it  is  visually  obvious  that  the  attenuation  for  T=4  and 
T=6  was  better  for  the  steep  rolloff  BPF  than  for  either  of  the 
other  two  filters  to  which  it  was  compared.) 

However,  the  previous  two  examples  both  used  the  same 
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Figure  10  :  ot-TL  BPF2500  L31  Transfer  Functions 
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example  of  a  steep  rolloff  filter  for  comparison.  It  could  be  just  a 
Alike  of  the  particular  coefficients  of  that  filter  that  the 
a- trimming  happened  to  have  less  of  an  effect  on  it.  So  another 
steep  rolloff  filter  was  designed. 

This  time  it  was  a  LPF  having  a  rolloff  similar  to  the  steep 
BPF,  but  a  much  narrower  pass  band.  In  this  case  the  rolloff  was 
0.04  normalized  frequency  with  a  pass  band  of  a  mere  0.01.  The 
designed  attenuation  in  the  stop  band  was  -28dB  for  the  31 
element  filter  and  -52dB  for  the  63  element  filter.  The  results  for 
this  filter  were  even  better  than  for  the  steep  BPF  as  can  be  seen 
in  Figures  11  and  12.  The  deviation  for  the  31  element  filter  is 
only  1  or  2dB  for  each  successive  pair  of  elements  trimmed  and 
the  deviation  for  the  63  element  filter  reaches  the  same  levels  only 
taking  a  part  of  the  stop  band  to  reach  that  point. 

SYMMETRY  PROPERTIES 


A  somewhat  unexpected  property  of  the  a- trimming  process 
was  the  high  degree  of  symmetry  it  maintained.  This  symmetry 
came  in  two  important  forms.  First,  the  graphs  of  each  of  the 
band  pass  filters  in  the  previous  section  showed  that  the  transfer 
functions  maintained  the  original  band  pass  filter's  symmetry 
about  the  pass  band.  The  attenuation  was  symmetric  in  each  of 
the  two  stop  bands. 

This  situation  was  not  necessarily  expected.  It's  existence  shows 
that  the  a-trimming  process  does  not  have  any  frequency 
selectivity.  The  trimming  in  the  a-TL  filter  does  not  modify  the 
set  of  linear  filter  coefficients  as  does  the  ct-trimmed  mean  filter. 


The  a-TM  filter’s  coefficients  are  modified  depending  on  how  many 
elements  are  trimmed  from  the  filter.  No  such  modification  occurs 
in  the  a-TL  filter.  Thus,  it  could  easily  be  expected  that 
a- trimming  a  linear  filter  may  affect  its  frequency  properties.  The 
simulation  shows  that  this  is  not  the  case,  however. 

Second,  the  other  important  symmetry  property  maintained 
was  one  dealing  with  similar  filters.  The  LPF  and  HPF  transfer 
functions  shown  in  Figures  13  and  14,  respectively,  show  how  the 
a- trimming  process  maintains  symmetry  across  the  frequency 
spectrum.  As  one  can  see  in  Table  1,  the  coefficients  of  each  of 
the  filters  have  the  same  magnitudes.  They  only  differ  in  sign. 

One  might  expect  that  the  the  filters  would  be  trimmed 
differently,  since  the  largest  positive  and  negative  elements  are 
trimmed  in  the  a- trimming  process,  causing  different  transfer 
functions. 

Table  1:  LPF1500  L9  and  HPF1500  L9  Coefficients 

LEE15QQ-L9  HEE25QQ.JL? 

-0.04739  =  h(l)  *  h(9)  =  -0.04739 

0.05110  =  h(2)  =  h(8)  =  -0.05110 

0.15307  =  h(3)  *  h(7)  =  0.15307 

0.25926  =  h(4)  =  h(6)  =  -0.25926 

0.30437  =  h(5)  =  0.30437 

The  explanation  for  this  symmetry  seems  fairly 
straightforward.  The  noise  used  to  calculate  the  transfer  functions 
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Figure  13  :  a-TL  LPF1500  LI  5  Transfer  Functions  plus  Ideal 

Frequency  Response 


is  zero  mean  and  Gaussian  and  is,  therefore,  equally  likely  to  be 
negative  as  positive.  Thus,  on  the  average,  the  effects  of  a 
difference  in  signs  of  the  coefficients  will  not  be  a  factor  in  the 
long  run  and  only  be  a  factor  in  determining  the  frequencies 
passed  in  any  given  interval,  high  frequencies  or  low  frequencies 
depending  on  the  filter.  An  interesting  test  would  be  to  give  the 
noise  a  dc  bias  so  that  every  data  point  was  the  same  sign.  The 
trimming  symmetry  that  these  two  filters  exhibit  may  suffer 
under  these  input  conditions. 

FILTER  GAIN 

A  troubling  problem  with  the  a-TL  filter  is  the  uncontrolled 
gain  of  the  system.  The  gain  calculation  outlined  in  Section  III  does 
a  good  job  of  estimating  the  gain  of  an  a-TL  filter  only  if  the 
linear  filter  is  a  low  pass  filter.  Experimentally,  the  gain  factor 
correction  was  very  good  for  the  LPF  case.  However,  it  was  totally 
unpredictable  for  any  filter  not  passing  dc.  This  gain  problem 
showed  up  in  the  transfer  functions  as  a  spike  at  0  Hz.  (  This  spike 
was  not  evident  when  the  median  was  used  to  select  the  transfer 
function,  only  the  mean.  The  distribution  of  values  at  the  dc  point 
has  a  much  larger  deviation  than  at  any  other  frequency  for  all  of 
the  filters  in  the  simulation.  See  Appendix  A  for  more  details.) 

Interestingly  enough,  though,  an  unsuspected  consistency  was 
found  while  attempting  to  calculate  the  error  to  evaluate  the 
performance  of  an  a-trimmed  BPF  in  rejecting  outliers.  It  was 
necessary  too  correct  the  gain  problem  in  order  to  calculate  the 
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MSE  between  an  ideal  signal  and  the  noise  corrupted  signal  filtered 
by  the  a-TL  filter.  So  the  output  was  modified  by  a  minimum 
MSE  calculation.  This  minimization  led  to  the  calculation  of  the 
required  normalization  factor.  Surprisingly,  the  normalization 
factor  (same  as  the  gain  factor  g)  for  a  given  a-TL  filter  remained 
fairly  constant  regardless  of  the  noise  on  the  input  signal.  The 
normalization  factors  for  two  of  the  band  pass  filters  are  given  in 
Figures  15  and  16.  The  values  in  these  figures  were  obtained  by 
filtering  the  ideal  signal— a  sinusoid  in  the  pass  band  of  the 
filter— with  the  a-TL  BPF  and  then  calculating  the  required 
normalization  factor  to  minimize  the  MSE  between  the  filter  output 
and  ideal  input. 

The  normalization  factors  calculated  when  there  was  noise  on 
the  input  signals  were  consistently  smaller  than  the  values  in  the 
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Figure  15  :  Gain  Factor  for  a-TL  BPF2500s 
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Figure  16  :  Gain  Factor  for  a-TL  BPF2500 


figures,  but  right  in  the  same  neighborhood.  This  would  be 
expected  since  the  noise  on  the  signal  would  almost  certainly  add 
slightly  to  the  MSE  between  the  ideal  and  the  corrupted  signal. 
OVERALL  FREQUENCY  PERFORMANCE 

One  of  the  items  that  has  not  been  stressed  very  much  is  the 
actual  measurement  of  the  frequency  performance  in  the  stop 
band  of  the  a-TL  filters.  Where  a  reference  has  been  made  to  a 
value  for  attenuation  in  the  stop  band  of  an  a-TL  filter,  the 
measurement  was  made  at  the  edge  of  the  filter  response  where 
the  attenuation  had  reached  its  best  point.  An  important 
measurement,  though,  is  how  good  the  attenuation  is  at  the  first 
siddobe  For  the  original  set  of  linear  FIR  filters  these  two 
measurements  were  the  same  since  the  attenuation  was  uniform 
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across  the  stop  band. 

This  is  far  from  the  case  when  measuring  the  attenuation  in 
the  first  sidelobe  of  an  a-TL  filter.  As  can  be  seen  in  Figure  17 
this  performance  measure  leads  to  some  very  poor  results. 
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Figure  17  :  First  Sidelobe  Level  of  BPF2500s 


This  figure  shows  that  the  most  serious  loss  of  attenuation 
occurs  for  the  first  two  elements  trimmed.  After  that  the  loss 
appears  to  be  fairly  linear  and  not  nearly  so  severe.  The  first 
sidelobe  is  most  affected  in  the  long  filter  lengths.  Since  the 
sidelobes  become  more  numerous  and  skinnier,  the  loss  of 
attenuation  shows  up  more  clearly  for  the  very  narrow  first 
sidelobes  of  the  long  filters. 
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This  makes  it  difficult  to  say  exactly  how  bad  the  degradation 
of  frequency  performance  is  for  the  a-TL  filter.  The  attenuation  is 
-  much  better  than  the  level  at  this  first  sidelobe  for  the  majority 
of  the  stop  band.  Yet  a  significant  amount  of  energy  is  passed 
through  the  filter  in  this  first  sidelobe  that  would  be  attenuated  in 
the  corresponding  linear  filter.  About  all  that  can  be  said  for  sure 
is  that  the  frequency  response  of  the  linear  filter  suffers  quite  a 
bit  when  the  filter  is  a- trimmed.  The  longer  the  linear  filter  and 
better  the  stop  band  attenuation,  the  more  a- trimming  degrades 
the  performance. 

B.  OUTLIER  RESISTANCE 


The  property  of  outlier  resistance  was  one  of  the  major  design 
considerations  in  all  of  the  filter  designs  discussed  in  Section  11.  One 
would  therefore  expect,  for  the  low  pass  filter  case  embodied  in  all 
of  those  designs,  that  adequate  noise  smoothing  with  impulse 
rejection  has  been  achieved  with  some  degree  of  success.  It  seems 
unlikely  that  a  new  filter  design  which  has  outlier  resistance  as 
only  half  of  the  major  filter  objectives— the  a-TL  filter— would  be 
able  to  compete  with  any  of  the  previously  mentioned  filters.  On 
the  other  hand,  there  is  nothing  to  compare  a  band  pass  or  high 
pass  a-TL  filter  to  except  the  corresponding  linear  filter.  In  this 
case  one  would  hope  that  the  a-TL  filter  could  outperform  the 
linear  filter  in  rejecting  impulses. 


VISUAL  EVALUATION 


In  visually  inspecting  the  performance  of  the  a-TL  filter, 
however,  it  is  not  practical  to  look  at  a  high  pass  or  even  band 
pass  signal.  No  meaningful  information  is  readily  derived  from 
trying  to  do  so.  In  order  to  determine  the  effectiveness  of  the  a-TL 
filter  in  rejecting  impulses  by  a  visual  inspection  it  is  necessary  to 
limit  the  investigation  to  the  LPF  case.  This  approach  also  allows  a 
comparison  of  the  a-TL  filter  to  several  of  the  other  filter  designs 
discussed  in  Section  II.  Specifically,  a  comparison  of  the  low  pass 
a-TL  filter  will  be  made  to  the  a-TM  filter  (and  the  mean  and 
median  filters,  noting  that  these  are  special  cases  of  the  a-TM 
filter  having  0  and  N-l  elements  trimmed,  respectively) . 

The  input  signal  that  best  showed  the  filters’  performance  was 
a  signal  similar  to  the  one  used  by  Lee  and  Kassam  in  one  of  their 
papers  on  order  statistic  filters.1  The  signal  has  several  edges  and 
several  monotone  regions.  The  signal  is  shown  in  Figure  18,  before 
and  after  it  is  corrupted  by  white  Gaussian  noise,  a2  =  0.2.  Figure 
19  shows  the  same  signal  corrupted  by  white  Laplacian  noise, 
o2  =  0.5.  This  signal  is  known  as  signal  1  (SIG1)  and  the  two 
corrupted  signals  are  known  as  S1G1G  (Gaussian)  and  SIGH 
(impulsive).  A  third  version  of  this  signal,  SIG11G,  is  SIGH  plus 
white  Gaussian  noise,  o2  =  0.1,  and  is  shown  in  Figure  20. 

The  filters  used  in  the  simulation  for  the  visual  evaluation  were 
LPF300  L9  and  LPF1500  L9.  The  9  element  filters  were  chosen 
because  they  were  short  compared  to  the  features  of  the  signal. 

The  LPF300  signed  is  particuleurly  interesting  because  it  was 
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Figure  19  :  SIG1I  -  signal  1  +  noise  ~L(0,0.5) 
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Figure  20  :  SIG1IG  -  SIGH  +  noise  ~N(0,0.1) 


designed  to  have  a  transfer  function  very  similar  to  that  of  a 
mean  filter,  having  a  similar  pass  band  width  and  rolloff. 

The  results  of  SIG1G  being  filtered  by  a  mean  and  median  filter 
of  length  9  and  a-TM  9T6  are  shown  in  Figures  2i,  22  and  23, 
respectively.  As  expected,  the  a-TM  is  a  very  good  compromise 
between  the  noise  smoothing  of  the  mean  filter  and  the  edge 
preservation  of  the  median  filter.  Figure  24  shows  that  the  LPF300 
L9  filter,  even  though  it  has  a  similar  transfer  function  to  the 
mean  filter,  does  not  do  as  good  a  job  smoothing  the  Gaussian 
noise,  but  still  introduces  the  same  smearing  of  edges.  This  is  due 
to  the  nearly  constant  coefficients  of  the  filter.  However,  if  the 
LPF300  L9  filter  is  a-trimmed  6  elements,  the  results  are  very 
comparable  to  the  a-TM  filter  of  the  same  length  and  trimming. 
The  a-TL  filter  fails  to  smooth  the  noise  as  well  near  the  end  of 
the  sinusoidal  section,  but  matches  the  flat  pulse  much  better  as 
shown  in  Figure  25.  The  LPF1500  L9  filter  lets  through  too  much 
of  the  noise,  even  when  a-trimmed  6  elements  as  shown  in 
Figures  26  and  27,  respectively. 

Of  a  more  primary  interest  is  how  these  filters  perform  in  the 
presence  of  impulsive  noise.  The  mean,  a-TM  and  median  filters 
perform  as  expected  in  Figures  28,  29  and  30,  respectively.  The 
important  spots  to  watch  in  all  of  these  figures  are  the  doublet 
type  impulses  in  the  first  flat  region,  the  single  and  double  width 
impulses  at  the  top  of  the  sinusoid,  the  pair  of  opposite  pulses  in 
the  middle  of  the  sinusoid  region,  the  doublet  type  impulses  at  the 
falling  edge  of  the  step  and  the  string  of  five  impulses  on  the  final 
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Figure  21  :  SJG1G  Filtered  by  Mean  9 
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Figure  22  :  SI616  Filtered  by  Median  9 
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Figure  23  :  SI61G  Filtered  by  a-TM  9T6 
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Figure  25  :  S1G1G  Filtered  by  a-TL  LPF300  9T6 
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Figure  28  :  SIGH  Filtered  by  Mean  9 
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Figure  29  :  SIG1I  Filtered  by  a-TM  9T6 
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Figure  30  :  SIGlf  Filtered  by  Median  9 


ramp.  Note  in  each  of  the  three  preceding  figures,  the  filters  dealt 
with  all  of  the  impulses  fairly  well  except  the  impulses  at  the  top 
of  the  sinusoid  and  those  on  the  final  ramp. 

Amazingly,  the  LPF300  L9  filter  takes  care  of  the  effects  of  all 
of  the  impulses  quite  well  as  shown  in  Figure  31.  However,  this 
filter  still  lets  through  too  much  small  level  noise  and  introduces 
quite  a  bit  of  distortion  to  the  signal.  It  can  be  seen  in  Figure  32 
that  when  this  filter  is  a- trimmed  6  elements  the  distortion  and 
low  level  noise  are  taken  care  of,  but  the  pair  of  impulses  on  the 
sinusoid  have  quite  an  effect.  The  real  power  of  a- trimming  is 
seen  on  the  LPF1500  L9  In  Figures  33  and  34  one  can  sec 
that  this  linear  filter  still  passes  too  much  high  frequency  noise, 
but  the  distortion  is  smaller  than  that  of  LPF300  due  to  the  larger 
pass  band  being  better  able  to  follow  the  edges  in  the  signal.  Also, 
trimming  six  elements  has  some  degree  of  success  against  the 
outliers.  In  fact,  this  is  the  only  filter  that  can  be  said  to  have 
taken  out  the  majority  of  the  effect  of  the  pair  of  impulses  on  the 
top  of  the  sinusoid. 

Although  these  graphs  give  a  feeling  of  what  can  be  expected 
from  an  a-TL  filter  compared  to  one  of  the  more  traditional 
filters,  it  is  still  difficult  to  sec  what  is  actually  happening  in  the 
a-TL  filter  One  example  that  delivered  some  really  tangible  visual 
results  is  shown  in  Figure  35  In  this  graph  the  input  signal  is 
S1G11G.  The  two  different  outputs  plotted  here  are  those  for  a 
median  filter  of  length  9  and  LPF300  L9  a- trimmed  eight 
elements,  in  effect  an  a-TL  median  filter  At  first  one  may  think 
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Figure  31  :  SIGH  Filtered  by  LPF300  L9 
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Figure  52  :  SIG1I  Filtered  by  o-TL  LPF300  9T6 
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Figure  33  :  SIGH  Filtered  by  LPF1500  L9 
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Figure  35  :  Results  of  SIG1I6  Filtered  by  Median  9  and  a-TL 

LPF300  9T8 


that  these  two  filters  should  produce  the  same  output  since  they 
are  both  nine  element  filters  with  eight  elements  trimmed. 

However,  one  must  remember  that  the  a-TL  filter  weights  the 
input  by  the  set  of  linear  coefficients  before  the  values  are  ordered 
and  trimmed.  So  even  though  an  element  in  the  input  window  of 
an  cr-TL  filter  is  the  largest  it  is  not  necessarily  going  to  be 
trimmed.  It  is  only  trimmed  if  its  value  times  its  weighting 
coefficient  are  large  enough  to  force  it  to  the  ends  of  the  window 
after  sorting.  As  can  be  seen  in  the  figure,  the  two  graphs  give 
the  same  outputs  in  some  areas  and  arc  significantly  different  in 
others.  The  most  notable  feature  is  that  the  two  filters  produce 
nearly  identical  outputs  in  all  areas  of  the  graph  where  the  signal 
is  changing  most  rapidly:  near  each  of  the  edges  and  along  the 
center  of  the  sinusoid.  The  areas  where  the  outputs  deviate  the 
most  are,  in  general,  the  flattest  parts  of  the  signal. 

This  suggests  that  something  could  possibly  be  done  to  improve 
the  performance  of  the  a-TL  filter  or  at  least  suggests  what  the 
real  differences  are  between  the  a-TL  and  a-TM  filters.  The  reason 
for  the  poor  performance  in  the  flat  areas  of  the  signal  is  most 
likely  due  to  the  fact  that  the  linear  filter  in  this  case  has  a  finite 
band  pass.  We  know  that  the  mean  filter  is  optimum  in  smoothing 
Gaussian  noise,  but  the  linear  LPF  must  pass  a  certain  band  of 
frequencies.  Note  that  the  oscillations  in  these  flat  areas  are  not 
nearly  so  great  as  they  were  in  the  original  signal.  The  high 
frequency  noise  has  been  filtered  out,  but  some  low  frequency 
noise  remains.  The  steep  areas  of  the  plot  are  changing  fast  enough 
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that  they  are  at  the  edge  of  the  filter  band,  at  least,  thus  the 
a-TL  filter  can  do  a  good  job  of  smoothing  noise  that  has  any 
higher  frequency  components  in  these  relatively  steep  areas. 

NMSE  COMPARISON 

A  more  effective  or  at  least  a  more  informative  method  of 
evaluating  the  performance  of  the  a-TL  filter  in  resisting  outliers  is 
with  an  error  analysis  between  the  filtered  output  signals  and  the 
original  "ideal*  signals  This  discussion  will  concentrate  on  three 
different  ideal  signals  each  corrupted  by  either  Gaussian  noise, 
Lapladan  (impulsive)  noise,  or  an  additive  combination  of  both 
types  of  noise.  The  first  signal  is  a  low  frequency  signal  with  edges 
in  it.  It  is  the  signal  described  above  as  S1G1  Following  the  same 
notation,  a  second  signal,  S1G2,  is  a  sinusoidal  signal  with  a 
frequency  in  the  center  of  the  pass  band  for  two  different  band 
pass  signals.  S1G3  is  a  similar  signal  with  a  frequency  that  will 
place  it  in  the  pass  band  of  both  of  the  LP  fitters  previously 
discussed.  These  three  signals  will  allow  several  different 
comparisons  and  evaluations  Signals  1  and  3  will  allow  different 
filters  to  be  compared  in  their  performance  with  respect  to  signals 
that  have  edges  versus  slowly  varying  signals  These  two  signals 
also  give  two  examples  for  comparing  the  a-TL  filter  to  the  a-TM 
filter.  Finally,  signals  2  and  3  will  show  the  relationship  between 
the  low  pass  and  band  pass  cases  of  the  a-TL  filter 

For  every  case  tested  two  errors  were  measured  normalized 
mean  square  error  (NMSE)  and  normalized  average  error  (NAE) 
These  errors  were  calculated  in  the  following  manner 


n  =  1,2,.  .  .N  (9) 


NMSE  =  (1/N)  I  [  |y(n)|  -  H(n)|  ]2, 

2  t(n)P 

where  y(n)  is  the  filter  output  and  i(n)  is  the  ideal  signal  and 

NAE  =  (1/N)  I  |y(n)  -  i(n)|,  n  =  1,2 . N.  (10) 

2  |i(n)| 


Experimentally,  as  evidenced  by  Figures  36  and  37,  NAE 
followed  the  same  pattern  as  NMSE  almost  without  exception.  For 
this  reason  NMSE  will  be  the  only  error  discussed  since  it  is  the 
more  common  measure  of  error. 

As  evidenced  by  these  two  graphs,  the  a-TM  filter  performs 
better  and  better  against  Gaussian  noise  on  a  signal  with  edges  in 
it  as  the  number  of  elements  trimmed  is  increased.  This  is  largely 
due  to  the  reduction  of  distortion  on  the  edges  by  the  reduced 
number  of  elements  in  the  averaging.  This  effect  is  so  acute  that 
the  median  filter  has  the  minimum  NMSE  for  each  filter  length. 

This  conclusion  is  backed  up  by  looking  at  the  NMSE  plot  of  this 
tame  set  of  filters  when  the  input  was  the  smooth  signal  3 
corruf^ed  by  Gaussian  noise  as  shown  in  Figure  38. 

Note  in  this  graph  that  the  vertical  axis  is  only  a  small  portion 
of  that  of  the  previous  one.  Note  also  that  trimming  has  no 
advent  again  n  effect  In  fact,  the  error  between  filter  output  and 
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Figure  36  :  SIG1G  Filtered  by  a-TM  -  NAE 
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Figure  37  :  S1G1G  Filtered  by  a-TM  -  NMSE 


the  ideal  signal  actually  increases  somewhat  due  to  the  decreased 
number  of  dements  used  to  smooth  the  noise  and  the  inferior 
ability  of  the  median  filter  to  smooth  Gaussian  noise. 
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Figure  38  :  SIG3G  Filtered  by  o-TM  -  NMSE 

The  questions  to  ask  now  are  1)  can  the  a-TL  filters  perform 
as  well  in  smoothing  Gaussian  noise  as  the  a-TM,  and  2)  how  do 
the  two  compare  in  the  impulsive  noise  case?  Figures  39  and  40 
show  some  interesting  results.  In  the  case  of  SIG1G  the  LPF300 
a-TL  filters  (hollow  point  markers)  perform  like  the  a-TM  filters 
shown  in  Figures  36  and  37.  The  distortion  introduced  by  the 
linear  filters  by  themselves  is  gradually  reduced  by  more  and  more 
trimming.  Whereas  the  LPF1500  filters  (solid  point  markers)  are 


able  to  follow  the  edges  of  the  signal  1  fairly  well  initially  due  to 
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Figure  39  :  SIG1G  Filtered  by  a-TL  -  NMSE 
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Figure  40  :  SIG3G  Filtered  by  a-TL  -  NMSE 


the  larger  pass  band,  but  trimming  causes  a  degradation  in  noise 
smoothing  ability  and  thus  an  increase  in  the  NMSE.  For  the 
smooth  S1G3G  the  results  are  turned  in  favor  of  the  LPF300  filters’ 
ability  to  block  out  more  noise  than  the  wider  pass  band  of  the 
LPF1500  filters.  Indeed,  the  LPF300  filters  perform  nearly  as  well 
as  the  a-TM  filters.  It  is  interesting  to  note  the  unexpected 
decrease  in  error  for  a-TL  LPF300  L9  in  going  from  T=2  to  T=4 
and  to  a  lesser  extent  the  same  thing  happening  to  LPF300  L15  in 
Figure  40.  Perhaps  some  sort  of  resonance  is  touched  in  the  filter 
to  cause  this. 

The  a-TM  filter  and  the  a-TL  filter  perform  nearly  in  the  same 
manner  in  the  presence  of  impulsive  noise  on  signal  1  as  in  the 
presence  of  Gaussian  noise.  This  is  most  likely  due  to  the  number 
of  edges  in  signal  1.  The  filters  that  can  best  deal  with  the  edges 
have  the  lowest  NMSE  values.  The  same  characteristics  which 
allow  these  filters  to  handle  the  edges  in  the  signal  are  also  the 
same  characteristics  that  help  the  filters  deal  with  outliers.  An 
interesting  point  to  note  here  is  that  LPF300  filters  have  the 
desired  performance  benefits  derived  from  trimming  elements  off  of 
the  filter  (NMSE  goes  down  as  T  goes  up),  while  the  LPF1500 
filters  have  a  better  initial  starting  point  in  both  the  impulsive  and 
Gaussian  noise  cases.  It  would  be  interesting  to  design  a  filter  with 
the  steep  rolloff  of  LPF300  having  the  pass  band  of  LPF1500  to  see 
if  a  better  performing  a-TL  filter  could  be  obtained. 

The  results  for  the  performance  of  all  of  the  filters  on  signal  3 
in  the  presence  of  impulsive  noise  are  somewhat  surprising.  The 


initial  NMSE  between  S1631  and  the  ideal  signal  before  filtering  is 
0.504,  while  for  S1G3G  the  NMSE  was  only  0  375  Yet,  nearly  all 
of  the  filters  more  effectively  suppressed  the  impulsive  noise  than 
the  Gaussian  noise  as  can  be  seen  when  Figures  41  and  42  are 
compared  to  the  figures  previously  shown  for  S1G1G.  Note  the 
different  scales  on  the  vertical  axes.  Notice  also  in  this  case  that, 
in  the  presence  of  impulsive  noise,  trimming  only  has  a  beneficial 
effect  for  filters  of  short  length  and  that  the  median  is  not  the 
best  filter  to  use  on  this  smooth  signal. 

The  error  calculations  for  all  of  the  filters  were  predictable  for 
the  case  of  the  two  signals  corrupted  by  Gaussian  and  impulsive 
noise  at  the  same  time.  Superposition  of  performance  in  the  two 
cases  seemed  to  hold  very  well. 


alpha-TM  on  Signal  31  - 
NMSE  -  0.504 


Number  Trimmed 

Figure  41  :  SIG3I  Filtered  by  ct-TM  -  NMSE 
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Figure  42  :  SIG3I  Filtered  by  o-TL  -  NMSE 

The  really  interesting  example  in  these  error  calculations  is  the 
case  of  signal  2  and  the  band  pass  filters.  As  was  pointed  out 
earlier,  the  ideal  signal  is  just  a  sinusoid  whose  frequency  lies  in 
the  pass  band  of  the  two  band  pass  filters  discussed  in  part  A.  As 
one  can  see  in  Figure  43  the  linear  filter  performs  an  adequate  job 
of  smoothing  Gaussian  noise.  Trimming  elements  off  of  either  filter 
docs  not  have  a  very  big  or  very  consistent  effect.  In  general, 
a- trimming  tends  to  help  the  shorter  filters  and  harm  the  longer 
ones.  This  is  probably  due  to  the  trimming  introducing  more 
distortion  in  the  longer  filters  since  the  signal  changes  are  short 
compared  to  the  lengths  of  the  longer  filters. 

However,  when  the  signal  is  corrupted  with  impulsive  noise, 
a- trimming  can  have  a  great  effect  as  evidenced  in  Figure  44.  For 
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this  signal,  the  kind  of  affect  the  a-TL  filter  was  supposed  to  have 
is  finally  realized.  In  this  instance,  the  choice  of  an  a-trimmed 
linear  filter  clearly  has  an  advantage  over  its  corresponding  linear 
filter,  especially  for  only  two  or  four  elements  being  trimmed. 

As  one  might  expect,  when  signal  2  is  corrupted  by  both  Gaussian 
and  impulsive  noise  it  maintains  the  basic  performance  shown  in 
the  previous  figure.  The  a- trimming  helps  reject  the  impulses  and 
thereby  decreases  the  NMSE  without  sacrificing  very  much  in  the 
way  of  Gaussian  noise  suppression. 
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VI.  SUMMARY,  FUTURE  RESEARCH  AND  CONCLUSIONS 


A.  SUMMARY 

The  a- trimmed  linear  filter  model  came  out  of  a  natural 
development  of  evolutionary  filter  designs.  It  is  another  step  in  the 
process  of  trying  to  marry  desirable  traits  from  different  families 
of  filters,  if  is  noteworthy  in  that  it  is  a  first  attempt  at 
combining  the  inherently  nonlinear  filter  characteristic  of  outlier 
resistance  with  the  definitely  linear  filter  characteristic  of 
frequency  selectivity.  A  viable  filter  design  having  these  two 
primary  characteristics  could  be  extremely  useful. 

The  model  has  a  solid  foundation  under  it.  Like  the  mean  and 
median  filters  are  a  subset  of  the  L  filter,  so  are  the  L  filter  and 
the  a-trimmed  linear  filter  special  cases  of  the  generalized  order 
statistic  model  outlined  in  Section  III.  As  with  any  new  design, 
however,  there  are  problems  to  be  overcome.  Indeterminate  gain 
for  some  of  the  filters  and  a  serious  degradation  in  frequency 
response  characteristics  for  nearly  all  the  examples  considered  here 
are  two  of  the  major  things  that  need  to  be  investigated  further  if 
the  a-TL  is  to  be  a  viable  filter  design.  As  this  simulation  showed, 
too  much  performance  is  given  up  in  the  frequency  response  of  the 
linear  filter  for  an  amount  of  outlier  resistance  already  achieved  or 
exceeded  in  other  filter  designs.  There  may  be  some  specific 
applications  where  an  a-TL  filter  would  be  the  best  choice,  but  it 
has  not  proven  worthy  as  a  general  filter  model,  vet. 


B  FUTURE  RESEARCH 


OPTIMIZING  THE  COEFFICIENTS 

Perhaps  the  most  remarkable  result  of  the  entire  simulation 
was  the  appearance  that  some  filters  were  affected  less  by  the 
a- trimming  process  due  to  their  having  a  steep  rolloff  design  in  the 
transition  band.  This  leads  to  the  conclusion  that  some 
characteristic  of  the  linear  filter  coefficients  can  be  exploited  in 
order  to  minimize  the  effects  that  a- trimming  has  upon  the  filter. 
No  conclusive  evidence  was  found  in  this  simulation  that  actually 
pinned  down  the  characteristic  or  characteristics  of  the  coefficients 
that  seemed  to  be  resistant  to  a- trimming.  An  area  for  further 
study  could  include  a  way  to  optimize  the  linear  filter  coefficients 
so  that  this  effect  may  be  achieved.  In  addition,  it  may  be 
possible  to  adaptively  modify  the  coefficients  to  optimize  their 
performance. 

HR  FILTERING 

Just  as  an  HR  filter  of  a  certain  length  can  easily  outperform 
an  FIR  filter  of  the  same  length,  perhaps  using  a  feedback  loop  in 
the  generalized  order  statistic  filter  model  could  improve  the 
frequency  performance  of  the  trimmed  linear  filter.  Using 
feedback,  the  I1R  linear  filter  can  be  designed  to  have  a  very  steep 
rolloff.  It  has  already  been  shown  that  some  relationship  exists 
between  the  steepness  of  the  designed  rolloff  and  the  immunity  of 
a  filter  to  a-trimming.  It  is  possible  that  an  HR  filter  may  prove 
to  be  more  resistant  to  the  damaging  effects  of  a-trimming.  This 


may  be  the  only  fix  necessary  to  adequately  preserve  enough  stop 
band  attenuation  during  the  trimming. 

The  GOS  filter  model  would  have  to  be  modified  slightly  in  order 
to  accommodate  an  UR  filter.  A  feedback  loop  on  the  time 
weighting  coefficients  is  obviously  necessary  in  order  to  have  an  HR 
filtering  operation. 

MODIFIED  TRIMMED  LINEAR  FILTER 

Looking  at  the  suggested  modification  outlined  above  makes  one 
wonder  what  would  happen  if  there  were  a  feedback  loop  on  the 
rank  weighting  coefficients  (bn's) .  This  modification  was  made  to 
the  o-trimmed  mean  filter  as  outlined  in  Section  II.  The 
modification  was  called  the  modified  trimmed  mean  (MTM)  filter 
and  performed  at  least  as  well  as  the  a-TM  filter. 

Recall  that  the  MTM  filter  provided  a  feedback  loop  that 
trimmed  a  number  of  elements  from  the  active  window  based  on  a 
range  parameter  q.  The  median,  mk,  was  selected  from  the 
window  and  then  all  elements  outside  of  the  range  [  mk  ±  q]  were 
trimmed  off.  Having  the  number  of  elements  trimmed  become 
dependent  on  the  data  allowed  the  filter  to  act  as  a  median  filter 
when  much  of  the  data  fell  outside  the  range  and  like  a  mean 
filter  when  there  was  very  little  deviation  from  the  median.  The 
range  parameter  was  chosen  based  on  a  priori  knowledge  of  the 
noise. 

This  scheme  may  prove  particularly  useful  in  the  case  of  the 
trimmed  linear  filter.  A  modified  trimmed  linear  (MTL)  filter  could 
outperform  an  a-TL  filter  for  two  reasons.  As  the  simulation 


82 

showed,  the  major  degradation  in  performance  occurred  when  just 
two  elements  were  trimmed.  If  the  signal  were  particularly  noisy 
with  many  impulses,  the  degradation  would  not  be  significantly 
greater  than  for  a  lower  fixed  number  of  samples  trimmed. 
However,  the  real  payoff  would  occur  during  the  times  when  no 
elements  were  trimmed.  The  MTL  filter  would  allow  for  the 
possibility  of  not  trimming  any  elements  during  portions  of  the 
signal  that  did  not  require  trimming.  This  would  allow  the  full 
effect  of  the  linear  filter  to  be  felt. 

By  choosing  q  such  that  an  element  was  trimmed  for  only 
extreme  cases,  a  filter  may  be  developed  that  has  some  resistance 
to  outliers  while  maintaining  a  good  frequency  response.  This  case 
would  be  extremely  interesting  to  pursue.  It  may  give  an  insight 
to  just  how  sensitive  a  linear  filter  is  to  trimming  of  any  kind. 

C.  CONCLUSIONS 

As  can  be  seen  from  the  development  of  this  model  and  hinted 
at  by  some  of  the  results  of  the  simulation,  the  trimmed  linear 
filter  holds  some  promise  in  living  up  to  its  expectations  as  a 
frequency  selective  filter  with  outlier  resistance  properties 
However,  as  the  simulation  also  showed,  the  a- trimming  process 
can  extract  a  fairly  high  price  from  the  frequency  performance  of 
a  linear  FIR  filter  for  a  little  resistance  to  outliers 

This  is  the  major  drawback  of  the  a  TL  filter  design  The 
section  on  outlier  resistance  showed  that  this  filter  teugr  ouki 
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have  a  comparable  degree  of  impulse  rejection  to  the  a- trimmed 
mean  filter  in  the  low  pass  case.  It  at  least  did  a  fair  job  at 
rejecting  outliers  in  all  of  the  cases  considered.  In  addition  the 
results  were  surprisingly  good  for  the  BPF  case.  Trimming  two  or 
four  elements  from  the  BPF  designs  in  the  simulation  produced  a 
marked  improvement  in  outlier  resistance  in  the  NMSE 
calculations.  Using  a  more  sophisticated  trimming  scheme  may  give 
even  better  rejection  of  impulses 

As  the  possibilities  discussed  in  section  B  above  show,  there  is 
still  hope  for  the  a-TL  filter  The  simulation  showed  several 
interesting  relationships  between  some  of  the  filter  desgpi 
parameters  and  how  the  a- trimming  process  affected  the 
performance  of  the  linear  filters  Some  of  these  mav  prr»v*  jseta 
in  helping  to  improve  the  performance  of  the  »  T*.  Mite-  * 
point  that  it  is  a  useful  des*n  for  more  thar  a  re*  we  * 
applications 
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APPENDIX  A 


This  appendix  contains  some  of  the  specific  details  used  in 
calculating  the  transfer  functions  for  the  simulation.  Several 
aspects  of  the  model  and  computer  system  require  some  detailed 
explanation  for  those  readers  interested  in  the  actual  development 
and  implementation  of  the  model.  The  details  described  in  this 
appendix  are  1)  Gaussian  and  Laplacian  random  number 
generation;  2)  the  histogram  of  calculated  data  and  its 
implications;  3)  window  selection  for  the  input  data;  and  4)  error 
calculation  to  give  an  idea  of  confidence  in  the  simulation. 

RANDOM  NUMBER  GENERATION 

The  Pascal  compiler  used  on  the  IBM  4341  system  provided  a 
random  number  generation  function  that  was  uniformly 
distributed.  In  the  GOSF  simulation  it  was  necessary  to  generate 
both  Gaussian  and  Laplacian  random  variables.  This  was  done  using 
a  transformation  of  variables. 

In  order  to  transform  variables  from  one  distribution  to 
another  it  is  necessary  to  set  the  probabilities  of  the  desired 
variable  equal  to  those  of  the  known  variable  as  shown  below  for 
the  case  of  a  known  uniform  rv  and  the  desired  Gaussian  rv: 

u  n 

J  dx  *  J  exp{-y2/2o2}dy/\/27to  (Al) 

0  -n 

where  u  is  a  uniform  rv  and  n  is  a  Guassian  rv  with  standard 
deviation  o.  However,  we  know  from  basic  probability  theory  class 
that  the  normal  distribution  cannot  be  solved  in  this  form. 


Therefore,  instead  of  solving  the  above  equation  we  solve 


u  r 

J  dx  =  J  (y/o2)exp{-y2/2o2}dy  (A2) 

o  o 

where  u  is  the  same  uniform  rv,  but  r  is  now  a  Rayleigh  rv.  This 
equation  can  be  manipulated  and  solved  for  r  as  a  function  of  u. 

The  Rayleigh  distribution  is  obtained  by  looking  at  a  two 
dimensional  space  that  is  Gaussian  along  x  and  y.  The  radial 
variable  of  this  space  has  a  Rayleigh  distribution  and  the  angle 
variable  is  uniformly  distributed.  Transforming  the  known 
uniformly  distributed  variable  into  a  Rayleigh  and  generating  an 
additional  uniform  rv  will  allow  an  indirect  transformation  to  two 
Gaussian  random  variables.  The  above  equation  is  solved  for  r  in 

t 

terms  of  u  by  substituting  variables,  integrating  and  solving.  The 
final  equation  becomes 

t 

i 

r  *  \/-2o2ln(l-u) .  (A3)  1 

{ 

\ 

1 

i 

Then  combining  this  Rayleigh  rv  with  a  uniformly  distributed  rv  0 

vy re  make  the  two  independent  Guassian  random  variables  ! 


The  same  approach  is  taken  for  the  Laplacian  (or  double¬ 
exponential)  case,  only  equating  the  probabilities  is  directly 
solvable.  The  equation  is 

u  d 

J  dx  *  J  (\Z2/2o)exp{-(>/2/o)fyl}<iy  (A6) 

o  -d 

where  u  is  the  uniform  rv  and  d  is  a  Laplacian  rv  with  standard 
deviation  o.  This  equation  is  solved  by  dividing  the  right  integral 
into  two  parts  to  get  rid  of  the  absolute  value  of  y,  integrating 
and  solving.  The  result  gives 

d  *  -(o/^2)ln(l-u).  (A7) 

The  results  for  either  case  can  be  quickly  checked  by  noting 
that  as  u— *1  then  r~*CO  and  d-»00  Likewise,  as  u  approaches  0 
then  r-*0  and  d-*0. 

A  test  case  was  run  for  the  Gaussian  distribution.  Twenty 
thousand  uniform  random  variables  were  generated  and  were 
transformed  to  two  sets  of  ten  thousand  Gaussian  random 
variables.  The  results  including  the  expected  Gaussian  percentages 
are  shown  in  Table  Al. 


Tabic  Al:  Random  Number  Transformation 


Cumulative  Cumulative  Cumulative 
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•uniform 

•Gaussian  1 

♦Gaussian  2 

Normal  % 

0.5 

2077 

3875 

3882 

0.3830 

1.0 

2000 

6834 

6881 

0.6826 

1.5 

1948 

8646 

8655 

0.8664 

2.0 

1905 

9539 

9526 

0.9544 

2.5 

1930 

9866 

9860 

0.9876 

3.0 

2080 

9978 

9968 

0.9974 

3.5 

2038 

9992 

9994 

0.9996 

4.0 

2044 

9999 

10000 

0.9999 

4.5 

1979 

10000 

10000 

1.0000 

5.0 

Total 

1999 

20000 

10000 

10000 

1.0000 

WINDOW  SELECTION 

One  of  the  most  difficult  problems  to  solve  in  setting  up  the 
simulation  was  the  method  for  calculating  the  transfer  functions 
for  the  a- trimmed  linear  filter.  As  was  explained  in  the  Chapter 
IV  section  on  transfer  functions,  it  is  not  a  simple  matter  to  find 
the  frequency  response  characteristics  of  a  nonlinear  filter  that  is 
designed  to  reject  impulses  on  the  input  signal.  An  attempt  to  find 
a  simple  deterministic  method  eluded  all  efforts.  For  this  reason  a 
probabilistic  measure  was  attempted. 

The  first  problem  encountered  in  using  Gaussian  noise  as  the 
filter  input  was  one  of  leakage.  The  transfer  functions  were 
adversely  affected  by  noise  near  the  ends  of  the  data  stream.  This 
problem  was  easily  alleviated  by  windowing  the  input  data 
Several  different  windows  were  tried.  The  Blackman  window  was 
chosen  based  on  its  performance  in  calculating  the  transfer 


function  of  a  linear  filter  with  none  of  the  elements  trimmed.  By 
visual  inspection  it  was  determined  that  by  using  the  Blackman 
window,  given  by 

w(n)  *  0.42  -  0.5cos(2im/(N-i))  +  0.08cos(4nn/(N-l))  (A8) 

for  0£n£N-l,  the  linear  filter  frequency  response  was  most  closely 
determined  with  the  random  input  transfer  function  calculation.  In 
fact,  as  was  discussed  in  Chapter  IV,  the  NMSE  would  go  to  zero 
when  using  just  two  trials  of  random  data  input.  This  is  due  to 
the  fact  that,  of  the  windows  chosen,  the  Blackman  window  had 
the  best  leakage  characteristics.  One  needs  to  take  into 
consideration,  however,  that  in  the  determination  of  the  transfer 
functions  for  the  a-TL  filters  the  relatively  poor  resolution 
characteristics  of  the  Blackman  window  function  may  adversely 
affect  the  results,  but  cannot  be  avoided. 

HISTOGRAM  DISTRIBUTION 

A  somewhat  more  interesting  problem  came  about  in  the  actual 
calculation  of  the  transfer  function.  All  of  the  data  about  transfer 
functions  presented  in  this  thesis  is  based  on  500  inputs  of  256 
samples  of  white  Gaussian  noise.  Table  A1  in  the  previous  section 
showed  how  well  the  noise  fit  the  Gaussian  distribution  and  Figure 
A1  shows  the  autocorrelation  of  2048  samples  of  this  noise.  The 
correlation  length  indeed  appears  to  be  very  near  zero. 

However,  the  first  attempt  to  calculate  the  transfer  function  of 
a  particular  filter  from  the  500  trials  was  done  by  averaging  each 


2048 


of  the  500  transfer  functions  on  a  point  by  point  basis.  This 
procedure  produced  transfer  functions  which  were  unexpectedly 
"noisy"  in  nature.  That  is  to  say,  even  though  the  basic  shape  of 
a  given  transfer  function  was  believable,  there  were  unexplained 
bumps  and  spikes  on  top  of  the  shape,  all  of  which  were  positive. 
Correcting  a  small  error  in  the  FFT  algorithm  and  fixing  a  problem 
with  the  number  of  significant  digits  in  the  data  generation 
program  only  alleviated  a  small  amount  of  the  apparent  noise  on 
the  basic  transfer  function  shape. 

An  amount  of  uncertainty  in  the  calculation  was  to  be  expected 
since  the  transfer  functions  were  generated  using  a  random  input. 
Indeed,  a  small  random  deviation  from  a  smooth  shape  was 
observed.  However,  on  top  of  this  fluctuation  were  unexplained 
spikes  at  various  frequencies  of  the  transfer  function.  Initially  this 
was  explained  as  a  frequency  selectivity  of  the  a-trimming  process 
or  perhaps  some  sort  of  resonance  was  introduced  by  the 
trimming.  But,  as  was  explained  in  Chapter  V,  the  a-TL  filter  was 
later  found  to  not  have  any  irregular  frequency  responses. 

The  answer  was  found  when  histograms  of  transfer  function 
values  at  given  points  of  the  function  were  generated.  This  showed 
that  the  distribution  of  values  at  each  point  was  not  very 
symmetric  over  much  of  the  transfer  function.  Only  in  the 
frequencies  of  the  pass  band  of  the  filter  was  the  histogram 
anywhere  near  symmetric.  For  most  of  the  rest  of  the  points  the 
histogram  had  the  majority  of  the  points  grouped  around  some 
relatively  small  value  with  an  exponentially  shaped  fall  off  of 


points  away  from  this  main  hump.  This  led  to  the  conclusion  that 
averaging  the  500  values  may  not  be  the  best  way  to  show  the 
central  tendency  of  all  of  the  values  at  each  point  of  the  transfer 
function.  This  conclusion  was  further  supported  by  more 
investigation  into  the  histograms  of  the  data  which  found  some 
extremely  large,  unexpected  values. 

Figure  A2  shows  a  transfer  function  for  a  band  pass  filter 
calculated  early  in  the  simulation  (the  function  is  symmetric  about 
0  Hz  at  sample  128)  using  point  by  point  averaging.  The  arrows 
point  to  individual  values  of  the  averaged  transfer  function  and 
identify  which  sample  they  are.  One  can  see  that  the  general 
shape  of  this  transfer  function  is  the  same  as  those  for  some  of  the 
band  pass  filters  in  Chapter  V  with  some  impulsive-type  noise 
added  to  it.  Figures  A3  through  A6  show  some  representative 
histograms  of  the  data.  Specifically  the  histograms  are  for  points 
91  though  94,  respectively.  The  vertical  scale  shows  number  of 
occurrences  in  a  particular  bin  with  the  horizontal  scale  being  ten 
times  the  actual  value  of  the  transfer  function  at  that  point  for 
that  trial,  (note:  the  transfer  function  in  Figure  A2  was  scaled  by 
a  constant  after  the  histograms  were  generated;  therefore,  the 
average  of  all  of  the  points  of  the  histogram  divided  by  ten  will 
not  necessarily  give  the  value  of  magnitude  shown  in  the  graph  of 
the  transfer  function.) 

All  of  the  histograms  of  the  data  in  the  stop  bands  are  of  this 
same  basic  shape.  They  show  that  the  data  is  skewed  to  the  right. 
Thus  the  central  tendency  of  these  distributions  may  be  better 


transfer  function 


Figure  A2:  Transfer  Function  of  BPF2500s  L31  (linear  scale) 


described  by  the  median  rather  than  the  mean.  However,  the  real 
evidence  for  this  conclusion  comes  from  noticing  the  trials  that  fall 
into  bin  100  in  the  data  for  samples  91,  92  and  93.  The  histogram 
generation  program  limited  the  value  of  any  individual  trial  to  10. 
Any  occurrence  greater  than  10  was  put  into  bin  100.  So  the 
actual  values  for  these  large  occurrences  were  printed  out.  For  the 
single  large  occurrence  in  sample  91  the  value  was  just  11.5. 
Likewise,  for  the  single  occurrence  in  sample  93  the  value  was 
13.1.  However,  for  sample  92  there  were  two  large  occurrences. 
One  was  only  11.5,  but  the  second  was  169.8,  nearly  15  times 
greater  than  the  next  largest  value!  Since  there  were  only  500 
trials  to  average  over  this  one  extraordinary  trial  was  able  to  add 
more  than  ten  percent  to  the  final  averaged  value  at  that 
particular  sample. 

This  was  the  main  reason  for  using  the  median  of  the  data  at 
each  point  to  determine  the  "most  likely"  value  for  the  transfer 
function.  Yet  this  brings  up  another  interesting  point;  what  kind 
of  certainty  does  this  calculation  provide?  Is  this  simple  approach 
to  finding  the  transfer  function  actually  valid  in  this  case?  These 
questions  are  discussed  in  the  next  section. 

ERROR  CALCULATION 

The  last  two  sections  bring  up  some  questions  about  the  validity 
of  the  transfer  function  calculations  used  in  this  thesis.  The 
windowing  of  the  input  data  is  necessary  but  forces  a  choice 
between  some  leakage/resolution  trade-offs  of  various  windows.  The 
somewhat  strange  histograms  of  the  calculated  data  along  with  the 


peppering  of  large  outliers  in  this  data  introduces  some  additional 
uncertainty  into  the  measurement  of  the  transfer  functions  of  the 
a-TL  filters.  The  question  is  how  far  off  are  the  measured  transfer 
functions? 

The  choice  of  the  window  to  use  on  the  input  data  was 
straightforward.  Of  all  the  windows  tried,  the  Blackman  window 
did  the  best  job  at  reducing  the  leakage  problems  encountered  in 
the  transfer  function  calculations.  The  different  degrees  of 
resolution  between  the  various  windows  made  no  real  visible 
difference  under  the  criteria  selected  for  choosing  the  best  window. 
As  stated  before,  the  criteria  was  to  match  the  impulse  response 
frequency  characteristics  of  a  linear  filter  with  the  random  input 
transfer  function  calculation  of  the  same  untrimmed  linear  filter. 
The  Blackman  window  performed  the  best  since  it  had  the  best 
leakage  characteristics  of  all  the  windows  tested.  Once  the  leakage 
was  reduced  enough  the  transfer  function  matched  almost 
identically  with  only  one  trial  of  random  input  to  the  expected 
frequency  response  of  the  linear  filters.  Therefore,  it  is  assumed 
that  it  would  be  difficult  to  obtain  better  results  for  the  transfer 
functions  with  respect  to  the  window  used. 

There  still  may  be  some  question  in  regard  to  the  method  of 
"averaging”  the  transfer  function  trials,  however.  In  order  to 
address  this  question  some  confidence  intervals  were  found.  During 
the  median  selection  it  was  a  relatively  simple  manner  to  choose 
other  values  at  each  point.  Various  “confidence  intervals”  were 
generated  in  this  manner  by  selecting  the  10  and  90  percent 


values  to  give  an  80  percent  confidence  interval,  for  example.  This 
process  gave  a  picture  like  that  shown  in  Figure  A7  for  the  skinny 
band  pass  filter  of  15  elements  with  2  trimmed,  BPF2500s  15T2. 
The  figure  shows  the  ideal  transfer  function  of  the  linear  filter  and 
the  transfer  function  of  the  filter  with  2  elements  trimmed 
selected  by  the  median  of  500  trials  of  Guassian  noise  as  the  input. 
The  figure  also  shows  10  and  90  percent  lines  for  the  transfer 
function  with  2  elements  trimmed.  These  are  the  50th  and  450th 
largest  values  on  a  point  by  point  basis  of  the  ordered  data  from 
the  500  individual  trials.  They  are  normalized  by  the  same 
constant  as  was  calculated  to  normalize  the  median  response  to 
OdB  in  the  pass  band.  This  gives  a  good  idea  of  the  range  of  values 
this  method  of  calculating  the  transfer  functions  covers.  One  can 
see  that  if  the  10  and  90  percent  lines  were  properly  normalized 
they  would  come  fairly  close  to  being  identical  to  the  median  line. 
The  only  major  difference  is  the  wide  variance  at  the  dc  point  as 
was  discussed  in  Chapter  V. 

A  similar  picture  is  shown  in  Figure  A8  for  the  same  filter  with 
6  elements  trimmed.  The  major  difference  between  the  two  figures 
is  the  expected  increase  in  variance  between  the  10  and  90  percent 
lines  for  the  filter  with  6  elements  trimmed.  One  may  also  note 
that  it  appears,  by  taking  into  account  the  possibility  of  getting 
the  response  of  the  10  percent  line  in  the  pass  band  and  that  of 
the  90  percent  line  in  the  stop  band,  the  frequency  response  for 
the  filter  may  be  nearly  flat!  In  practice,  however,  one  finds  that 
this  is  not  true.  For  an  individual  trial  of  Gaussian  noise  input  the 
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Figure  A7:  80%  Confidence  Interval  for  BPF2500s  15T2 
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transfer  function  always  maintained  the  same  basic  transfer 
function  shape.  The  shape  was  simply  much  more  "noisy”  than 
the  relatively  smooth  response  given  by  the  median  of  500  trials. 
These  confidence  intervals  merely  show  the  magnitude  of  possible 
variations  within  a  transfer  function  for  a  single  input. 

This  explanation  makes  sense  when  one  takes  into  account  the 
nonlinearity  of  the  a-TL  filter.  Just  as  for  the  simple  median 
filter,  it  is  not  possible  to  predict  what  the  exact  frequency 
response  of  the  filter  will  be  for  any  arbitrary  input.  It  can  only 
be  said  that  the  filter  will  have  a  general  frequency  response  most 
likely  represented  by  the  shape  given  by  the  median  of  500  trials 
of  random  input  as  shown  in  this  thesis.  This  is  just  the  same  as 
saying  that  the  median  filter  is  a  low  pass  operation.  The  median’s 
exact  frequency  response  cannot  be  determined  without  knowing 
the  input.  It  is  only  known  in  general  that  it  will  act  in  a  low 
pass  fashion  with  an  approximate  cutoff  and  somewhat  variable 
stop  band  performance. 


APPENDIX  B 


In  order  to  investigate  the  characteristics  of  the  generalized 
trimmed  linear  filter  I  developed  a  program  that  would  implement 
the  various  types  of  filters,  a  program  to  generate  input  signals  to 
filter,  and  some  programs  to  compare  performances  of  the  various 
filters  on  a  set  of  input  signals.  The  filter  model  program  is  simply 
an  implementation  of  the  model  described  in  the  previous  section. 
All  of  the  filters  we  are  interested  in  using  as  comparisons  are 
subsets  of  this  generalized  model. 

LANGUAGE  SELECTION  AND  MODEL  STRUCTURE 

I  chose  to  write  the  program  in  Pascal  instead  of  the  perhaps 
more  obvious  choice  of  the  more  computationally  efficient  FORTRAN 
for  two  reasons.  First  was  that  the  model  structure  seemed  most 
logically  implemented  using  pointers  with  a  linked  list  architecture. 
FORTRAN  does  not  have  a  pointer  data  type  and,  therefore,  a 
linked  list  structure  would  be  difficult  to  implement.  Second,  1  like 
the  stronger  data  type  checking  in  Pascal  or,  rather,  1  dislike  the 
"weak  typing*  of  FORTRAN.  This  was  an  important  consideration 
due  to  the  moderate  complexity  of  the  model  combined  with  the 
slight  degree  of  complexity  of  the  pointer/linked  list  structure  to  be 
implemented. 

The  reason  I  chose  a  linked  list  structure  was  that  it  most 
efficiently  allowed  me  to  add  and  delete  elements  from  the  active 
filter  window  and  provided  a  relatively  efficient  sorting  routine. 
Another  good  choice  might  have  been  to  use  an  array  structure  for 
the  windows,  but  1  discarded  this  as  less  efficient  in  the  two  main 
manipulations  of  the  windows:  adding  and  deleting  elements  and 


107 


sorting.  The  efficiency  of  the  sorting  routine  becomes  important  in 
this  model  since  the  window  must  be  re-sorted  at  each  step 
through  the  input  signal.  The  insertion  sort  routine  for  a  linked  list 
structure  takes  the  first  value  and  calls  it  the  smallest  value.  Each 
successive  value  to  be  sorted  is  inserted  into  the  proper  place  in 
the  sequence  taking  care  to  keep  track  of  the  smallest  value. 
PROGRAM  SET  DESCRIPTION 

The  following  is  a  brief  description  of  each  of  the  main 
component  programs  used  in  the  simulation  as  well  as  a  discussion 
as  to  how  they  were  combined  to  form  the  large,  repetitive 
program  MEGAMEDIAN.  The  descriptions  here  are  limited  to  the 
major  Pascal  programs.  Many  smaller  usually  machine  dependent 
programs  were  written  to  produce  graphical  outputs  and  to  better 
control  the  data  file  structure  on  the  IBM  4341.  Since  these 
programs  would  be  of  little  use  to  the  general  programmer  they 
have  been  omitted  here. 

MAKEDATA:  This  program  allows  the  user  to  generate  an  input 
sequence  that  is  to  be  filtered.  The  input  may  have  any 
combination  of  the  following  possible  components:  constant  levels, 
ramps,  steps  or  impulses,  sinusoids  or  blocks  of  sinusoids,  Gaussian 
noise,  or  Laplacian  noise.  In  addition,  the  user  controls  many  of 
the  variables  associated  with  each  of  these  signal  components 
including  location  in  the  signal,  magnitude,  duration,  and 
variance,  to  name  a  few.  The  user  also  picks  the  number  of 
samples  in  the  signal.  This  value  is  stored  as  the  first  data  point  in 
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the  data  file  and  is  used  by  the  filtering  and  Fourier  transform 
programs  to  properly  handle  the  signal.  In  addition,  once  the 
desired  signal  has  been  generated  the  user  may  store  it  in  one  of 
two  different  files  so  that  one  may  be  used  as  a  reference  while 
the  other  is  varied.  The  user  also  has  the  option  to  go  back  to 
either  of  these  two  files  in  the  future  and  add  onto  the  existing 
signal. 

FREQ:  This  program  takes  the  discrete  Fourier  Transform  (DFT) 
of  the  input.  The  DFT  is  implemented  with  the  decimation-in- 
frequency  fast  Fourier  transform  (FFT)  algorithm.  The  user  is 
prompted  for  which  data  file  to  use  as  the  source  for  the  Fourier 
transform.  The  options  are  either  of  the  two  signal  files  created  by 
MAKEDATA,  the  output  of  the  autocorrelation  program,  Rxx,  or 
the  output  of  the  filter  program.  If  the  input  signal  is  not 
factorable  by  2  then  the  remaining  elements  are  padded  with 
zeros.  After  the  FFT  is  taken,  the  real  and  imaginary  parts  are 
converted  to  magnitude  and  phase  since  this  simulation  is 
primarily  interested  in  the  magnitude  of  the  transform.  The  user 
is  then  given  the  option  to  have  the  output  converted  to  a  log 
scale.  Finally,  the  user  has  several  options  to  store  the  output  in 
different  files  for  later  manipulations  and  calculations. 

FILTER:  This  program  obviously  implements  the  filter  model 
discussed  in  Chapter  III.  It  is  only  important  here  to  note  the 
inputs  the  user  may  make.  The  user  controls  the  window  size  to 
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be  used  and  then  selects  the  type  of  filtering  operation  to  perform. 
Once  this  choice  has  been  made  the  program  prompts  the  user  for 
the  necessary  inputs  to  complete  the  filtering  operation. 

MSE:  This  program  calculates  not  only  the  NMSE  between  two 
inputs,  but  also  the  NAE.  In  addition  the  program  generates  two 
output  files  to  be  graphed  The  first  merely  combines  the  two 
inputs  into  a  single  file  so  that  both  graphs  may  be  plotted  at  the 
same  time.  The  second  file  is  a  graph  of  the  difference  signal 
between  the  two  inputs.  These  two  graphs  are  useful  not  only  for 
simulation  insights,  but  also  for  troubleshooting.  The  user  has  a 
wide  range  of  choices  to  select  which  data  files  to  calculate  the 
errors  between.  The  program  can  handle  files  stored  in  log 
magnitude  as  well. 

MEGAMEDIAN:  This  file  creates  transfer  functions.  MAKEDATA, 
FREQ,  and  FILTER  are  all  procedures  in  this  very  large  program. 
The  main  program  loops  through  these  three  programs  as  well  as  a 
procedure  to  calculate  the  transfer  function  for  each  of  the  5 00 
trials  in  order  to  generate  a  transfer  function  for  a  particular 
filter. 

The  listings  for  these  five  programs  are  Appendix  D. 


Table  Cl  :  Filter  Coefficients  for  LPF300  L9 


************************************  *  *  *  t  %  *  $  ****  *  ******  r%  t  ; 

FINITE  IMPULSE  RESPONSE  (FIR  > 
LINEAR  PHASE  DIGITAL  FILTER  DESIGI 
REMEZ  EXCHANGE  ALGORITHM 


BANDPASS  FILTER 


FILTER  LENGTH  = 


*****  IMPULSE  RESPONSE  ***** 


H  v 

1  > 

= 

.  2 3 3 2 8 2 9 0 E  •*-  0  0 

- 

H 

H  •; 

. '  '» 

= 

. 45690290E-0 1 

= 

H  *. 

C 

H  < 

3  • 

= 

.  484  55°40E— ■ :>1 

ar 

H  • 

H  < 

A  . 

- 

.  50 17917 OE— “ 

= 

H  * 

L. 

H  ( 

5  > 

= 

.  So^LBSOOE— 0 1 

— 

H  ( 

er 
„  * 

LOWER  BAND  EDGE 
UPPER  BAND  EDGE 
DESIRED  VALUE 
WEIGHTING 
DEVIATION 
DEVIATION  IN  DB 


BAND  1 

. 000000000 
.  0 1  ooooooo 
1 . 000000000 
2. 000000000 

.211 253*00 


-13. 503900000 


BAND.  2 

. 050000000 
50C^000'  )/*•*’■ 
. 000000000 
1 • 000000000 

.  42250*^300 
-7.48730700^ 


EXTREMA  FREQUENCIES 

.  0100000  .  050000-..' 

. 5 00000 C 


*.  ********  *  *  *  *  *  *  *  *  *  %  *  *  *  *  *  *.  *  *  *  *  *r*  *  *  *»**  *  *  * :+  *  *  *  1  *  ' 


Table  C2  :  Filter  Coefficients  for  LPF300  L15 


***********************  ********************  ************** 

FINITE  IMPULSE  RESPONSE  (FIR) 
LINEAR  PHASE  DIGITAL  FILTER  DESIGN 
REME2  EXCHANGE  ALGORITHM 


BANDPASS  FILTER- 
FILTER  LENGTH  =  15 


*****  IMPULSE  RESPONSE  ***** 


H  ( 

1 ) 

= 

.  1 1714120E+00 

=  H  ( 

15) 

H  ( 

2) 

. 40280550E— 0 1 

=  H  ( 

14) 

H  ( 

3 ) 

= 

. 45621000E-01 

=  H  ( 

13) 

H  ( 

4) 

. 50370080E-01 

=  H  ( 

12) 

H  ( 

5) 

= 

.  543 1 0600E— 0 1 

=  H  ( 

11) 

H  ( 

6) 

= 

. 57236030E—01 

=  H  ( 

10) 

H  ( 

7) 

= 

. 589853 10E-01 

*  H  ( 

9) 

H  ( 

8) 

= 

. 595809 10E-01 

=  H  ( 

8) 

BAND 

1  BAND 

r> 

LOWER  BAND  EDGE 
UPPER  BAND  EDGE 
DESIRED  VALUE 
WEIGHTING 
DEVIATION 
DEVIATION  IN  DB 


. 000000000 
. 010000000 
1 . 000000000 
1 . 500000000 
. 131174300 
“17. 6430 1 0000 


. 050000000 
. 500000000 
. 000000000 
1 . 000000000 
. 196761400 
-14. 121190000 


EXTREMA  FREQUENCIES 

. O 1 00000  . 0500000  . 085 1 562  .1515625  .21 796S 

.2882813  .3585938  .4289063  .5000000 


******************************************************* 
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Table  C3  :  Filter  Coefficients  for  LPF300  L31 


****************************** * * * t * * *  *  *  * *  * *  *  *  * %  * **** *  # #  * t 

FINITE  IMPULSE  RESPONSE  (FIR) 
LINEAR  PHASE  DIGITAL  FILTER  DESIGN 
RE ME 2  EXCHANGE  ALGORITHM 

BANDPASS  FILTER- 

FILTER  LENGTH  =  31 


IMPULSE 

h  (  i  :■  = 

RESPONSE  ***** 
-. 1 1461630E—01 

— ■ 

U  ( 

T  4  V 

HI  2) 

=5 

. 1 4994630E—0 1 

= 

H  ( 

30) 

H  (  3) 

as 

. 1 40 1 0050E-0 i 

as 

H  ( 

H  (  4  • 

= 

•  1  o f-  J  3  <£>  5  i  J  E — o  1 

= 

H  ( 

28) 

H  (  5) 

s= 

. 1 9599590E-0 1 

= 

H  ( 

27) 

H  (  6) 

= 

.23948930E-01 

as 

H  ( 

26) 

H  (  7) 

=r 

. 287341 10E -01 

sr 

H  ( 

25> 

H  (  8) 

= 

. 3368587 OE-0 1 

= 

H  ( 

24) 

H  (  9) 

= 

. 386 1 7050E-0 1 

as 

H  ( 

0*7;  \ 

H  ( 2. 0 ) 

= 

. 433 1 B480E— 0 1 

= 

H  ( 

OO  'j 

H  (11) 

. 47660 100E-01 

= 

H  ( 

1  s. 

H  ( 12) 

= 

.51 456090E-0 1 

= 

H  ( 

20) 

H<  13) 

= 

. 54554B90E-03 

s 

H  ( 

]  ?) 

h  <  14  :■ 

= 

. 56882300E— 01 

St 

H  ( 

18) 

H  <  15) 

= 

. 5S2B5B60E-01 

= 

H  ( 

17) 

H  ( 1 6 ) 

=s 

. 587401 80E-02 

Hi 

16) 

LONER  BANE!  EDGE 
UPPER  BAND  EDGE 
DESIRED  VALUE 
WEIGHTING 
DEV I AT I ON 
DEVIATION  IN  DB 


BAND  1 
. 0  jOOOO OOO 
. 0 1 0000000 
1 . 000000000 
1 . 000000000 
. 039385330 
■•g .  093170000 


BAND  2 
. 050000000 
. 500000000 
. 000000000 
1 . 000000000 
.  7JOJO -J'C- 

— 2  S . 0 9317 O 0 0 0 


EXTREMA  FREQUENCIES 

. 0000000  . 01 00000  . 0500000 
.126171?  .1 593750  . 1 9453 1 2 
.2960938  .3312500  .3644532 
. 4660 1 56  . 5000000 


06367 1 9 
2277344 
3976563 


A  —*■  O  -<  ■— i  P 

•  *+  w  -i.  CD  J  ■. 


********************************  *  *  *  *  %  %  #  %  %  %  t  *  *  * t  *  *  *  *  *  *  #  *  t  *  t  8 


j£ 


$ 

§s 

n*. 

1 

4 


NVWW.W\ 


.v.v.v.v.v.v.v 
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FINITE  IMPULSE  RESPONSE  (FIR) 
LINEAR  PHASE  DIGITAL  FILTER  DESIGN 
REMEZ  EXCHANGE  ALGORITHM 

BANDPASS  FILTER 

FILTER  LENGTH  »  63 

*****  IMPULSE  RESPONSE  ***** 


H  (  1 ) 

* 

-. 1 B686920E— 02 

x 

H  ( 

63) 

H  <  2) 

X 

1 5462900E— 02 

x 

H  ( 

62) 

H  (  3) 

x 

-.21 003040E-02 

m 

H  ( 

61) 

H  (  4) 

x 

-. 270696S0E— 02 

x 

H  ( 

60) 

H  <  5) 

* 

-. 33382 1B0E-02 

* 

H  ( 

59) 

H  (  6) 

* 

-. 395672 10E-02 

x 

H  ( 

58) 

H  (  7) 

s 

- . 45 1 67 690E— 02 

m 

H  < 

57) 

H  (  8) 

= 

-. 49655480E— 02 

x 

H  ( 

56) 

H  (  9) 

* 

— . 52443800E-02 

s 

H  ( 

55) 

H  ( 1 0 ) 

X 

— . 52912720E— 02 

x 

H  ( 

54) 

H  ( 1 1 ) 

X 

-.  504377 10E-02 

x 

H  ( 

53) 

H  ( 12) 

= 

- . 444 1 6870E— 02 

= 

H  ( 

52) 

H  ( 1 3) 

SB 

-. 34307 100E-02 

x 

H  ( 

51) 

H  <  14) 

as 

-. 1 9657400E— 02 

x 

H  ( 

50) 

H  ( 15) 

X 

-. 13794770E— 04 

m 

H  ( 

49) 

H  ( 16) 

x 

. 24416860E— 02 

x 

H  ( 

48) 

H  ( 17 ) 

■ 

. 540 1 6280E— 02 

x 

H  ( 

47) 

H  ( 18) 

* 

.  88445020E-02 

s 

H  ( 

46) 

H  ( 19) 

X 

.  1 273 1 530E— 0 1 

x 

H  ( 

45) 

H  (20) 

e 

. 1 7 0006B0E— 0 1 

X 

H  ( 

44) 

H  (21 ) 

as 

.21 57 1 820E— 0 1 

X 

H  ( 

43) 

H  (22) 

= 

. 263476 10E-01 

X 

H  ( 

42) 

H  (23) 

x 

. 31215140E— 01 

X 

H  ( 

41) 

H  (24) 

s 

. 3605 1 430E— 0 1 

X 

H  ( 

40) 

H  (25) 

= 

. 40727 130E-01 

X 

H  ( 

39) 

H  (26) 

X 

. 45 1 1 0720E— 0 1 

X 

H  ( 

38) 

H  (27) 

x 

.  49074520E— 01 

X 

H  ( 

37) 

H  (28) 

s 

. 52500080E— 0 1 

X 

H  ( 

36) 

H  (29) 

x 

. 55282880E— 01 

X 

H  ( 

35) 

H  (30) 

= 

.  57335970E— 01 

X 

H  ( 

34) 

H  (31 ) 

as 

.  58594040E— 01 

X 

H  ( 

33) 

H  (32) 

x 

.  59018030E-01 

X 

H  ( 

BAND  1  BAND  2 

. OOOOOOOOO  . 050000000 

. 0 1 0000000  . 500000000 

1 . OOOOOOOOO  . OOOOOOOOO 

1 . OOOOOOOOO  1 . OOOOOOOOO 

. 002450794  . 002450794 

52.213810000  -52.213810000 

EXTREMA  FREQUENCIES 

.0067708  .0100000  .0500000  .0536458  .0630208 

.0755208  .0890625  .1036458  .1187499  .1343751 

.1494794  .1651046  .1807298  .1968759  .2125011 

.2281263  .2442724  . 259B973  .2760428  .2921884 

.3078132  .3239588  .3401044  .3557291  .3718747 

.3880203  .4041659  .4197907  .4359362  .4520818 

. 4677066  . 4838522  . 5000000 


LOWER  BAND  EDGE 
UPPER  BAND  EDGE 
DESIRED  VALUE 
WEIGHTING 
DEVIATION 

DEVIATION  IN  DB  - 


Table  C5  :  Filter  Coefficients  for  LPF1500  L9 


**********************************  ************************ 

FINITE  IMPULSE  RESPONSE  (FIR'* 
LINEAR  PHASE  DIGITAL  FILTER  DESIGN 
REME2  EXCHANGE  ALGORITHM 

BANDPASS  FILTER 

FILTER  LENGTH  =  9 

*****  IMPULSE  RESPONSE  ***** 


H  ( 

1  '  - 

— . 47385260E— O 1 

= 

H  - 

Q 

H  ( 

2)  = 

.  51 096300E-01 

= 

H  - 

C; 

H  ( 

)  — 

.  1530725 0 E + 0 0 

= 

H  « 

“7 

H  ( 

4)  * 

. 25925620E+00 

— 

H  < 

£> 

H  ( 

5 )  = 

.  30436960E+00 

sr 

H  < 

cr 

LOWER  BAND  EDGE 
UPPER  BAND  EDGE 
DESIRED  VALUE 
WEIGHTING 
DEVIATION 
DEVIATION  IN  DB 


BAND  1  BAND  2 

. OOOOOOOOO  . 200000000 

. 1 00000000  . 500000000 

1 . OOOOOOOOO  . OOOOOOOOO 

1 . OOOOOOOOO  1 . 300000000 

. 136449300  . 104961000 

-17. 300560000  -19. 579420000 


EXTREMA  FREQUENCIES 

.0000000  .1000000  .2000000  .  2633737  .  ’'efo-j 

. 500000C 


*.*#*  *  *  ********  *  ****  *  *  ***  *  *  *  *  *  *  *  *  *  *  *  *  ****** *  *  *  *  a 
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Table  C6  :  Filter  Coefficients  for  LPF1500  L31 


*  *  *  *  *  *  *  *  *  *  *  *  *  *  *  ♦  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  * #  *  *:  *  *  *  *  *  *  *  *  *  *  *  *  *  *  * 

FINITE  IMPULSE  RESPONSE  (FIR) 

LINEAR  PHASE  DIGITAL  FILTER  DESIGN 
REME2  EXCHANGE  ALGORITHM 

BANDPASS  FILTER 

FILTER  LENGTH  =  31 


I MPUL 

SE 

RESPONSE  ***** 

H  <  1  / 

= 

. 1 895&240E— 02 

- 

H  C 

31  ‘ 

H  (  2) 

* 

. 1679o450E-02 

H  t. 

30  ? 

H  (  3) 

= 

-. 1 34&7520E— 02 

* 

H  < 

3*9 

H  (  4 > 

= 

o2655380E-02 

=r 

H  ' 

28  i 

H  (  5  > 

= 

-  .  7<s84040oE— 02 

= 

H  < 

■  ■ 

H  <  6 ) 

= 

3389007 OE— 03 

= 

H  ( 

26  ' 

H  (  7) 

= 

. 1351 1050E-01 

= 

H  ( 

25 ) 

H  (  8 ) 

= 

. 21470 150E-01 

s 

H  ( 

24 

H  '  9) 

= 

. 946698 10E-02 

at 

H  ( 

23) 

H  *  1  0  > 

= 

2250991 OE-01 

as 

H  • 

'  ‘  ■  *  » 

H  ( 1 1  ) 

= 

-.51 376640E-0 1 

as 

H  ( 

21  ) 

H  ( 12) 

* 

4 1 209470E-0 1 

as 

H< 

20 ) 

H  (1 3  > 

= 

.29889860E-01 

as 

H  ( 

19) 

H  ( 1 4  > 

= 

. 1 46 1 6060E+00 

s 

H  ( 

18) 

H  (  1 5 ) 

= 

. 25564390E+00 

s 

H  ( 

17) 

H  (  1  6  ) 

= 

. 30050370E+00 

s 

H  ( 

1  6 ) 

LOWER  BAND  EDGE 

BAND  1 

.  OOoOOoOO'- ' 

BAND  2 

•  200000000 

UPPER  BAND  EDGE 

•.  1  •!*  •*  *  0  (")  O  «« 0 

DESIRED  VALUE 

1 . OoOOOOOOO 

. OOOOOOOOO 

WEIGHTING 

1 0 . OOOOOOOOO 

1 0. OOOOOOOOO 

DEVIATION 

. 00 \ 523 1 95 

. 001523195 

DEV In T ION  IN  DB 

-5o.  344830! >00 

-56. 34483000'. 

E  X 1 REMa  FREQUENC I ES 

. 0000000  . 0437500  . 0 

718750  .0927033 

1  Mi  li  It  |i  if  • 

.  200000*. 

.207291“  .2 

260416  .2520833 

. 2302089 

. 31 041 “4 

. 34062&2  . 3 

■'29182  .4041687 

.  4 3 o 4  e»OA 

.  467"’  |  11 

. 5000000 

**********  *  *  *  t  ********************  *  *  *  *  *  $ :  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  * 


■  «  ».«  lu  «.»  <•■«  «-»  «.»  j»L«iLlitliil 


FINITE  IMPULSE  RESPONSE  (FIFO 
LINEAR  PHASE  DIGITAL  FILTEP  DESIGN 
REME2  EXCHANGE:  ALGORITHM 


BANDPASS  FILTER 
FILTER  LENGTH  =  63 


1 ITPUL 

,.3E 

RESPONSE  * *  *: # * 

H  (  1  ) 

= 

1 14972 10E-O4 

= 

H  (  63 .1 

H  (  2  > 

ss 

— .  40 1 77320E— < >5 

= 

H  (  62) 

H  (  3) 

r: 

. 40 1 04340E-04 

= 

H  (  6 1  > 

H  (  4 ) 

= 

.  9473 0 3 2 0 E — 0  4 

= 

FI  <  60  > 

Hi  5  ’■ 

=: 

. 6  J  52972 OE— 04 

ss 

H<  59) 

H  (  6  > 

s 

-- .  i  3574090E— 03 

= 

H  *  55 ) 

H  7  ■ 

= 

-• .  3846042DE  -  03 

= 

-!<  57  • 

H  •  8  ■ 

ss 

3538 1 9  30E -03 

= 

Hi  56  > 

11,  r~.  . 

n  •  1  ■’ 

= 

.  21 45867  0E--03 

= 

Pi 1  55,' 

H  ( 1 0  > 

r= 

. 1048450OE-O2 

= 

M  (  54  > 

H  <  .1  1  • 

S 

1 2  f? 2 0 2 6 0 E  --  0  2 

Fi  ■'  53  ■ 

H  <  1  2 ) 

as 

.  643651  1  OE--04 

-• 

Fh'  52) 

H  (  i  3  "> 

=r 

-  ..  2 1 220690E-  02 

i|  (  5 1  > 

H  i  1  4  ) 

= 

- . 3333S030F-02 

ss 

H  i  50  > 

H  :  1  5) 

-  .  :!  4 7594 2 OE -02 

~ 

FI  <  49) 

H  ( 1 6  > 

= 

.321 931 00E-02 

= 

H<  48  » 

H  ( 1 7  ;■ 

ss 

. 70782 /60E— 02 

= 

Hi  47) 

h  ( 1  a :» 

= 

-  52989580E--02 

s= 

H  •:  46 ) 

H  (  1  9  ■. 

s 

-.31 939230E-02 

= 

H  (  45) 

H  •;  20  > 

= 

12623200E-01 

ss 

H*  44) 

H  ( 2: 1  > 

= 

1  3 3 0 247 0 E - 0 1 

= 

H(  43) 

H  ( 22 ) 

= 

22721 440E-03 

= 

H  (  42) 

H  t  23  > 

s 

. 1 94077 1 OE-0 1 

= 

H  (  4  1 

H  < 24 

= 

. 28 1 24250E-O  3 

= 

H  i.  40  ) 

H  (25) 

zz 

. 1 J 444060E— 0 1 

ss 

H  (  39 ) 

H  ( 26  * 

== 

— . 26 1 39 S 1 OE ~0 1 

ss 

H '  78' 

H <  2“  > 

r:: 

- .56691 01 OE~oi 

- 

J„  !  f  "7  i 

-J  /  •—  —  •*.  .1 

— 

•  •  ““ 

1-1  46 

=- 

*"  «  T  -•  C?  i-"'  w  w  .1 

F  ( 29 ) 

=r 

.31 1 4443 OE —01 

H  f  35- 

H  • 

•  j  469  v4  *  oE  '  * 

H  (  3  4  •• 

H  ( 3 :!  > 

n 

. 2565237 OE  -Ou 

ss 

H  (  73  • 

H  *  72  • 

ss 

. 300334  50E+00 

= 

H  C  32  • 

BAND  1 

LOWER  BAND  EDGE  .  OOOOO  "OO0 
UPPER  BAND  EDGE  . 1 00000000 
DESIF'EI)  VALUE  J  .  000000000 
WE I GHT I MG  J o . OOOOOOOOO 
DE  I  AT  I  ON  .  0Gm0080~0 
DEV I AT I ON  I N  DB  -  1 0 1 . 962200 000 


BANE,1  A 

.  2  O 0 0 0 0 0 0 0 
, 5000000 00 
,  OOuOOOO'  n~i 
I  0. 000000'  ;0'0 
.  fi  ‘)0009<  90 
10  ;  .  S £>22000 00 


FREQUENCIES 


oOOO  >’  >0 
0708333 
200000*0 

,  241 1464 
31 30205 
,  3°0 1 C 27 
,  46974  ~3 


0151 "42 
08229  I  6 
2020834 
2546881 
3231245 
4057274 
4843*T2 1 


. 032 2033 
„  O 9 1 66 6 6 
. 2083334 
. 2687504 
. 343749: 
. 4213522 
. 5000000 


OF'—'-  i  hr 


4369' 


. j  r  ooooo 

,,  22  36462 
. 3979 T6t 
..  3-44  ~?p 


& 


& 


w 

jy 

I 


Ci.f 

M 


1 
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Tabic  C8  :  Filter  Coefficients  for  BPF2500s  L9 


t**Xtt***t*t*%%*tt*tt%%tttt*t**t%*tZt*%*****%t*t**%*tt**tttii 


FINITE  IMPULSE  RESPONSE  (FIR) 
LINEAR  PHASE  DIGITAL  FILTER  DESIGN 
REME2  EXCHANGE  ALGORITHM 


BANDPASS  FILTER 
FILTER  LENGTH  =  ? 

*****  IMPULSE  RESPONSE  ***** 


H  ( 

1 ) 

= 

.  2 0 3846 0 0  E + 0  0 

5= 

H  ( 

?  > 

H  ( 

2 ) 

= 

— • 3 066 9 860 E — 07 

= 

H  ( 

8) 

H  ( 

3 ) 

= 

1 4731 930E+00 

as 

H  ( 

7 ) 

H  ( 

4) 

= 

25204350E-09 

= 

H  ■: 

6*  / 

H  ( 

5) 

= 

. 161 203 1 OE+OO 

H  < 

5) 

LOWER  BAND  EDGE 
UPPER  BAND  EDGE 
DESIRED  VALUE 
WEIGHTING 
DEVIATION 
DEVIATION  IN  DB 


BAND  1 

ooooooooo 
1 80000000 
ooooooooo 
ooooooooo 

11. 268580000 


1 


BAND  2 
. 230000000 
. 270000000 
1 . OOOOOOOOO 
1 . 400000000 
. 195183200 
■14.1911 40000 


BAND  3 
. 320000000 
. 500000000 
. OOOOOOOOO 
1 . OOOOOOOOO 
. 273256500 
-11. 268580000 


EXTREMA  FREQUENCIES 

. 0000000  . 1 100000 
. 3899999 


. 1900000 


■  2  i  oo  o  o  o 


■  OOOOO 


l] 


s 


Table  C9  :  Filter  Coefficients  for  BPF2500s  L15 


****************** ************** *******$  1 1  *  ******  * : 


FINITE  IMPULSE  RESPONSE  (FIR) 
LINEAR  PHASE  DIGITAL  FILTER  DESIGN 
REMEZ  EXCHANGE  ALGORITHM 

BANDPASS  FILTER 


ER  LENG 

;th 

=  15 

*  I MPUL 

□  L 

REbPUNSE  *  *  *  X  it 

Hi  i  > 

= 

. 1 32674 10E-04 

= 

H  •' 

15) 

H  (  2) 

= 

60674 130E -01 

= 

H  i 

14) 

H  (  3  • 

= 

.  3 1 5 6 2 B  0  0  E  -  0  4 

— 

H  i 

13) 

H  i  4  i 

= 

.  1 3  o  3  S  7  7.QE  t-  0  0 

= 

H  •, 

12) 

H  (  5) 

= 

. 3391 3420E-04 

— 

H  < 

.1. 1 ) 

Hi  6  / 

= 

201 1 3000E+00 

H  i 

1  0  > 

Hi  7) 

=: 

. 33453400E-04 

H( 

9 ) 

H  (  G) 

. 23S48360E+00 

= 

H  ( 

3) 

LOWER  BAND  EDGE 

UPPER  BAND  EDGE 

DESIRED  VALUE 

WEIGHTING 

DEVIATION 

DEV  I ATI UN  IN  DB 


BAND  1 

. ooooooooo 

. 150000000 

. ooooooooo 
1 o. ooooooooo 

. 03992 1 9 1 o 
•L  /  •  975  /  40000 


EXTREMA  FREQUENCIES 

.■  .rmcrn^ ■  i  •*,  ••  e-\ 

•  v  JO  J  l'  •— 1 0  •  i  i  v  /  ( — / 

. 2700000  «  3500000 


*  BAND  2 
. 230000000 
. 27000000 0 
1 . OOOOOOOOO 

1 0 . ooooooooo 

Q C2 /*?  1  Cl  ;m\ 


i  so  oooo 


1C 


. 377343/  . 439S43 / 

I.  ’A‘  .!/  U.-  w  vif  \l/  w  U/  ■.(<  -At  \L>  tb  tv  J.1  \i.-  vl*  il/  ')/  \l/  'j:  »l»  ■.!.  '-v  vk  vl/  \i .•  \i  \b  \U  J.  \V  w  w  »!•  •/  >1/  tl-  J  vl  •- 1  -A:  •A.'  -J 

+  *■  +  +  *•  +  *  *■  k  ft-  ft-  <f-  4-  Jf.  *t-  4-  ft-  ft-  r  4-  4--  4  4  T  T-  4 


BAND  3 
350000000 
5 0 0 0 0 O  0  0  0 
>0 0000000 
>00000000 
>3992  3  9 1 C 


H  4  ■*!  '  V  ' 


! 
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Tabic  CIO  :  Filter  Coefficients  for  BPF2500s  L31 


************************************************* ******** 


*****  I 

MPULSE 

RESPONSE  ***** 

H  ( 

i )  = 

.  ci  7  !  j 3924 0 E — v  o  — 

H  i 

3 1  . 

I-!  <. 

= 

. 9 u 4917 3 0 E— 02  — 

H  < 

30  • 

H  ( 

3  \  ~ 

- . 29628 1 20E-05  * 

hi 

29 ;» 

H  i 

4)  = 

18620890E-01  = 

H  i 

28) 

H  ( 

5  )  = 

-. 97091 370E-06  = 

H  i 

.  ~7  • 

+L  /  / 

H  ( 

6 )  = 

.  1S94U40E-01  = 

H  i 

26) 

H  ( 

7 )  = 

.  118781 30E-05  = 

Hi 

25) 

Hi 

8)  = 

.  64614740E-02  = 

H  i 

’7  Al 

Hi 

9)  = 

.  62477670E— 06  — 

H  i 

/ 

H  ( 1 0 )  » 

- . 64797790E-0 1  = 

H  i 

22 ) 

H  ( 

11)  = 

. 2387930 0 E - 0 5  = 

H  i 

21  > 

H  •» 

12)  = 

. 1 437 1 060E+00  = 

H  ( 

'-.  I  v 

H  ( 1 3 )  = 

-. 195491 10E-05  = 

H  ( 

i  9 ) 

H (1 4 )  = 

-.21 299060E+00  = 

H  i 

1 S ) 

H  i 

.B2123430E-Q6  = 

H  i 

1 7  '' 

H  i 

16)  = 

.  2 4 0 o 4 20 0 E  •+  0 0  * 

i  1  / 

n  \ 

1 6 ) 

BAND 

1  BAND 

2 

r  i  \  i  ir  — 

n  ■ 

LOWE H  3 rt N L) 

EDGE 

. 000000 

000  . 2300000 

0  0 

3 5 0 0 0 0 0 C  • 

UF'F'EF,  3ANL 

EDGE 

. 150000 

OOO  . 2700000 

i  *i 

500000000 

DESIRED  VALUE 

. oOoooo 

00 0  1 . 00 00000 

0  0 

0  O  O  :j  v  O  O  v 

WEIGHTING 

10 

. OOOOOO 

0 0 0  1 0 .  0 0 0 0 0 0 o 

Cj  C) 

i  m 

OOOOOOC;:jO 

DEV I  AT I ON 

. 004 148 

392  . 004 1 433 

004148392 

DEV I AT  I ON 

IN  DB  -47 

. 642350 

0 0 0  —47.6423 5 0 0 

■JO 

-47. 

6 42 3 5 O0 00 

EXTREMA  FREQUENCIES 

. 0 0 0 0 0 O 0  .03125 o 0  . 0 6250 0 0 
. 1 406250  .1500000  .2300000 
. 3500000  . 3578 1 25  . 3782969 
. 469 1 406  . 5000000 


406640. 


*  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  * .  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  X  *  *  *  x  *  *  *  *  * 


a  %j.  »!•  ti.  •>* 


rt 

i 


i 


& 


baa 


FINITE  IMPULSE  RESPONSE  OFIP  • 
LINEAR  PHASE  DIGITAL.  FILTER  DESIGN 
REME2  EXCHANGE  ALGORITHM 


BANDPASS  FILTER 
FILTER  LENGTH 


=  t'O 


***** 

I MPU 

_SE  RESPONSE  ***** 

H 

f  1 ) 

=  - . 3553 1 050E— 06 

= 

H  ( 

£>  •->  } 

H 

(  2) 

=  . 23697640E— 03 

SS 

H( 

O-L  ) 

H 

(  3) 

=  .  13197250£=£5 

= 

H  ( 

r  a  \ 

o  i  ■■ 

H 

(  4) 

=  61 86754 QE -03 

s 

H  ( 

60 ) 

H 

s'  5) 

=  20625860E-05 

Sr 

H  ( 

5°  ' 

H 

(  6 ) 

=  .  62493 140E-03 

= 

H  t 

53 ) 

H 

(  7 ) 

=  .  1 048 1 090E—05 

= 

H  ( 

57) 

H 

O 

=  .664 0 0 0 4  0 E - 0 3 

s= 

H  • 

56) 

H 

(  9) 

=  *  278387 20E-05 

= 

H  ( 

55) 

H 

( t  o  > 

=  - . 354S8550E-02 

H  ( 

54  > 

H 

(ll) 

=  793946S0E-05 

= 

H  ■' 

rr-r\ 

wJ  0 

H 

(  1 2 ) 

=  . 63765240E-02 

=r 

H  ( 

52) 

H 

(  13) 

=  .  1 0  o  0  6 1 5  V  E — O  4 

= 

H  ( 

53  ) 

H 

<  1  4  ) 

=  53608550E-02 

3= 

H  ( 

50) 

H 

( 1 5  > 

=  -- .  4509 1 6S0E-05 

= 

l-l  * 

49) 

H 

(16) 

=  30799540E-02 

— 

l-i  ( 

48) 

H 

(17) 

*  - . S2562560E-05 

= 

H  ( 

47 ) 

H 

(13) 

=  . 18017490E-01 

= 

H  < 

46) 

H 

(19) 

=  .21 076480E— 04 

= 

H  ( 

45) 

H 

<  20 ) 

-  -. 30320 170E-01 

= 

HO 

44 

H 

(21 > 

=  235912 10E--04 

= 

H  ( 

43) 

H 

/  \ 

\  46 

*  . 26S72950E-01 

= 

H  ( 

42) 

H 

(23) 

=  .  10135640E-04 

= 

H  ( 

41  > 

H 

( 24 ) 

=  .  66S50530E-02 

= 

H  ( 

40) 

H 

(25) 

=  .  13716950E-04 

= 

H  ( 

39) 

H 

( 26 ) 

=  -.  70965640E-01 

= 

!-!  < 

TO 

...<  U.I  •' 

H 

'  27  > 

=  333 1 2280E— 04 

= 

(..!  ‘ 

.l’1  * 

I-: 

( 28  ) 

-  _  ‘  5;  m  )S(  J.-'**"  <E ■+  cj( ! 

s= 

H  t 

•  -O  ' 

H 

>;  29 ) 

“  .  4 4 6  /  4 6 O  ti  0  4 

=•* 

H  < 

N 

i  7  *) 

~  — . 2 1 5S7720E+U0 

H 

ft 

H 

(31  ) 

=  - . 3 4526260E-04 

= 

H  < 

7"^  I 

H 

( 32 ) 

=  . 24 14871 OE+OO 

5= 

H  * 

“3  r* 

1 


LOWER- 

UPPER 


BAND  EDGE 
BAND  EDGE 
DESIRED  VALUE 
WE I GHT I NG 
DEVIATION 
DEVIATION  IN  DB 


BAND 

. 000000000 
. 150000000 
. 000000000 
1 . 000000000 
-  000060967 
-84.  297950000- 


BAND 


1 


270000000 
OOOOOOO'X) 

ooooooooo 
00006096“' 
-84.  29795000*"' 


EXTREMA  FREQUENCIES 
„ 0000000 
. 076 1 " 1 9 
. 1406250 
•  2397 2  56 
. 3500000 
39492 1 9 
. 469 3 406 


0 1 56250 
090S203 
1 4746  09 


r-  ■—>  r>  '**■ .  <-•>  - 


4095703 

4847656 


0302734 
1044922 
1 50  <0000 
2602735 
3597656 


0458 
1  1  5 1 


5000000 


><  J'  K.  > 


•  i  )OU 


.x  u 


OOOOc  ’ 

,* 


1  6HH 
'"■“'■■DO 


.  —  04*- 


Table  Cl  2  :  Filter  Coefficients  for  BPF2500  L9 


*********************************************************** 

FINITE  IMPULSE  RESPONSE  (FIR) 

LINEAR  PHASE  DIGITAL  FILTER  DESIGN 
REMEZ  EXCHANGE  ALGORITHM 

BANDPASS  FILTER 

FILTER  LENGTH  -  9 

*****  IMPULSE  RESPONSE  ***** 


H  ( 

1) 

m 

. 39850840E— 01 

m 

H  ( 

9) 

H  < 

2) 

m 

. 83494020E— 08 

m 

H  ( 

8> 

H  ( 

3) 

m 

- . 29925420E+00 

m 

H  ( 

7) 

H  ( 

4) 

m 

. 14653130E-07 

m 

H  ( 

6) 

H  ( 

5) 

m 

. 38029830E+00 

m 

H( 

5) 

LOWER  BAND  EDGE 
UPPER  BAND  EDGE 
DESIRED  VALUE 
WEIGHTING 
DEVIATION 
DEVIATION  IN  DB 


BAND  1 
. OOOOOOOOO 
. 100000000 
. OOOOOOOOO 
1 . OOOOOOOOO 
. 098508380 
20. 130310000 


BAND  2 
. 200000000 
. 300000000 
1 . OOOOOOOOO 
1 . OOOOOOOOO 
. 098508380 
-20. 13051000c 


BAND  3 
. 400000000 
. 500000000 
. OOOOOOOOO 
1 . OOOOOOOOO 
. 098508380 
-20. 130510000 


EXTREMA  FREQUENCIES 

.0000000  .1000000  .2000000  .2500000  .3000000 

. 4000000 


*********************************************************** 
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Tabic  C13  :  Filter  Coefficients  for  BPF2500  L15 


********************************************************** 

FINITE  IMPULSE  RESPONSE  (FIR) 
LINEAR  PHASE  DIGITAL  FILTER  DESIGN 
REMEZ  EXCHANGE  ALGORITHM 

BANDPASS  FILTER 


FILTER  LENGTH 

*****  IMPULSE 

H  (  1)  - 
H  <  2)  - 
H (  3)  « 

H  (  4)  ■ 

H (  5)  - 
H <  6)  • 
H  <  7)  ■ 

H  (  8)  * 


-  15 

RESPONSE  ***** 
. 57887090E— 08 
. 38867450E— 01 
6B377540E— 08 
. 71470 100E-01 
1 498 121 OE— 07 
288B6740E+00 
4381 1290E— OB 
.41 0735 1 OE+OO 


LOWER  BAND  EDGE 
UPPER  BAND  EDGE 
DESIRED  VALUE 
WEIGHTING 
DEVIATION 
DEVIATION  IN  DB 


BAND  1 
. 000000000 
. lOOOOOOOO 
. 000000000 
1 . ooooooooo 

. 053675240 
>25.404490000 


BAND  2 
. 200000000 
. 300000000 
1 . OOOOOOOOO 
1 . OOOOOOOOO 
. 053675240 
■25. 404490000 


BAND  3 
. 400000000 
. 500000000 
. OOOOOOOOO 
1 . OOOOOOOOO 
. 053675240 
-25. 404490000 


EXTREMA  FREQUENCIES 

.0000000  .0666667  .1000000  .2000000  .2500000 

. 3000000  . 4000000  . 4333333  . 5000000 

********************************************************* 


o| 

>3 1 


i 

H  ( 

15) 

V 

V 

H( 

14) 

H  ( 

13) 

s 

H  ( 

12) 

H  ( 

11) 

*> 

H  ( 

10) 

v: 

H  ( 

9) 

V 

H  < 

8) 

a 
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Table  C14  :  Filter  Coefficients  for  BPF2500  L31 


****************  ****************************************  *  * 

FINITE  IMPULSE  RESPONSE  (FIFO 
LINEAR  PHASE  DIGITAL  FILTER  DESIGN 
REMEZ  EXCHANGE  ALGORITHM 

BANDPASS  FILTER 

FILTER  LENGTH  =  31 


IMPULSE 

RESPONSE  ***** 

H  ( 

1  ) 

-. 421S3190E-05 

— 

Hi 

31  ) 

H  ( 

2) 

= 

- . 36066390E-02 

= 

H  ( 

30) 

h  •; 

3) 

= 

. 227 1 1 920E-05 

= 

H  ( 

29) 

H  ( 

4) 

= 

. 1 1 444820E-0 1 

— 

H  ( 

28) 

H  ( 

5) 

= 

. 537566 10E-05 

= 

H  ( 

27) 

H< 

6) 

= 

. 21266940E-02 

= 

H  ( 

26) 

H  ( 

7) 

= 

— . 66096S60E—05 

= 

H  ( 

25) 

H  ( 

8) 

= 

-. 42699350E-01 

* 

H  ( 

24) 

H  ( 

9) 

= 

. 6961 1780E— 06 

= 

H  ( 

23) 

H  ( 1 0 ) 

= 

. 42655 130E-01 

= 

H  ( 

W  ) 

H  <  1 1 ) 

s 

. 422375 10E-05 

as 

H  ( 

21) 

H  ( 1 

2) 

S 

. 84 1 09580E— 0 1 

'as 

H  ( 

20) 

H  < 1 

3> 

= 

-. 29676600E-05 

= 

H  ( 

19) 

H  <  1 4 ) 

as 

-.291 1 7750E+00 

= 

H  ( 

18) 

H  ( 1 

5) 

s 

. 1231 3500E-05 

= 

H  ( 

17) 

H  ( 1  6 ) 

= 

. 396 1 0960E+00 

= 

H  ( 

16) 

LOWER  BAND  EDGE 
UPPER  BAND  EDGE 
DESIRED  VALUE 
WEIGHTING 
DEVIATION 
DEVIATION  IN  DB 


BAND  1 

0  0  O 0 0 0 0 o  o 
. 1 oooooooo 
. 000000000 
1 0 . ooooooooo 

. 00 18151 04 
-54.82191 0000 


BAND  2 

. 200000000 
. 300000000 
1 . OOOOOOOOO 
1 0 . OOOOOOOOO 
. 00 18151 04 
■54 .8219 1 0000 


BAND  3 

u  4  *  J  0  Ci  Ci 0 0 0 C- 
. 500000000 
. OOOOOOOOO 

1 0 » ooooooooo 

. 00 1815 1 04 
-54 .82191 0000 


EXTREMA  FREQUENCIES 

.0000000  .0390625  .0703125  .0917969  .100000*; 

.  2000000  .  2078 1 25  .  2253906  .  25078 1 2  .  2”’ 42 1  S'" 

.  29 1 7969  .  3000000  .  4000000  .  4078 125  .4  3 1  250C. 

. 4625000  . 5000000 

* * * * *  * * * % 1 1 %  * 1 1 * *  *  # * ********************** * * * * *  * * * * * * * * * 


Table  C15  :  Filter  Coefficients  for  HPF2500  L31 


C****#*» ************************** ******************* 


FINITE  IMPULSE  RESPONSE  (FIR) 
LINEAR  PHASE  DIBITAL  FILTER  DESIGN 
REMEZ  EXCHANGE  ALGORITHM 


FILTER  LENGTH  = 


X  X  X  X  X 


IMPU! 

LSE 

RESPONSE  ***** 

H  ( 

1  / 

s 

. 2 1 0 428 0 0 E - 0 2 

= 

H  ( 

315 

H 

2; 

=r 

1 9346590E— 05 

= 

H  ( 

30; 

H  < 

3) 

=5 

-.  464d>5720E— 02 

= 

H  ( 

29 ) 

H  ( 

n 

/ 

as 

.  7303379 0 E  -- 0 6 

= 

H  ( 

23) 

H  ( 

5) 

55 

. 942 1 9370E-02 

= 

H  ( 

27 5 

H  < 

6) 

ss 

— . 23609320E— 05 

= 

H< 

26) 

H  ( 

7) 

s 

-. 1 7204930E— 0 1 

s 

H  ( 

25) 

H  ( 

8) 

= 

. 202 1 2450E-05 

= 

H  ( 

24) 

H  ( 

9) 

4= 

. 29778450E-01 

= 

H  < 

23) 

H  UO) 

= 

— , 5 1 4595 1 OE— 05 

= 

H< 

22) 

H  ( 1 1 ) 

= 

-.51521 BOOE-O 1 

= 

H  < 

21  ) 

H  ( 12) 

= 

.41 184310E-05 

= 

H  ( 

20) 

H  ( 1 3 ) 

.9841 637 0E-01 

= 

H  ( 

195 

H  ( 

14) 

= 

-. 47544330E-05 

= 

H< 

IS) 

H  ( 15) 

= 

3 1 56787 OE+OO 

- 

H  ( 

17) 

H  ( 

1  6) 

3= 

.  5 0 0 0 0 4 3 0 E  -i-  00 

= 

H  ( 

1  w  • 

LOWER  BAND  EDGE- 
UPPER  BAND  EDGE 
DESIRED  VALUE 
WEIGHTING 
DEVIATION 
DEVIATION  IN  DB 


BAND  1 

. ooooooooo 
. 200000000 
. ooooooooo 
1 o . ooooooooo 

. 001349337 
-57 . 397530000 


our.!) 

. 300000000 
. 500000000 
1 . OOOOOOOOO 
10. OOOOOOOOO 
.001349337 


-57.397530000 


EXTREMA  FREQUENCIES 

. 03 1 2500  . 0625000  . 09 1 7969  .121 0933 
.1738281  .1933594  .2000000  . 3000000 
.3253906  .3507813  .3781250  .4074219 
. 46992 1 9  . 5000000 


1 4S4375 

-T  .••••. 7n  ■*  . er 

f  u  i  ^ 

,  45&fe / 1 ~ 


************************************************* * * * * * * * 
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a 

$ 

f 

s 

a 
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Table  C16  :  Filter  Coefficients  for  HPF3500  L9 


a*****######***##****##**##**#*##****#******#*#*#**#*****:! 


FINITE  IMPULSE  RESPONSE  (FIR) 
LINEAR  PHASE  DIGITAL  FILTER  DESIGh 
REMEZ  EXCHANGE  ALGORITHM 


BANDPASS  FILTER 
FILTER  LENGTH  =  9 


*****  IMPULSE  RESPONSE  **:**: 


H  ( 

1  > 

= 

-. 473S5250E-01 

= 

H  ( 

9) 

H  < 

2) 

= 

51096310E-01 

= 

H  < 

S) 

H  ( 

3 ) 

- 

. 1 5307250E+00 

= 

H  ( 

7) 

H  ( 

4) 

= 

-. 25925620E+00 

= 

H  ( 

O  ) 

H  ( 

5  > 

= 

.  3 0 4 3 6 9 60 E + 0  0 

sc 

H  ( 

5  \ 

LOWER  BAND  EDGE 
UPPER  BAND  EDGE 
DESIRED  VALUE 
WEIGHTING 
DEVIATION 
DEVIATION  IN  DB 


BAND  1 

. ooooooooo 

. 300000000 

. ooooooooo 

1 . 300000000 
. 104961000 
■19.579420000 


BAND  2 
. 400000000 
. 500000000 
1 . OOOOOOOOO 

1 . ooooooooo 

. 136449300 
-17. 300560000 


EXTREMA  FREQUENCIES 

.0000000  .12 

. 5000000 


X  JL 


3000000 


.  4  00' 


t ************* * ****** * *  *  *  * * ************** * * * * *  *  *  *  * **** *  * ; 


VV 


faS4J!fiS5Sfl£53 


awes 


Lrj 

fef 

L*  » 


* 

k 

s 
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Table  C17  :  Filter  Coefficients  for  HPF3500  L15 


I-********************************#**************  ********** 

FINITE  IMPULSE  RESPONSE  (FIR) 
LINEAR  PHASE  DIGITAL  FILTER  DESIGN 
REMEZ  EXCHANGE  ALGORITHM 

BANDPASS  FILTER 

FILTER  LENGTH  =  15 

*****  IMPULSE  RESPONSE  ***** 


H  ( 

1  \ 

-. 1 32848 10E-01 

= 

H  ! 

1 5 ) 

H  \ 

: 

= 

2272337 OE— 0  3 

S 

H  < 

<  n 

H  ( 

•.  i 

SK 

.  4 4 7 6 4 2 4 0 E — 0 1 

= 

H  <. 

3  7. 

H  ! 

4) 

= 

-.38 0 6 7 3 0 0 E - 0 1 

I-!  ( 

1  2  ! 

hi  < 

5 ) 

= 

—  .27  oS  7  c<80E — 0 1 

H( 

.  3.  ) 

H  ■ 

o  > 

. 141 937 1 OE+OO 

8S 

H  ( 

1 0 ) 

H  < 

— r  v. 

-. 25439 1B0E+00 

sr 

H< 

9) 

H  ( 

2) 

= 

. 30128990E+00 

= 

H  ( 

8) 

BAND  1 

LONER  BAND  EDGE  . OOOOOOOOO 
UPPER  BAND  EDGE  .300000000 
DESIRED  VALUE  •  .OOOOOOOOO 
WE I GHT I NG  1 . OOOOOOOOO 
DEV I AT I  ON  . 0364273 1 0 
DEVIATION  IN  DB  -28.771430000 


BAND  2 

.  400000000 
. 500000000 
1 . OOOOOOOOO 
1 .  OOOOOOOOO 
. 0364273 1 0 


-28 . 77 1 430000 


EXTREMA  FREQUENCIES 

. 0 oO  C '  0 0  •  •  . 0 7 0 3  3 25 
. 3O0000C  • 4000000 

*  ******  *  *  *  *  *  *  *  *  ****  *  *  *  *  *  *********  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  S'  *  *  *  *  N: : 


. 1406250  .2109375 

. 4312500  . 5000000 


Table  C18  :  Filter  Coefficients  for  HPF3500  L31 


i 


****$*#***#*#**##**#* 


FINITE  IMPULSE  RESPONSE  (FIR) 
LINEAR  PHASE  DIGITAL  FILTER  DESIGN 
REME2  EXCHANGE  ALGORITHM 


BANDPASS  FILTER 
FILTER  LENGTH  =  31 


t**rt  IMPULSE  RESPONSE  ***** 


H  (  1  • 

= 

1 8 9 6 8 2 4 0 E — 0  2 

= 

H  ( 

•T  +  . 

1  • 

h  •: 

2 ) 

= 

. 1 6820350E-02 

= 

H  ■. 

30 ) 

H  ( 

3) 

s 

. 13447050E— 02 

= 

H  ? 

2:9  '> 

H  ( 

4  :• 

= 

62657850E-02 

5= 

H  f- 

nr.  % 
X.  u.'  • 

H(  5) 

= 

. 76874290E— 02 

= 

H 

27> 

H  (  6> 

= 

345 1 6920E-03 

S 

H  < 

26> 

H  (  7) 

1 3507 420E-0 1 

= 

H  ( 

25) 

H  (  8) 

= 

. 21473700E-01 

= 

H  < 

24 ) 

H  (  ?) 

= 

9476359OE-02 

H  ( 

'-.‘7 

H  ( 10) 

= 

22502630E-0 1 

= 

H  ( 

. 

H  (  1 1  ) 

s 

. 5 1 378460E— 0 1 

= 

H  ( 

21 ) 

H  <  12) 

= 

-.4121 8850E— 0 1 

= 

H  ( 

20) 

H  (  1 3 ) 

= 

- . 29878980E-0 1 

= 

H  ( 

19) 

H  (14) 

— 

. 1 46 1 57 1 OE+OO 

= 

H  ( 

18) 

H  ( 1 5 

= 

-. 25565 100E+00 

se 

H< 

17) 

H  1 16 .) 

= 

. 30051 6° OE+OO 

= 

H  ( 

16) 

LOWER  BAND  EDGE 
UPPER  BAND  EDGE 
DESIRED  VALUE 
WEIGHT  INC 
DEVIATION 
DEVIATION  IN  DB 


BAND  1 

.  0 0 0 0 O 0 0 0 0 
. 300000000 


. 000000000 


1 O . 000000000 
.  00 1 522 1 6. 1 
-56 . 350730000 


BAND  2 

4  0  (")  (")  0  f)  O  O  Q 

.  500000000 
1 . 000000000 • 
1 0 . 000000000 
. OC 1 522 1 6 1 


EXTREMA  FREQUENCIES 


0000000 

. 0312500 

. 0644531 

. 095703 1 

1582031 

. 1894531 

. 220703 1 

. 2480469 

29296SS 

. 3000000 

.  4000000 

. 4078125 

4566406 

. 5000000 

******  * *  *  * * * *  * * *  *  *  *  * * *  * *  * *  * *******  *  * *  *  * *  *  *  *  * * *  *  *  * 


**  * 


a.  a.  *».-  \i. 
•r  .-T  •'!' 
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Table  Cl  9  .  Filter  Coefficients  for  HPF3500  L63 


FINITE  IMPULSE  RESPONSE  (FIR- 
LINEAR  PHASE  DIGITAL  FILTER  DESIGN 
REMEZ  EXCHANGE  ALGORITHM 


BANDPASS  FILTER 
FILLER  LENGTH  =  63 


*A  IMPU 

r-'C 

woe;,  r,  c 

SPONSE  ***** 

H  (  1  ) 

s 

1 1 472290E— 04 

= 

H  < 

o3 ) 

H  <  2) 

=■  — . 

398077 0 0 E — 0 5 

= 

H  ( 

62) 

H  (  3) 

as  —  * 

40 0 9 8 0 0 0 E - 0 4 

= 

H  ( 

6 1  / 

H  (  4  > 

ss 

94604C 10E-04 

as 

H  ( 

60  > 

H  t  5  > 

=  —  * 

6  3  302900E-04 

= 

H  ( 

cr  rz>  • 

H  (  6  J 

as  — . 

1 35883 10E-03 

= 

H  ( 

58) 

H  (  7 ) 

= 

334  36840E-03 

= 

l.J  f 

rr“7  ' 

H  i  S'- 

=r  -  „ 

35313870E-03 

= 

:-i  ( 

‘7,-. 

H  i  ?  : 

=  -  „ 

21528440E-03 

as 

h  < 

err  \ 

H  (  1 0 ) 

= 

1 0483380E —02 

= 

l-i  ( 

Zi  ? 

i-l  <  J  3  • 

—  ~  „ 

1 2 6 0 7 5 4 0 E — 0  2 

= 

H  ( 

cr~r  \ 

h  >; :  2 1 

- 

625 1 8 08 OE- 0*1 

= 

{-I  ( 

vJ--  i 

h  ;  i  3 ) 

21223730E-02 

= 

H  ( 

53. ) 

H  ( 1 4 : 

-- 

333221 40E-02 

= 

H  ( 

50) 

K  1 5  ) 

= 

1 4724660E-02 

= 

H  < 

49) 

H  '  1  6  ) 

7  7  7  7  (*)  7  7  f")  rr  — .  fm)  7 

= 

H  < 

40 

H  (  1  7  , 

= 

7 0 7 7 2 9 3 0 E — 0  2 

= 

H  ( 

47) 

H  (  1 S 

= 

52939430E-02 

= 

H  ( 

46) 

H  f  1 9  > 

3 1 994880E— 02 

H  < 

45) 

F  • 20  > 

SZ  ~  . 

1262433 0E-03 

= 

H  ( 

44 ' 

H  ( 2 1  > 

= 

1 3296960E— V 1 

= 

PN 

43) 

H  <  22 > 

—  —  . 

2 1 857700E-03 

= 

H  < 

42) 

H  (23) 

= 

1941 2380E— 0 1 

= 

H  ( 

41  ) 

H  (24) 

= 

231201 10E-01 

= 

H  ( 

40) 

H  (25) 

=S  -  # 

1 1 433360E— 01 

=. 

H( 

3?) 

H  < 26) 

=  — „ 

26 1 4S700E-0 1 

= 

H  ( 

38) 

u  ( 7  7 

= 

5 6 6 9 0 2 6 OE — 0  3 

H  ' 

7*7  i 

H  :  ?S ) 

a= 

4 3 6 7 9 2 0 0 E — 0  3 

= 

I-!  < 

7  6 

U  ;  7C  -f 

=:  -  , 

31 1 56300E— 0 3 

= 

H  ( 

-»cr 

H  *■  30 ) 

— 

1 4S49800E+00 

= 

H  < 

“T  ft  \ 
*T  * 

H  •'  1  ) 

- 

25652060E400 

5S 

H  \ 

77;  ■ 

H  •  32 ) 

- 

3 0 0 3 2 0 8 0 E  4  0  0 

= 

H  ( 

LOWER  BAND  EDGE 
UPPEF  BAND  EDGE 
DESIRED  VALUE 
WEIGHTING 
DEVIATION 


BAND  1 
.  (X > 00 0 0 0 0 0 
. 300000000 
. 000000000 
1 000000000 
-  00000804? 


BAND  2 

. 400000000 
, 500000000 
1 . 000000000 
1 . 000000000 
. 00000804? 


DEV I  AT  I  ON  I N  DB  - 1 0 1 . 884900000  -  1 0 1 . 884900000 


EXTREMA  FREQUENCIES 

.  OOOOOOO  .  0 1. 56250 
.0782250  . 0937500 
.1562500  .1713750 
.2314453  .24511?2 
.  29 1  '•'922  .2978516 
. 40878? 1  . 4 1 75781 
.4703125  .4849609 


03 1 2500 
1  OfyiTTgO 
1865234 
2587891 
3000000 
4292969 
5000000 


. 0468750 
. 1 250000 

.  U  ,.L  ]  *t  O  *+ 

27  «  47;^.^ 
.  400001)0 
. 44 1 ?922 


m  rrtr*  _  n  -i 
.  -t  x.j  u  r  n-r  a 


mwmm 


mwwm: 


I 


I 


i 

g 


& 

I 


1  _«g  ,1 1 


FILE:  MAKEOATA  FULL  A1  FTD  COMPUTING  FACILITY 


program  MAKEOATA  FILE  *) 

•) 

This  program  allows  tha  usar  to  maka  a  data  flla  that  contains  *) 
a  signal  with  many  different  charactar istics  possibia  or  soma  *) 
combination  of  charactar istics.  Thasa  posslbl 1 itlas  ara  llstad  *) 
in  tha  first  part  of  tha  main  program  and  include  various  tlma  *) 
wavaforms  and  two  diffarant  nolsa  distributions.  Tha  program  Is  *) 
daslgnad  to  glva  tha  usar  maximum  flaxlblllty  In  selecting  tha  *) 
das  1 rad  signal  charactar istics.  *) 


program  MAKEDATA_F I LE ( INPUT , OUTPUT ) ; 


PI  >  3.1415926535898: 
MAX_NUM_SAMPLES  •  16384; 


CH  :  CHAR ; 

DATE. 

TIME  :  ALFA; 

DATA  SAMPLES  :  TEXT; 

SEE01 , 

SEED2. 

ANS. 

I. 

J  :  INTEGER; 

A. 

U. 

R. 

0. 

X. 

DUMMY  NUM, 

INITIAL  VAL, 

SIGMA  SO. 

LEVEL? 

PROBIMP. 

IMPULSE  VALUE. 

STEP  OFFSET. 

NOISE  AMPLITUDE  :  REAL; 

FILE  CHOICE. 

FREQ? 

CHOICE. 

START  FREQ. 

STOP  FREQ, 

START  SAMP. 

STOP  SAMP. 

NUMSAMP . 

NUMSTEPS . 

SIGN  :  INTEGER; 

DONE  :  BOOLEAN; 

SIGNAL  :  ar ray( . 1 . .MAX_NUM_SAMPLES. )  of  REAL; 
XINCLUOE  CMS 


(•  function  SGN 

(• 

(•  This  function  returns  tha  sign  of  tha  TEST_VALUE. 


function  SGN  (  TEST_VALUE  :  REAL  )  :  INTEGER; 
bag  in 


’■V 


& 


I 


I 


B 


$ 

$ 

I 

r 

D. 


FILE:  MAKEDATA  FULL 


A 1  FTO  COMPUTING  FACILITY 


132 


if  TEST  VALUE  <  O  then 
SGN  -1; 

and; 


begin  (•  of  MAKEOATA_FILE  •) 

torn 1n( INPUT); 
tormout ( OUTPUT ) ; 

CMS ( ' CLRSCRN ' , I ) ; 

wr1ta1n( 'Would  you  Ilka  to  1)  add  to  tha  existing  signal'); 
wrltaln('  2)  add  to  tha  ideal  signal'); 

writa1n('  or  also)  naka  a  now  signal'); 

raadln(ANS); 

If  (ANS  »  1)  or  (ANS  >2)  than 
bagin 

casa  ANS  of 

1  :  rasat(OATA  SAMPLES. 'nana-SIGNAL.POATA. *') ; 

2  :  reset ( DAT A~SAMPLES , ' name* IDEAL . POATA .  •  '  ) ; 
and; 

raad1n(0ATA  SAMPLES. I .DUMMY  NUM); 
while  not (eof( DATA  SAMPLES)?  do 

raadln(0ATA  SAMPLES. I . SIGNAL( . I . )  )  : 

NUMSAMP 

wrltaln( 'Thara  ara '. NUMSAMP : 6 , '  sanplas  of  data'); 

and 

also 

rapaat 

wrlta1n('How  many  sanplas  to  taka  (1024  racommandad ) ? ' ) ; 
raadln( NUMSAMP); 

If  (NUMSAMP  >  MAX_NUM_SAMPLES )  or  (NUMSAMP  <  1)  than 

wr1ta1n('The  number  of  sanplas  must  ba  between  1  and  4096'); 
until  (NUMSAMP  >  0)  and  (NUMSAMP  <•  MAX_NUM_SAMPLES ) : 
datet 1 na( DATE. TIME ) ;  (*  ganorata  random  sead  for  •) 

roadstr ( str (TIME ) . I : 2.CH, SEE01 : 2.CH.SEED2) ;  <•  random  number  •) 

I  :■  I *SEE01 *SEE02+I+SEED1+SEED2 ;  (•  generator  •) 

U  : •  random( I ) ; 
wrlteln( I ,U) ; 
wri tain; 

wr1ta1n( 'BUILD  A  SIGNAL'); 

rapaat 


wrlta1n( 'What  would  you  like  in  it?'); 

wri tel n( ' 

D 

Constant  level'); 

wri teln( ' 

2) 

Monotone  ramps'); 

writeln( ' 

3) 

Sinusoids' ) ; 

wrlteln( ' 

4) 

Stops' ) ; 

wri ta1n( ' 

S) 

Gaussian  noise'); 

wri tel n( ' 

6) 

Leplacian  noise'); 

wrlte1n( ' 

T) 

HR  filter  the  signal' 

writaln( ' 
raadl n( CHOICE); 
casa  CHOICE  of 

8) 

That' 's  all'); 

1  ;  begin 

wrlta1n( 'What  shall  the  constant  lavsl  bo?'); 
raadl n( LEVEL ) ; 
if  LEVEL  <>  0  than 

for  I  :•  1  to  NUMSAMP  do 
SIGNAL(.I.)  :■  LEVEL; 
and; 

2  :  bagin 

writa1n( 'What  is  the  starting  samp la  for  tha  ramp?'); 

raadl n( START_SAMP ) ; 

writeln( ' And~the  ending  sample?'); 

raadl  n(  ST0P_SAMP ) ; 

wrlta1n( 'To“what  relative  value  should  the  ramp  r iso(f al 1 )?' ) ; 
raadl n< LEVEL); 

INITIAL  VAL  :■  SIGNAL( . START  SAMP.); 
for  I  START  SAMP  to  STOP  SAMP  do 

SIGNAL(.I.)  :■  INITIAL  VAL  ♦  LEV£L/obs( STOP  SAMP-START  SAMP)* 


FILE:  MAKEDATA  FULL 


A 1  FTD  COMPUTING  FACILITY 


( I  -  ST  ART_SAMP ) ; 


3  :  begin 

writaln! *  Input  frequency  of  slnuelod  in  hz(0  if  group  desired)  ); 
raadln(FREQ); 

If  FREQ  <>  0  then 
begin 

wr 1  tel n( 'starting  sample(0  If  ell  samples)?'); 
read 1 n ( ST ART_SAMP ) ; 

If  START_SAMP  ■  0  then 
begin 

START  SAMP  :•  1; 

STOP_SAMP  :•  NUMSAMP; 

end 

else 

begin 

writaln! 'ending  sample?'); 
read1n(ST0P_SAMP); 
end; 

for  I  :■  START  SAMP  to  STOP_SAMP  do 

SIGNAL( . I . )  : ■  SIGNAL ( . I ■ )  ♦  S ln( 2*PI *FREQ*I/NUMSAMP ) ; 

end 

else 

begin 

wr 1  tel n< 'Whet  la  the  starting  f nequency? ' ) ; 
read! n( START  FREQ): 

wr 1  tel n( 'Whet  is  the  last  f requency? ' ) ; 
read! n( STOP  FREO); 

for  FREO  :■  START  FREO  to  STOP_PREO  do 
for  I  :«  1  to  NUMSAMP  do 

SIGNALS  I.)  :•  SIGNAL(.I.)  ♦  S1n( 2*PI •FREQ*I/NUMSAMP ) ; 

end; 

end; 

4  :  begin 

wrlteln('How  many  steps  in  the  signal?'); 
readln(NUMSTEPS); 
for  U  :•  1  to  NUMSTEPS  do 
begin 

writeln( 'Input  the  starting  sample  for  step'.U); 

reed! n( START_SAMP ) ; 

wr 1 teln( 'And  the  last  sample?'); 

readln(STOP_SAMP); 

wr 1 teln( 'And  the  offset?'); 

read! n( STEP  OFFSET); 

for  I  START_SAMP  to  STOP_SAMP  do 

SIGNAL( . I . )”: •  SIGNAL(.I.)  ♦  STEP_OFFSET; 

end; 

end; 

5  :  begin 

wrlteln( 'What  shall  the  noise  sigma  squared  be  (gaussian  dlst.)?  ;; 
writaln! 'Sigma  squared?'); 
read1n(SIGMA_SQ) ; 
for  I  :■  1  to  NUMSAMP  do 
begin 

U  :■  random(O); 

R  :•  sqrt<-2.0»  SIGMA_SQ  *  1n(U)): 

U  :■  random(O); 

0  :■  2.0  •  PI  •  U; 

X  :■  R  *  cos(O); 

SIGNAL( . I . )  :•  SIGNAL! . I . )  ♦  X; 
end; 


6  :  begin 

writaln! 'What  shall  the  probability  of  an  impulse  be?'); 
readln(PROBIMP) ; 

If  PROS  IMP  >  0  then 
begin 

writeln('Do  you  want  1)  all  positive  impulses  or  2)  both  negative 


P, 

$ 


B 
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'  and  posltlva  impulsas. ' ): 
raadln(SIGN); 

If  SIGN  -  2  than 
SIGN  -1; 

wr 1 ta)n( 'What  shall  tha  impulsa  standard  davlation  ba?'); 
raad1n( IMPULSE  VALUE); 
for  J  :•  1  to  NUMSAMP  do 
bag  In 

U  :■  randon(O); 

If  U  <■  PROBIMP  than 
bagln 

U  :«  random(O); 

X  :•  -(IMPULSE  VALUE/sqrt ( 2) )"1n( 1-U) ; 

SIGNAL (.J.)  :-_SlGNAL( . J. )  +  X  •  SGN( randon(O)  + 
(O.S  •  SIGN)); 

and; 


and; 


and; 


7  :  bagln 

writaln( 'What  is  tha  paramatar  ALPHA?'); 
raadln(A); 

for  I  :■  NUMSAMP  downto  2  do 

SIGNAL(.I.)  :•  A*SIGNAL( .1-1.)  ♦  SIGNAL(.I.); 
and; 


8 

DONE 


:■  trua; 


OTHERWISE  (* 

wr1taln( 'That  Is  not  an  option'); 


and; 

CMS( 'CLRSCRN'  .1); 
unt 1 1  DONE ; 


do  nothing 

of  caaa  statamant 


wrlta1n('ls  this  an  1)  idaal  signal  or  a  2)  signal  to  ba  flltarad?'); 
raad1n(FlLE  CHOICE); 

If  FILE  CHOICE  ■  1  than 

rawr~ta(OATA  SAMPLES, 'nama*IOEAL .POATA . * ' ) 


Mi 


rawrl tatDATA  SAMPLES. ' nama»SIGNAL . POATA . * ' ); 
(•  wr 1 tal n( DAT A^SAMPLES ,0. NUMSAMP/ 1000.0: 15:8);  •) 
for  I  :-  1  to  NUMSAMP  do 

wr1ta1n(DATA  SAMPLES, I :4,SIGNAL( .1 . ) : 15 : 8 ) : 


and; 


C 


Of  MAKEDATA  FILE 


; 


VI 

vj 

vj 

>1 

\1 

V] 
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program  FREQ 

This  program  takas  tha  diacrata  Four  tar  transform  (OFT)  of  a 
signal.  Tha  signal  may  ba  up  to  16384  elements  long,  but  that 
langth  is  not  racommandad  dua  to  tha  langth  of  tima  to 
ealculata.  Tha  OFT  is  Implamantad  using  a  decimatlon-in- 
fraquancy  FFT  algorithm. 


program  FREQ  (INPUT, OUTPUT); 
%INCLUDE  CMS 


PI  -  3. 1415926535898; 
MAX_ARRAY_SIZE  -  16384; 


CHOICE. 

DUMMY, 

I. 

0. 

M. 

N. 

L. 

NV2, 

NM1. 

IP. 

LE. 

K  :  INTEGER; 
LARGEST . 

TEST. 

LEI. 

UR. 

UI. 

WR, 

WI. 

TR. 

TI. 

TMR. 

TMI  :  REAL; 
MAG. 

PHASE. 

XR. 

XI  :  array! . ‘ 

FILEA  :  TEXT; 


array! . 1 . .MAX  ARRAY  SIZE.)  of  REAL; 
TEXT; 


function  RAISE 

This  function  raises  A  to  an  integer  power  N.  This  function  is 
not  tntnnsical  ly  available  in  the  Pascal  compiler 


function  RAISE  (  A.N  :  INTEGER  )  INTEGER; 


TEMP. 

I  :  INTEGER; 


TEMP  •  1; 

for  I  *  1  to  N  do 

TEMP  •  TEMP  •  A: 


f.v./.v 
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RAISE  :•  TEMP; 


function  LOG 


Thla  function  laplementa  tho  baao  10  logarithm  of  X.  Tha  Pascal 
compilar  only  givaa  an  intrinalc  natural  logarithm.  Tha  baas  10 
logarithm  la  cal cuts tad  by 


log  x  •  In  x 
in  10 


•) 

•) 

•) 

•) 

•) 

•) 

•) 

•) 

•) 

•) 

•) 


function  LOG  (  X  :  REAL  )  :  REAL; 


If  X 
X 


LOG 


■0.0  than 
:•  O  00000001; 

•  1n(x)/1n( 10.0); 


begin 


(  •  of  FREQ 


•) 


term «n( INPUT) ; 
tarmout ( OUTPUT ) ; 
CMS( ' CLRSCRN' . I ) ; 
ur«teln( 'Of  which 
wr i tel  n( ' 
wr i tal n( ' 
wrl taint ' 
wr 1 to1n( • 


file  ara  wa  taking  tha  OFT' 

1)  IOEAL '  )  ; 

2)  SIGNAL' ); 

2)  GRAPH' ); 

4)  Rxx ' ); 


raadtn(CHOICE); 

M  •  O; 

if  CHOICE  <>  2  than 


); 


begin 

case  CHOICE  of 

1  :  reeet (F Ilia. -name- IDEAL.POATA.*’ ); 

2  reeetCFlLEA. -naae-SIQNAL  POATA  •' ); 

4  reset ( Fill A . - nama-Rxx . POATA  • ' ) ; 

otnarwiaa 


ana; 

if  CHOICE  «>  4  tnen 

raaein(FiLIA,OUMMv,xR(  i  )).  ( «XR( i )-»Oummy  read  f i la  a i*a» ) 

repeat  (■  read  in  data  to  fft  •) 

H  •  N  ♦  1  ; 

raadi n( FI LEA. DUMMY ,XR(  N  )). 

XI(  N  >  •  0. 

until  eof(FlLlA); 


aiaa 

begin 

reeetIFiLCA.  name-GRAPH  poata  •). 
m  •  m  ♦  i . 

reed  i  n  ( f  i  l  E  A .  dummy  ,  x  i  (  u  l.xai  n  )).  (•  xkn>  is  S  Oummy  read*) 
III  N  )  •  0. 

until  eof ( F I L I A ) . 


a  •  o. 

TEST  • 


TIST/2 
1  . 


( *  calculate  m 


•  1 
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until  TEST  <■  1; 
if  TEST  <  1  than 

for  I  :■  N+1  to  RAISE(2,M)  do  (*  length  pad  with  0 

XR( . I . )  :■  O.O; 

N  :■  RAISE(2.M); 
wrl taln(n) ; 

for  L  :•  1  to  M  do  (*  calculate  FFT 

begin 

LE  RAISE(2.M+1-L): 

LEI  :•  LE  /  2.0; 

UR  :■  1; 

UI  :«  O; 

WR  :•  coe(PI/LE1); 

WI  :■  -e1n(PI/LE 1 ) ; 

for  J  :»  1  to  round(LEl)  do 

begin 

for  K  :■  round((d-1)/LE)  to  round(N/LE)  do 
begin 

I  :•  K»LE  ♦  d: 
if  I  <•  N  then 
begin 

IP  :■  I  ♦  round(LEI); 

TR  :•  XR( .1 . )  «•  XR(  .IP.  ); 

TI  :•  XI( .1 . )  ♦  XI( .IP.  ); 

TMR  :-XR( .1 . )  -  XR( .IP.  ); 

TMI  :•  XI(.I.)  -  XI (.IP.); 

XR( . IP .  )  TMR*UR  -  TMI»UI; 

XK.IP.)  :•  TMR*UI  ♦  TMI*UR; 

XR( . I.  )  :■  TR; 

XI ( . I . )  :•  TI; 
end; 

end:  (•  NEXT  K  •) 

TR  UR*WR  -  UI*WI; 

UI  :•  UR*WI  *  UI*WR; 

UR  :•  TR; 
end;  (•  NEXT  d  •) 
end;  (•  NEXT  L  •) 

NV2  :•  N  dlv  2; 

NM1  :•  N-1; 
d  :•  1; 

for  I  :■  1  to  NM1  do 
begin 

If  I  <  d  then 
begin 

TR  :•  XR( . I . ) ; 

XR( .1 .  )  :•  XR( .d.  ); 

XR( . d . )  :•  TR; 

TI  XI ( .  I  .  ) ; 

XI(.I.)  :•  XH.d.); 

XI(.d.)  TI; 
end; 

K  :■  NV2; 
while  K  <  d  do 
begin 

d  :•  d-K; 

K  :•  K  dlv  2; 
end; 

d  :•  d*K; 
end;  (•  NEXT  I  •) 

for  I  :■  1  to  N  do  (*  calc  mag  ft  phaae  of  OFT 

begin 

MA8(.I.)  :•  sqrt (aqr (XR( . I . ) )  ♦  eqr(XI( .1 . ))) ; 

PHASE ( . I . )  :■  0.0; 
if  XR( .1 . )  <>  0.0  then 

PHASE ( . I . )  :■  arctan(XI( . I . )/XR( .  I . ) )  •  180  /  PI; 

If  XR( .1 . )  <  0.0  then 
if  XI(.I.)  >  0.0  then 

PHASE ( . I . )  :■  180  ♦  PHASE(.I.) 

else 

PHASE ( . I . )  :•  - 


180  +  PHASE ( . I . ) 
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wrlta1n( 'Would  you  Ilka  tho  output  in  LOG  magnltuda  ( l-yas.O-no)?' ); 
raad)n(CHOICE); 

If  CHOICE  -  1  than 
bag  in 

LARGEST  : •  O.O; 
for  I  :•  1  to  N  do 

LARGEST  : ■  max(MAG( . I .). LARGEST) ; 
for  I  :>  1  to  N  do 

MAG( . I . )  :■  20* 1 og(MAG( . I . ) /LARGEST ) ; 

and; 

wrlta1n('Ia  thla  1)  an  X(w)'); 
wrltaln('  2)  a  Y(w)'): 

»rtt»ln('  or  alaa)  a  ragular  fraq  output'): 
raadln(CHOICE) ; 
casa  CHOICE  of 

1  :  rawrf ta(FILEA, 'nama*FREOX . POATA . * ' ) : 

2  :  rawr1ta(FILEA, 'nana*FREOV. POATA.*' ); 
otharwiaa 

rawr1ta(FILEA, 'nama*FREO. POATA.*' ); 

and; 

for  I  :■  1  to  min(N,l024)  do 

wr1taln(FILEA,I :4.MAG( .1 . ) : 15:8); 


and. 


of  FREO  •) 
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program  FILTER 


Author : 

Oat* 


Doug  Rldor 
10  July  86 


Purpoaa 


This  program  was  daslgnad  to  implement  a  generalized 
digital  flltar  according  to  tha  following  schama. 


Onca  a  window  slza  has  baan  chosan  tha  usar  is  asked  to 
supply  tha  raquastad  weighting  values  or  choose  a 
predefined  sat  for  either  tha  time  or  rank  weights. 


tinm  ordered  — 
1 nput  — >  | 

1  1  •• 

■  1  ! 

t 

t  t .  . 

t 

A(1)-0 

-0  -0.. 

. A(n)-0 

t 

t  t 

t 

|  SORT  | 

t 

t  t. . 

t 

B(1)-0 

-0  -0 . . 

.B(n)-0 

t 

T  VT 

t 

1 

♦  ♦ 

♦  |  •  ( 

The  program  then  opens 

the  external 

file  where  the 

values  of 

the  sampled 

data 


signal  to  be  filtered  are  and  initializes  tha  window. 
Initializing  consists  of  filling  tha  first  half  of  the 
window  with  the  first  data  value  and  then  reading  in 
subsequent  values  until  the  window  is  fll'led.  Next 
the  window  is  weighted  by  the  W( )  set  of  coefficients. 
The  window  is  then  rank  ordered.  After  this  the 
second  set  of  coefficients  are  applied  which  trim  off 
the  smallest  and  largest  values  of  the  data. 

After  this  last  trimming  the  remaining  samples  are 
summed  and  gain  adjusted.  Finally,  the 
window  is  stepped  over  one  value  in  time  by  deleting 
the  oldest  value  and  inserting  the  new  value  in  the 
proper  rank  order.  The  filtering  and  stepping 
continue  until  the  end  of  the  data  file  Is  reached. 


•) 

•) 

•) 

•) 

•) 

•) 

•) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

•) 

•) 

*) 

•) 

•) 

•) 

•) 

•) 

•) 

•) 

*) 

•) 

•) 

*) 


program  FILTER  ( INPUT, OUTPUT) ; 
^INCLUDE  CMS 


const 

MAXJXINDOW_SIZE  -  64; 


type 


DATA  ELEMENT  POINTER  ■  PDATA  ELEMENT; 
OATA  ELEMENT 


WEIGHT  ARRAY  TYPE 


record 

ACTUAL  VALUE 
MODIFIED  VAL 
NEXT  LARGEST 
NEXT~IN_TIME 
end; 

»  array ( . 1 . . MAX_WINDOW_SIZE . )  of  REAL; 


REAL; 

REAL; 

DATA  ELEMENT  POINTER; 
data"el  EMENT_PO INTER 


var 


a 

ft 


i 


* 

y 

¥ 

£ 

v 


s 


I 


V, 

V, 


4 


$ 

§ 

s 


P 
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SMALLEST  VALUE . 
OLOEST  VALUE 
DATA  FILE  SIZE. 

I. 

WINDOW  SIZE 
DUMMY." 

COOED  FILE  SIZE. 
FILTER  INPUT, 
FILTER- OUTPUT 
GRAPH  FILE. 

DATA  SAMPLES 
END  OF  OATA 
TIME  WEIGHT  ARRAY. 
RANK  WEIGHT  ARRAY 


DATA_ELEMENT— POINTER; 
INTEGER; 

REAL; 

TEXT; 

BOOLEAN; 

WEIGHT  ARRAY  TYPE; 


. . . . * . . 

C  •) 
(*  procedure  GET  PARAMETERS  •) 
(*  ”  •) 
(•  This  procedure  returns  the  user  input  persmeters  thet  will  be  •) 
(•  used  during  the  progrsm  run.  *) 
(•  •) 
. . . . . . * . . 


procedure  GET  PARAMETERS  (  vsr  WINDOW  SIZE  :  INTEGER; 

ver  TIME  WEIGHT  ARRAY  :  WEIGHT  ARRAY  TYPE; 
vsr  RANK- WEIGHT  ARRAY  :  WEIGHT- ARRAY  TYPE); 


ver 


OC. 

ANS. 

CHOICE. 

I 

CHECK. 
VALIO  SIZE 
GAIN. 

DUMMY 
COEFF  FILE 


:  INTEGER; 

:  BOOLEAN; 

;  REAL; 

:  TEXT; 


•) 

•) 

•) 

•) 

•) 

•) 


procedure  MANUAL.INPUT 

This  procedure  lets  the  user  input  the  weighting 
coefficients  to  be  used. 


procedure  MANUAL  INPUT<  ver  INPUT  ARRAY  :  WEIGHT  ARRAY  TYPE; 

WINOOW  SIZE  :  INTEGER); 


ver 

I. 

ANS  :  INTEGER; 
begin 

wr 1  tel n( ' Input  the  desired  filter  coef f icients  ' ) ; 
for  I  :■  1  to  round! WlNOOW_SI ZE/2 )  do 

begin 

write1n( 'A( ' , I ;2. ' )  •  7'); 
reedlnl input_aRRay<  I) ) ; 

end; 

writein(‘!s  the  filter  syeeetric  ( i*yes.O*no)?' ) ; 
reedln(ANS); 

for  I  :•  round! WINDOW  SIZf/2)*i  to  WINOOW  SIZE  do 
if  ANS  •  1  then 
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INPUT_ARRAY( .1 . )  : «INPUT_ARRAY( . WINOOW_SIZE  -1+1.) 

•  IM 

begin 

wrtte1n( 'A( ' ,1 :2, ' )  ■  ?'): 
rtadln( INPUT_ARRAY( . I . ) ) ; 

•nd; 


procedure  DO_ALPHA_TRIM 

Tnia  procedure  calculates  the  number  of  elements  to  be 
trimmed  off  of  each  and  of  the  rank  weight  array  and  0's  them. 


procedure  DO  ALPHA  TRIM  (  var  RANK  WEIGHT  ARRAY  :  WEIGHT  ARRAYJTYPE: 

WIND0W_SI2E~:  INTEGER); 

var 

I. 

NUM  TRIMMED. 

TRIM. 

NUMJ.EFT  :  INTEGER; 
begin 

writeln( ' Input  how  many  elements  to  trim  (must  be  even  a))'); 
read) n( NUM  TRIMMED); 

TRIM  :•  NUM  TRIMMED  dlv  2; 

NUM  LEFT  :•  WINDOW  SIZE  -  NUM  TRtMMEO; 
for” I  :■  1  to  TRIM  do 

RANK  WEIGHT  ARRAY ( . I . )  :•  O; 
for  I  :•  TRIM*"  to  TRIM*NUM  LEFT  do 
RANK  WEIGHT  ARRAY ( . I . )  :■  1; 
for  I  :■  TRIM* NUM  LEFT*1  to  WINOOW  SIZE  do 
RANK  WEIGHT  ARRAY ( . I . )  :•  0; 


procedure  CALCULATE 

This  procedure  calculates  the  gain  factor  for  the  filter 
coefficients  so  that  the  output  is  properly  scaled. 


•) 

•> 

•) 

•> 

•) 

*) 

•) 

•) 


procedure  CALCULATE  (  var  GAIN  :  REAL; 

WINDOW  SIZE  :  INTEGER; 

TIME  WEIGHT  ARRAY  :  WEIGHT  ARRAY  TYPE; 
RANk"wEIGHT— ARRAY  :  WEIGHT  ARRAY  TYPE); 


var 


A. 

d . 

1  :  INTEGER; 

OROER  :  WEIGHT  ARRAY  TYPE;  (•  this  is  the  TWA  ordered  •) 


OROBR(  .I  )  :•  TIME  WEIGHT  ARRAY (  . 1 .  ) ; 
for  I  :•  ]  to  WINOOW  SIZE~0O 

If  TIME_WEIGHT_ARRAY< . I . )  <  0RDER( . 1 . )  then 

begin 
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for  A  :■  (1-1)  downto  i  do 
ORDER (  .  A+1 . )  : ■  ORDER ( . A . ) ; 

ORDER ( . 1 . )  :■  TIME  WEIGHT  ARRAY ( . I . ) ; 


begin 

J  :■  2; 

while  (OROER( . J. )  <>  0)  and  (TIME  WEIGHT  ARRAY ( . I . )  > 

0RDER( . J ■ ) )  do 

J  :•  J  ♦  1; 

for  A  :■  (1-1)  downto  J  do 
0RD£R( .A+1 . )  : ■  0RDER( .A. ) ; 

0RDER( . J .  )  :•  TIME_WEIGHT_ARRAY( . I . ) ; 
and  ■ 

GAIN  0.0; 

for  1  :■  1  to  WINDOW  SIZE  do 

GAIN  :■  GAIN  ♦  0RDER(.I.)  •  RANK  WEIGHT  ARRAY ( . I . ) ; 


and; 


. . . 

(*  *) 

("  procedure  SHOW  *) 

(•  •) 

(•  This  procedure  clears  the  screen  and  than  provides  a  list  *) 

(*  of  the  coefficients  to  the  user.  *) 

(•  *) 

. . . 


procedure  SHOW  (  TIME  WEIGHT  ARRAY  :  WEIGHT  ARRAY  TYPE; 

RANK~WEIGHT“aRRAY  :  WEIGHT- ARRAY- TYPE ; 
WINDOW  SIZE-:  INTEGER); 


var 


I  :  INTEGER; 
begin 

CMS('CLRSCRN'.I); 

wr1teln< 'TIME  COEFF  VALUE  RANKCOEFF  VALUE'); 

for  I  :•  1  to  WINOOW  SIZE  do 

wr 1te1n( '  A('.I:2.')  -'.TIME  WEIGHT  ARRAY( . I . ) : 12:5. 

'  B(  '  .  1 : 2,  '  )  ■ '  ,  RANK_WEIGHT_>ARRAY (  .  I  .  ) :  12 : 5 ) ; 

end; 


. . . 

C  •) 

(*  procedure  MODIFY  *) 

(*  *) 

(*  This  procedure  modifies  either  the  whole  weight  array  or  *) 

(*  individual  coefficient  values  that  the  user  desires.  *) 

<«  *) 

(• - - - *) 


procedure  MOOIFY  (  var  MODIFY  ARRAY  :  WEIGHT  ARRAY  TYPE; 

WINOOW- SIZE  :  INTEGER?; 


var 


I  :  INTEGER; 

VAL  :  REAL; 

begin 

wrlteln( 'Which  element  (  0->  all  elements  )'); 
rescind ); 
if  I  ■  0  then 

MANUAL_INPUT ( MODIFY  ARRAY. WINDOW  SIZE) 

else 

begin 

wr 1  tel n( 'What  is  the  new  value'); 
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raadln(VAL) ; 

MOOIFY_ARRAY( .1 . )  :»VAL; 
and; 
and; 


bagin  (•  Of  GET_PARAMETERS  * 

CMS ( ' CLRSCRN ' . I ) ; 
wrltaln; 

wrlta1n( 'This  program  implamants  a  digital  flltar.'): 
wrltaln; 

rapaat 

wrlta1n( 'Input  tha  daslrad  window  alza(muat  bo  an  odd  intagar).') 
wrltaln; 

raad(WINOOW  SIZE); 

VALID  SIZE  :■  truo; 

If  not(odd(W!NOOW_SIZ£) )  than 
bag In  ~ 

wrlta1n('Tha  window  alza  must  bo  an  odd  intagar.'); 
VALID_SIZE  falsa; 
and: 

If  (WINOOW_$IZ£  >  MAX_WINDOW_SIZE )  or  (WINOOW_SIZE  <  0)  than 
bagl  n 

wrlto1n('Tha  slzo  must  ba  posltlva  but  lass  than 
MAX  WINDOW  SIZE); 

VALI0_S1ZE  :•  falsa? 
and;  " 

.  until  VALID.SIZE; 

wrlta1n( 'Would  you  Ilka  to  Initial iza  a  pradatarm 1 nad  flltar'); 
wr1ta1n( ' ( 1 .a.  a  moan  or  madlan  flltar)  or  manually  Input  tha'); 
w r1ta)n( 'fl it»r  coafflcianta?'); 
wrlta1n( '  i)  moan  flltar'); 

writo1n('  2)  madlan  flltar'); 

wr1to1n('  3)  alpha-tr immad  maan  flltar'); 

wrlta1n('  4)  alpha-tr Immad  llnaar  flltar'); 

wrlto1n('  5)  othar  coafflclant  input'); 

roadln(CHOICE); 

CMS( 'CLRSCRN' .1); 
casa  CHOICE  of 

1  :  for  I  :•  1  to  WIN00W_SIZE  do 
bogln 

TIME  WEIGHT  ARRAY ( . I . )  1/WIN00W  SIZE; 

RANK?WEIGHT?ARRAY(.I.)  1; 

and;  ”  ~ 

2:  bagin 

for  I  :>  1  to  WIND0W_SI2E  do 
bagin 

TIME  WEIGHT  ARRAY ( . I . )  :•  1; 

RANK_WEIGHT?ARRAY ( . I . )  : •  0; 
and; 

RANK_WEIGHT_ARRAY< . round( WIN00W_SIZE/2 ) . )  :•  1; 
and; 

3:  bagin 

for  I  :•  1  to  WINDOW  SIZE  do 

time_weight_array7.i . )  :■  1; 

00_ALPHA_TR IM( RANK_WE I GHT_ARRA Y . W I NOOW_S I Z  E ) ; 

CALCULATE(GAIN. WINDOW  SIZE. TIME  WEIGHT  ARRAY. RANK  WEIGHT  ARRAY); 
for  I  :•  1  to  WINOOW  SIZE  do 

RANK_WEIGHT_ARRAY7.I. )  :■  RANK_WEIGHT_ARRAY( . I . )/GAIN; 

and; 

4:  bagin 

CMS( 'CLRSCRN' .1); 
wrltaln; 
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wr1te)n( 'Alpha-trimmed  Itnur  filter'); 
wrltaln; 

wr1te1n('Are  tha  daslrad  flltar  coefficients  savad( 1-yes, O-no)?' ) ; 
read! n( ANS ) ; 
if  ANS  ■  1  than 
bag  In 

rasat(COEFF  FILE. ' naae»FIR . COEFF . * ' ); 
for  I  :■  1  to  WINDOW  SI2E  do 

raad1n(C0EFF_FILE7TIME_WEIGHT_ARRAY( .  I .  )  ) ; 
and  ” 

alsa 

begin 

CMS( 'CLRSCRN' , I ) ; 

MANUAL_INPUT(TIME_WEIGHT_ARRAY . WINDOW_SIZE ) ; 

wrlteln('Do  you  want  to  aava  this  aat  (1-yes,0-no)?'); 

raadln(ANS); 

if  ANS  -  1  than 

bagin 

rawrlta(COEFF  FILE. 'name-FIR .COEFF .•') ; 
for  I  :■  1  to  WINDOW  SIZE  do 

wr1te1n(C0£FF_FILE.TIME_WEIGHT_ARRAY( .1. )); 

and; 

and; 

CMS ('CLRSCRN '.I); 

00_ALPHA_TRIM( R ANK_WE I GHT_ ARRAY . WINDQW.SIZE ) ; 

CALCULATE (GAIN. WINOOW_SIZE . TIME_WEIGHT_ARRAY . RANK_WEIGHT_ARRAY ) ; 
wr4te)n('The  gain  for~this  sat  of  coefficients  ls',GAIN:7 :4) ; 
wrltaln; 

wr 1 ta1n( ' Is  this  a  flltar  that  paosaa  dc  ( 1-yas.O-no)?' ); 

>-aadln(OC); 
if  DC  ■  1  than 
bagin 

writaln('Tha  rank  walghtlng  eoafflciants  hava  baan  divided'); 
writeln('  by  tha  gain  in  order  to  normal iza  tha  de  level '); 
for  I  :■  1  to  WINDOW  SIZE  do 

RANK_WEIGHT_ARRAY(\l . )  :•  RANKJrfEIGHT_ARRAY( . I . )/GAIN; 

and; 

and; 

S:  begin 

wrlte1n('Ara  tha  daslrad  flltar  coefficients  saved( 1 -yes, 0-no )?') ; 
read ln( ANS); 
if  ANS  ■  1  than 
begin 

reset ( COE FF  FILE . ' nama-FIR . COEFF . * ' ) ; 
for  I  :•  1  to  WINDOW  SIZE  do 

readln(C0EFP_FILE7TIME_WEIGHT_ARRAY( .1. )); 

and 

else 

begin 

CMS( ' CLRSCRN' , I ) ; 
wrltaln; 

wr1ts1n('Now  input  tha  ' ,WINOOW_SIZE :2, '  weight  parameters  for  the'); 
wr1ta1n( 'time  ordered  buffer”); 

MANUAL.I NPUT ( T I ME_WE I GHT_ARR A Y . WI N00W_S I ZE ) ; 
wrlta1n('0o  you  want  to  save  this  «et”( l-yas,0-no)?' ) ; 
read)n(ANS) ; 
if  ANS  ■  1  than 
begin 

rewr1te(C0EFF  FILE. ' nama-FIR . COEFF .•') ; 
for  I  1  to  WINOOW  SIZE  do 

wrtte1n(C0EFF_FILE.TIME_WEIGHT_ARRAY( . I . )); 

and; 

and; 

CMS ('CLRSCRN'. I); 
wrltaln; 

wrlteln('Now  the  same  thing  for  the  rank  ordered  buffer'); 

MANUAL  INPUT (RANK  WEIGHT  ARRAY. WINOOW  SIZE); 
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CALCULAT E ( GA IN , WINDOW  SIZE. TIME  WEIGHT  ARRAY , RANK  WEIGHT  ARRAY); 
for  I  :«  1  to  WINDOW  SIZE  do 

RANK_WEIGHT_ARRAyT. I . )  :■  R ANK_W£ I GHT_ARRAY ( . I . )/GAIN; 

•nd; 


•nd;  (•  of  ease  statement  *) 

wr1t*ln( 'Would  you  Hkt  to  so*  tho  eooff fcl*nts( i-yes.O-no)?' ); 

roodln(I); 

if  I  ■  1  then 

SHOW ( T I ME„W£ I GHT_ARR A Y . RANK_WE I GHT_AR  R A Y . W I NOOW_S I Z  E ) : 

wr1t*ln( 'Would  you  11k*  to  chang*  anythlng( 1 -yes, O-no)?' ) : 

r*ad1n(ANS); 

whll*  ANS  ■  1  do 

begin 

wrlt*1n( 'Which  s*t  of  coefficients  1->tlme,  2->rank'); 
read1n( CHOICE ) ; 
case  CHOICE  of 

1  :  MODIFY (TIME  WEIGHT  ARRAY. WINDOW  SIZE); 

2  :  MODIFY ( RANK_WE IGHT~ARRAY ,WINDQW~SIZE ) ; 
end; 

SHOW  (TIME_WEIGHT_ARRAY . RANK_WEIGHT_ARRAY . WIND0W_SI2E ) ; 
wrlteln( 'Change  anything  *1s*( i-y*s7o-no)?' );  " 

readln(ANS) ; 
end; 

end;  (•  of  GETJ>ARAMETERS  *) 


(• 

(• 

<• 

(• 

(* 

(• 

<• 

(« 

(• 

(• 


procedure  INITIALIZE 

This  procedure  initializes  the  buffer  window  by  reading  the 
first  data  value  and  using  it  to  fill  the  first  half  of  the 
window  and  then  filling  the  rest  of  the  window  with  the 
subsequent  data  values  read  out  of  the  data  file. 


procedure  INITIALIZE 
var 


(  WINDOW  SIZE  :  INTEGER; 
var  OLDEST_VALUE  :  DATA_ELEMENT  POINTER): 


•) 

*) 

•) 

•) 

•) 

•) 

*) 

*) 

•) 

•) 


ELEMENT. 

NEXT  ELEMENT  :  DATA  ELEMENT  POINTER; 

DUMMY  :  REAL? 

I  . 

MIDDLE_OF_WINDOW  :  INTEGER; 

begin  <*  of  INITIALIZE  *) 

n*w( ELEMENT ) ;  (*  allocate  dynamic  variable  *) 

OLDEST  VALUE  :■  ELEMENT; 

r*ad(OATA_SAMPLES.OUMMY. ELEMENT*. ACTUAL  VALUE);  <*  read  1st  value  •) 
MIDDLE_OF_WINOOW  round(WINOOW_SIZE/27;  (*  find  middle  of  window  *) 

for  I  :■  2  to  MIDDLE J)F_WINDOW  do  (•  fill  1st  half  of  window  •) 
begin 

n*w(NEXT  ELEMENT); 

ELEMENT*7nEXT_IN_TIME  :■  NEXT  ELEMENT;  (*  link  data  values  •) 
ELEMENT  :■  NEXT  ELEMENT; 

ELEMENT*. ACTUAL.VALUE  OLDEST_VALUE* . ACTUAL_VALUE ; 
end; 

for  I  :•  (MIDDLE_OF_WINDOW  *  1)  to  WINDOW_SIZE  do  (*  read  In  data  *) 
begin  (*  for  last  half  of  window  •) 
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new(NEXT  ELEMENT); 

read (DATA  SAMPLES, DUMMY ,NEXT_ELEMENT*. ACTUAL  VALUE); 
ELEMENT#. NEXT  IN  TIME  :■  NEXT  ELEMENT; 

ELEMENT  :■  N£XT_ ELEMENT ; 
and; 


and;  (*  of  INITIALI2E  •) 

. . . . . . . 

(•  *) 

(•  procadura  WEIGHT  TIME  BUFFER  •) 

(•  ”  •) 
(*  This  procadura  walghts  aach  of  tha  elements  In  tha  tlma  odared  *) 

(*  buffar  by  tha  usar  input  values  storad  in  tha  TIME  WEIGHT  ARRAY*) 

(•  *) 

. . * . * . . 


procadura  WEIGHT  TIME  BUFFER  (  WINDOW  SIZE  :  INTEGER; 

OLDEST-VALUE  :  DATA  ELEMENT  POINTER: 

TIME  WEIGHT  ARRAY  :  WEIGHT  ARRAY  TYPE  ); 


var 


I  :  INTEGER; 

ELEMENT  :  OATA_ELEMENT_POINTER ; 

bag In  (•  of  WEIGHT_TIME_BUFFER  •) 

ELEMENT  :•  OLDEST  VALUE; 
for  I  1  to  WINOOW_SIZE  do 
bagtn 

ELEMENT#. MODI FI ED  VAL  :■  ELEMENT# . ACTUAL  VALUE  • 

TIME  WEIGHT  ARRAY( .1 .  ); 

ELEMENT  :•  ELEMENT# . NEXT_IN_tImE ; 
and; 


and; 


(•  Of  WEIGHT  TIME  BUFFER  •) 


procadura  RANK_ORDER 

This  procadura  taaas  tha  window  that  has  baan  modified  by  tha 
tlma  walghts  G  usas  an  insertion  sort  to  rank  order  aach  of  tne 
a laments.  At  tha  and  of  this  procedure  aach  element  in  tha 
window  is  ordered  by  size  of  tha  value  as  wall  as  in  time. 


procadura  RANK  OROER  (  WINDOW  SIZE  :  INTEGER; 

var  OLOEST_VALUE  :  DATA  ELEMENT  POINTER: 
var  SMALLEST_VALUE  :  DAT A_ELEMENT_POINTER ) ; 

var 

ELEMENT. 

NEW  ELEMENT  ;  DATA  ELEMENT  POINTER; 

I  :  INTEGER; 

begin  (*  of  RANK_ORDER  •) 

SMALLEST  VALUE  :-  OLDEST  VALUE;  (*  Start  with  first  value  as  •) 

SMALLEST**VALUE».NEXT  LARGEST  :*  nil;  (*  smallest  and  insert  each  •) 
NEW_ELEMENT  :•  OLDEST_VALUE# . NEXT_IN_TIME ;  (•  1st  new  valu  to  mart*) 
for  I  :■  2  to  wlNOOW_SIZE  do  '  ”  (•  loop  through  aach  value  •) 
begin  ~  (*  check  if  new  va1u<smanast* ) 

ELEMENT  :•  SMALLEST  VALUE;  (•  value  in  order  •) 

if  new_element#.modTfieo_val  <  SMALLEST_VALUE#.MOOIFIED_VAL 
than  begin  ~ 
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NEW  ELEMENT*. NEXT  LARGEST  :>  SMALLEST  VALUE; 
SMALLEST.VALUE  :•  NEW_ELEMENT ; 
end  ” 

■  1  a*  find  proper  order  •) 

begin 

while  (ELEMENT*. NEXT  LARGEST  <>  nil)  and 
(NEW  ELEMENT*. MOOIFIED  VAL  > 

ELEMENT*. NEXT  LARGEST*. MODIFIED  VAL)  do 
ELEMENT  :•  ELEMENT* . NEXT_LARGEST ; 

(•  and  ineert  in  place  *) 

NEW  ELEMENT*. NEXT  LARGEST  :>  ELEMENT*. NEXT  LARGEST; 
ELEMENT*. NEXTJ.ARGEST  :•  NEW_ELEMENT ; 
end; 

NEW_ELEMENT  NEW_ELEMENT* . NEXT_IN_YIME ;  (•  th*n  go  on  to  *) 

end;  ”  ~  T*  Tnaert  next  elmnt  in  order*) 

end;  (•  of  RANK_ORDER  •) 


. . * . * . . 

(•  *) 

(*  procedure  WEIGHT  RANK  BUFFER  *) 

(*  *) 

(*  Thle  procedure  weignts  each  of  the  elements  in  the  time  odered  *) 

(•  buffer  by  the  user  input  values  stored  in  the  TIME  WEIGHT  ARRAY*) 

(•  •) 

. . . . . . ••> 


procedure  WEIGHT  RANK  BUFFER  (  W1NOOW  SIZE  :  INTEGER; 

SMALLEST  VALUE  :  DATA  ELEMENT  POINTER; 
RANK  WEIGHT  ARRAY  :  WEIGHT  ARRAY  TYPE  ); 


var 


I  :  INTEGER; 

ELEMENT  :  OATA_ELEMENT_POINTER ; 


in 


(* 


Of  RANK  WEIGHT  BUFFER 


*) 


ELEMENT  :•  SMALLEST  VALUE; 
for  I  :•  1  to  WINOOW_SIZE  do 
begin 

ELEMENT*. MOOIFIED  VAL  ELEMENT*. MOOIFIED  VAL  • 

RANK  WEIGHT  ARRAyT.I): 
ELEMENT  :•  ELEMENT* . NEXT_LARGEST ; 
end; 


(•  of  RANK_WEIGHT_BUFFER  *) 


. . * . * . . 

(•  *) 

(•  procedure  FILTER  *) 

(•  *) 

(*  This  procedure  finds  the  middle  value  of  the  time  sequence.  *) 

(•  which  is  the  input  to  the  filter,  and  sums  the  values  of  the  *) 

(•  order  sequence,  the  filter  output,  and  sends  the  two  values  *) 

(*  back  to  tne  main  program  to  be  stored  in  a  graph  file.  *) 

(•  •) 

. . * . . 


procedure  FILTER  (  WINDOW  SIZE  :  INTEGER; 

ver  OLOEST'vALUE  :  DATA  ELEMENT  POINTER; 
var  SMALLEST  VALUE  :  DATA  ELEMENT  POINTER; 
ver  FILTER  INPUT  :  REAL; 
var  FILTER^OUTPUT  :  REAL); 

var 

I , 

MIOOLE_QF_WINDOW  :  INTEGER; 


a  -I  ,  ir 


4(m  .1  .( 
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INPUT  ELEMENT, 

OUTPUT_£LEMENT  :  DATA_ELEMENT_POINTER ; 


begirt 


Of  FILTER 


MIOOLE  OF  WINDOW  :■  round( WINDOW  SIZE/2); 

INPUT  ELEMENT  :>  OLDEST  VALUE? 

OUTPUT_EL£M£NT  :•  SMALLEST  VALUE; 

FILTER_OUTPUT  :•  0; 

for  I  :■  1  to  WINOOW_SIZE  do  (*  stop  thru  window  to  calc  *) 

begin  ”  (•  output  •) 

if  I  <  MIOOLE  OF  WINDOW  then  (•  filter  input  is  middle  val*) 

INPUT  ELEMENT”: •  INPUT  ELEMENT* . NEXT  IN  TIME; 

FILTER  OUTPUT  :•  FILTER  OUTPUT  ♦  OUTPUT~ELEMENT* . MOOI F I  ED  VAL; 
OUTPUT"ELEMENT  OUTPUT_ELEMENT* .NEXT_LARGEST ; 
end;  “  ” 

FILTER_INPUT  :■  INPUT_ELEMENT* . ACTUAL_VALUE ; 

9;  (•  Of  FILTER  •) 


procedure  STEP_ONE 

Step  the  window  over  one  data  value  by  deleting  out  the  oldest 
value  from  the  time  sequence  and  then  creating 
a  new  dynamic  variable,  reading  In  its  new  value,  and 
inserting  it  at  the  and  of  the  time  sequence. 


procedure  STEP_ONE  (var  OLDEST  VALUE  :  DATA  ELEMENT  POINTER; 

ENO_OF”OATA  :  BOOLEAN); 


DUMMY  :  REAL; 

ELEMENT. 

TARGET  ELEMENT  :  DATA  ELEMENT  POINTER; 


begin 


of  STEP  ONE 


(•  •  •  •  delete  oldest  value  from  window  •  •  •  •) 

TARGET  ELEMENT  :•  OLDEST  VALUE; 

OLO£ST?VALUE  :•  OLDEST_VALUE* . NEXT_IN_T IME ;  <•  delete  from  time  seq*) 


d 1 SpOSe( TARGET  JELEMENT ) ; 


(*  de-at locate  dynamic  varbl  *) 


(•  *  •  •  put  new  element  into  window  •  •  •  •) 
new( ELEMENT); 

TARGET_ELEMENT  :■  OLDEST_VALUE ;  (•  Start  at  beginning  to  find*) 

repeat  ~  (•  the  most  recent  value  •) 

TARGET  ELEMENT  :■  TARGET  ELEMENT* . NEXT  IN  TIME: 
until  TARGET  ELEMENT*. NEXT  In  TIME  •  nil;” 

TARGET  ELEMENT*. NEXT  IN  TIME  ?«  CLEMENT;  <•  add  new  value  •) 

ELEMENT*. NEXT_IN_TlMf  :•  nil;  (•  and  put  nil  value  on  end  •) 

if  not( ENO_OF_DATA )  then  (•  if  still  data  values  left  •) 

read(OATA_SAMPLES. DUMMY. element*. ACTUAL.VALUC )  (•  read  new  value*) 
else  (*  sTse  put  last  data  val  in  *) 

ELEMENT*. ACTUAL_VALUI  :•  TARGET_EiEMENT*  ACTUAL_VALUE ;  (•  buffer*) 


Of  STEP  ONE 


K 


£ 

$ 


K! 


I 

E 


) 
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MAIN  F ROOM AM 

Tha  Min  progran  First  initial  izm  tna  torotna !  For  input  an a 
output  and  than  gata  tha  nacaaaary  paraaatars  to  run  tna  prograo 
Froai  tha  uaar.  Tha  data  Fiia  ia  opanaa  and  ttitaraa  ana  tna 
raaulta  ara  storad  in  a  Fiia  to  ba  grapnaa  latar 


bag  in 

taraln(INPUT) ;  (•  opan  tarainai  For  input  • I 

taraout ( OUTPUT ) ;  (•  opan  tarainai  For  output  *1 

GET  PARAMETERS (WINDOW  SIZE. TIME  WEIGHT  ARRAY, RANK  WEIGHT  ARRAY) 
rasot ( DAT  A  SAMPLES. 'naaa«SIONAL-POATA . •' ) ; 
raadlnCDATA  SAMPLES. DUMMY. COO ED  FILE  SIZE); 

DATA_FILE_sTZE  round( 1000  •  C00ED~FILE_SIZE ) ; 
rawr7ta(GRAPH  FILE.  'naMaGRAPH.POATA~*  '  )  ;~ 

INITIALIZE ( WINDOW  SIZE. OLDEST  VALUE); 

ENO  OF  DATA  :•  FALSE; 

For”l  “•  1  to  OATA_FILE_SIZE  do 
bog  in 

WEIGHT  TIME  8UFFER(WIND0W  SIZE. OLDEST  VALUE. TIME  WEIOHT  ARRAY); 
RANK  ORDER (WINDOW  SIZE. OLDEST  VALUE . SMALLEST  VALUE); 

WEIGHT  RANK  BUFFER( WINDOW  SIZE , SMALLEST  VALUf.RANK  WEIGHT  ARRAY), 
filterTwindow  SIZE. OLDEST  VALUE. SMALLEST  VALUE. FILLER  INPUT. 
FILTER- OUTPUT  ) ; 

wr1taln(GRAPH— FILE. I . '  '.FILTER  INPUT.'  '. FILTER  OUTPUT ) ; 

IF  I  »  (DATA  FILE  SIZE  -  round(VINOOW  SIZE/2)  *-l )  than 
ENO  OF_DATA  :■  TRUE; 

STEP  ONE (OLDEST  VALUE . ENO  OF.DATa ) ; 


MLl  HU 
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»rarw  •*( 

Tnta  pry—  etiewiatM  tna  nontltiM  —  apuara  prror  ( NMSE  ) 
mnm  npr— 1  i—  avarapa  prror  (MAI)  bat upon  toe  data  flip*  Carp 
auat  ap  t— an  to  < naura  —ten  a  up  ia  tna  tdaai  flip  to  to  uaad 
—  tna  rpoprpne#  far  nor— I  t za 1 1 on  Tn i ■  prpy  —  also  panaretaa 
two  ftipa  to  y apn  lator  Ona  flip  eontatna  both  ngnai*  that 
top  orrpr  —a  ea'euiataa  apt naan  anp  tho  otnpr  ft  ip  eontatna 
thp  ptffpnpnpp  aotaoan  tnp  two  atgnpta 


pray—  HU  (  INPUT  .  OUTPUT  )  . 


xiNCmot  cat 


var 


I  . 

ouawt . 

CHOI ci. 

AMt . 

MUM  SOWLI!  INTEGER  ; 

ouaav  3 

IDEAL  VALUE  . 

ACTUAL  TALUS . 

oiff. 

•Ml  NOM. 

MAI  MOM. 

ma«T 

MMtl  IUL ; 

filea . 

f I  Lie . 

filed. 

I I  Lit  TEXT. 


■  tn 


tpnatnf INPUT) ; 
torpput ( OUTPUT ) ; 
CMS(  'CLRSCRM' . I ). 
wr « to i n( ' Ca i cu i ata 
**r  <  tpl n( ' 

»r i tpln( ' 
art  taint ' 
mr t  tal n( * 
art  taint ' 
rppplnt  CHOICE ) ; 

caao  choice  of 

t :  toogm 


tha  nmsc  ana  NAE  batwoon  thp'); 

1)  taoai  etgnai  ana  notaa  corruptod  atgnal'); 

2)  topal  atonal  and  ftitor  output'): 

3)  ftitar  input  and  output'); 

4)  idoai  and  taodtfiad  fraquancy  roaponaa'); 

5)  X(w)  and  Y(w)'  ); 


roaottFILEA,  na*a»IDEAL  POATA  .  •  '  ) ; 
raaat ( PI LEO .  naaa*S IQMAL  POATA . *' ) ; 


3 


boom 

raaat ( FI  LEA 
raaat (FILES 


'  nano* IDEAL . POATA . 
'  naMO*GRAPH . POATA . 


); 

); 


3 


4 


bao'n 

roaott  filea 
raaat (FILES 
and; 


na»a» SIGNAL . POATA . * 
'  na— •GSAPH  .  PDATA  .  •  ' 


) ; 


Oagin 

raaat ( FILEA 
raaat ( F I lEB 
and; 


nataa-IOEALSES  POATA  .  ); 

'  naaa-FREO. POATA  •' ); 


9 


tn 


raaattFlLEA 
raaat (FILES 


'  natM-FREOX  .  POATA  *  '  )  ; 
'  namasFREOY . POATA .  •  '  ) ; 


and; 

othprwtaa 

wrtta1n( 'That  ta  not  an  option'); 


mi:  MSI 
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NMSI  :■  0; 

NMSE  NORM  :»  0; 

NAE  T-  O: 

NAI  NORM  :•  0; 
if  CHOICE  <  4  than 

PMdl  n(  FI  LEA,  OUMMY  .DUMMY  2); 
NUM_SAMPLES  :■  round ( 1000  *  DUMMY_2); 


if  CHOICE  -  1  than 

raad l n<  F I LES . DUMMY . DUMMY_2 ) : 

rawr  1 ta( FILEC , ' nama*CHANG£ . PDATA . * ' ); 
rawr 1ta( FIlED. ' nama«DUALRESP . PDATA .• ' ); 
eaaa  CHOICE  of 

1:  for  1  :•  1  to  NUM_SAMPLES  do 
bagin 

raad1n(FILEA. OUMMY, IDEAL  VALUE); 
raadl h( F I LEB . OUMMY . ACTUAL  VALUE  ) ; 

OIFF  :•  IOEAL  VALUE  -  ACTUAL  VALUE; 

NMSE  :•  NMSE  ♦  sqr(OIFF); 

NMSE  NORM  :■  NMSE  NORM  ♦  sqr( IDEAL  VALUE); 

NAE  T-  NAE  *  abs(OIFF); 

NAE  NORM  :>  NAE  NORM  ♦  aba( IDEAL  VALUE); 
wr 1 ta 1 n( F I LEO . DUMMY : 4 . IDEAL  VALUE : 1 S : B . ACTUAL  VALUE 
wrltalni FILEC. DUMMY: 4. OIFF: 15:8); 

and; 

2 .  ( •  2  or  3  •) 

3:  for  I  :■  1  to  NUM_SAMPLES  do 
bagi  n 

raadl n< FI LEA. DUMMY, IOEAL  VALUE); 

raadl n( F I LEB . DUMMY . DUMMY  2 . ACTUAL, VALUE ) ; 

wr 1 ta 1 n<  F I LEO . DUMMY : 4 . IDEAL  VALUE : 15 : 8 . ACTUAL  VALUE 

DIFF  :•  IDEAL  VALUE  -  ACTUAL,VALUE ; 

wr 1 ta1n(FlLEC. OUMMY :4, OIFF: 15:8); 

NMSE  :■  NMSE  «■  sqr(OlFF); 

NMSE  NORM  :■  NMSE  NORM  ♦  sqr( IDEAL  VALUE); 

NAE  :■  NAE  ♦  abs(OIFF); 

NAE.NORM  :•  NAE_NORM  ♦  ab>( IDEAL_VALUE ) ; 
and; 

4 , ( •  4  or  5  •) 

5:  bagin 

wrltaln('Ara  tha  fraquancy  fllas  In  d8( 1-yaa,0-no' ) 
raadl n( ANS ) ; 

whlla  not (aof ( FILEA) )  do 
bagin 

raad1n( FILEA .DUMMY , IDEAL  VALUE); 
raadl n( F I LEB . DUMMY . ACTUAL  VALUE ) ; 
wr 1 ta 1 n( F I LED , OUMMY : 4 . I DE AL_V ALUE : 1 5 : 8 , ACTUAL.VALUE 
if  ANS  ■  1  than 
bagin 

IDEAL  VALUE  :■  axp< IDEAL  VALUE* 1 n( 10)/20 . 0) ; 
ACTUAL_VALUE  :■  axp( ACTUAL_VALUE*ln( 10)/20.0) 


:  1 5 : 8  )  ; 


15:8); 


OIFF  :•  IDEAL  VALUE  -  ACTUAL  VALUE; 

NMSE  :-  NMSE  *  sqr(DIFF); 

NMSE  NORM  :■  NMSE  NORM  *  sqr( IDEAL, VALUE ) ; 
NAE  ?-  NAE  *  abs(DXFF); 

NAE  NORM  :■  NAE  NORM  *  aba (IDEAL  VALUE); 
wr1taln( FILEC. DUMMY: 4, DIFF: 15:8); 
and; 
and; 


NMSE  :•  NMSE  /  NMSE  NORM; 

NAE  :■  NAE  /  NAE,NORM; 

wrltaln('Tha  normal Izad  MSE  is' .NMSE. 15:8); 
wr1ta)n('Tha  normal Izad  avaraga  orror  Is' .NAE : 15:8) ; 


aw 


Lv 

I 


I 

■ 
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( . . . . . . 

( . . . . . . . . . . 

(•  •) 

(*  program  MEGAMEDIAN  •) 

(•  *) 

(*  This  program  implements  a  sari as  of  procaduras  which  win  •) 

(*  craata  a  sat  of  data,  flltar  it  and  than  ganarata  tha  transfar  •) 

(*  function  for  tha  flltar.  Tha  program  win  loop  through  several  *) 

(*  runs  and  avaraga  tha  results  to  gat  a  battar  approximation  •) 

(*  to  tha  transfar  function  if  a  random  Input  to  tha  flltar  is  •) 

( *  usad .  • ) 

(*  •) 

(*  Tha  program  has  tha  following  flltar  typas  alraady  •) 

(•  implamantad:  •) 

( *  -  maan  f 1 1  tar  • ) 

(*  -  mad Ian  flltar  •) 

(*  -  alpha-tr immad  maan  flltar  (ATMP )  •) 

(*  and  tha  naw  flltar  wa  ara  trying  to  charactariza:  tha  *) 

(*  -  alpha-tr Immad  llnaar  flltar  (ATLP)  •) 

(•  •) 

<•  •) 

(•  PROGRAM  VARIABLES  :  •) 

(•  •) 

(*  SIGNAL  :  This  array  holds  tha  valuas  of  tha  signal  to  pa  •) 

( •  f 1 1 tarad .  • ) 

(*  •) 

(*  FFTX_MAG,  •) 

(*  FFTY_MAG  :  Tha  magnitudes  of  tha  PPT  of  tfia  input  signal  (X)  •) 

(•  ""  and  tha  flltar  output  (Y).  •) 

(•  •) 

(•  Tyx_MAG  :  Tha  transfar  function  of  tha  flltar  obtained  •) 

(•  ”  using  the  cross  power  spactrum  ■  Pxy/Pxx.  •) 

(.  .) 

<•  Tyx_SUM  :  Array  usad  to  sum  up  tha  above  variable  In  order  •) 

(*  to  average  tha  result  over  a  number  of  random  •) 

(•  inputs.  •) 

<•  •) 

(*  SIGNAL_SIZE  :  The  number  of  samples  taken  for  each  signal.  •) 

(•  •) 

(*  TRIALS  :  Tha  number  of  random  trials  to  avaraga.  •) 

<•  •) 

(«•  I.  •) 

("  J  :  Index  variables.  •  ) 

(•  •) 

(•  •) 

(*  In  this  program  tha  main  procedures  are  sat  off  by  a  single  line  *) 

(■  of  asterisks.  Sub-procedures  are  sat  off  by  a  single  line  of  •) 

(*  dashes.  Functions  or  other  small  procaduras  are  sat  off  by  a  •) 

(■  line  of  alternating  dashes  and  spaces.  *) 

(•  •) 

(.  main  procedures  ->  .....................  •) 

(•  •) 

(.  sub-proceures  ->  - - •) 

(*  •) 

(•  functions  or  other  •) 

(.  small  procedures ->  -  --  --  --  --  --  « ) 

(•  •) 

. . . . . . . . * . .....) 

( . . . . . . . . ...) 

program  MEGAMEDIAN  ( INPUT .OUTPUT  )  ; 

%INCLUDE  CMS 
type 

SIGNAL_TYPE  ■  array ( . -31 . . 1056 . )  of  REAL; 

MEOIAN_TYPE  ■  array( . 1 . . 128, 1 . .500. )  of  REAL; 

var 

SIGNAL, 

FFTX_MAG, 
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(try  MAS 

Tyx  MAS  *  :  SIGNAL  TYPE; 
Tyx~SUM  :  MEDIAN  TYPE; 
SIGNAL  SIZE . 

TRIALS. 

SEED  1 , 

SEE02. 

J. 

I  :  INTEGER; 

DUMMY  :  REAL ; 

DATE. 

TIMS  :  ALFA; 

CH  :  CHAR ; 


procedure  SORT 

This  procedure  is  used  to  sort  sn  arrsy  of  SIZE  eleeents  into 
order  It  is  found  here  at  the  beginning  of  tne  prograe  since 
it  is  used  in  both  tne  CALCULATE  procedure  ( m  FILTER)  to  sort 
tne  eleeents  of  tne  TIMS  jfRIGMTjkRRAV  into  order  so  that  tne 
gam  factor  eay  be  calculated  and  in  tne  MED  I  AN_AND_PR  I  NT 
procedure  to  sort  tns  SOO  elanents  of  tha  transfer  function  st 
each  of  2S4  points. 


procedure  SORT  (  SIZE  INTEGER; 

var  SORT_ARRAY  :  SIGNAL.TVPE  ); 


A. 

I  . 

J  :  INTEGER; 

OROER  :  SIGNAL JTYPS; 


for  I  : ■  1  to  SIZE  do 
OROER(  I  )  -  00; 

OROERC  1  )  •  SORT_ARRAV( . 1  ) ; 

for  I  •  2  to  SIZE  do 

If  SORT_ARRAV(  I  )  <  ORDER!  1)  then 
begin 

for  A  ■  (I-i)  down to  1  do 

0R0ER<  A* 1  )  <  ORDER (  A); 

ORDER!  1  >  : ■  SORT  ARRAY (  l  ) : 


else 

begin 

d  :  •  2: 

while  (ORDER*  d  )  <>  0)  and  (SORT  ARRAY ( .  I  )  > 

ORDER!  d  ) )  do 

d  :  •  d  ♦  1 ; 

for  A  .  •  (I-i)  down to  d  do 
OROER (  A* 1  )  a  0R0ER<  A  >; 

OROER (  d  )  : •  SORT_ARRAY ( . I . ) ; 


for  I  ■  1  to  SIZE  do 

SORT  ARRAY (I  )  •  ORDER)  I) ; 


procedure  makeoata 

This  procedure  generates  tne  input  sequences  to  be  filtered 
Many  different  signal  character  1st ics  are  possible  in  this 


K 


iO 

i 
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(*  procedure  They  arc  outlined  in  the  first  port  of  the  eotn  body  *) 


(•  of  tmo  procedure .  •) 
(•  •) 
. . I 


procedure  MAKIDATA  (  var  SIGNAL  SIGNAL_TYP|  ); 
oonet 

PI  •  3 . 14 1GGGGG3GGGO ; 

MAX_NUM_SAMPLCS  •  1034; 


var 


A. 

U. 

A. 

0. 

X. 

DUMMY  NUM 
INITIAL  VAL. 

SIGMA  SO. 

LIV1L , 

PP0G1MP . 

INPUL Sf  VALUC. 
stip  ofFsit . 

NO I Si  AMPLITUOC  RCAL; 

ANS. 


I  . 


0. 

FIlC  CMOICC, 
FRIO. 

CMOICC. 

START  FRIO. 

stop  Frio, 
star?  samp. 
stop^Iamp . 


NUMSTIPS . 

SIGN  INTIQIR; 


OONI  POOL I AN; 

0ATA_SAMPLIS  TIXT ; 


function  SON 

Thie  function  returns  the  sign  of  the  TIST_vALUI.  This 
is  not  avaiiadlo  as  an  intrinsic  Fsseat  function. 


function  SON  (  TtST  VALUI  :  RCAL  )  :  INTIQIR; 


SON  :>  1; 

if  TIST  VALUC  <  O  than 
SON  ;•  -1 ; 


in 


Of  MAKCDATA 


*) 


CMS( 'CLRSCRN' . I ) ; 

repeat 

writoln('How  many  samples  to  taka  (1024  recommended)? ' ) ; 
raadln(NUMSAMP); 

If  (NUNS AMP  >  MAX_NUM_SAMPLCS )  or  (NUMSAMP  <  1)  then 

vritoln('Tho  number  of  samples  must  do  oetveen  1  and  40S6'); 
until  (NUMSAMP  >  0)  and  (NUMSAMP  <■  MAX  NUM  SAMPLES); 

SIGNAL  SIZC  :•  NUMSAMP; 
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For  X  :»  >31  to  10M  do 
SIGNAL!. I.)  :■  0; 
wrlteln; 

writotal 'GUILD  A  SIGNAL' ); 

toot 

wrlteln! 'vnat  would  you  like  in  it?’); 


writotal ' 

1 ) 

Constant  lava)'); 

writotal ' 

2) 

Monotone  reaps ' ) ; 

writotal ' 

3) 

Sinusoids' ) : 

writotal ' 

4) 

Stops' ); 

writotal ' 

8) 

Gaussian  noise'); 

writotal ’ 

t) 

Lap loci an  noise'): 

wrl total ' 

?) 

IIR  filter  the  signal'): 

writotal ' 

•  ) 

That' ‘s  alt'); 

raadtal CHOICE ) ; 

eaee  CHOICE  of 

i  :  begin 

( *  eons  tant 

wri total 'vnat  snail  tho  constant  lava)  be?'); 
read! n! LEVEL): 
if  livel  <>  0  than 


fOT  I  :•  1  tO  NUMSAMP  00 
SIGNAL! . I . )  :•  LIVEL; 
and; 


2  :  Bag in  (•  r aap 

wrlteln! -wnat  ta  tna  starting  saapla  for  tha  raap?'); 
road 1 n( START_SAMP ) : 
wrlteln!  And'the  analog  saapla?'); 


reeaintSTOP.SAMP); 

writotal 'To  what  ralativa  waluo  should  tho  raap  r isa( fal 1 )? ' ) ; 
raadln(LEVCL); 

INITIAL  VAL  :•  SIGNAL! . ST ART_SAMP  ) ; 
for  I  :•  START  SAMP  tO  STQP_SAMP  do 

SIGNAL!  I)  ;•  INITIAL  VAL  ♦  Lf Vf L/ePSl STOP  5AMP-START_SAMP)» 

(I-START^SANP); 

and; 


•) 


3  :  bag in  (•  slnas  and  blocks  of  sinas  •) 

wrl total ' Input  fraquancy  of  slnuslod  in  HztO  if  group  daslrad)'); 
raadta(FRKO); 
if  FRIQ  <>  o  than 
begin 

writotal 'starting  saaplolO  if  all  saaplaa)?'); 

raadtal START  SAMP); 

if  START_SAMP  •  0  than 

begin  “ 

START  SAMP  :■  1; 

STOP_SAMP  :■  NUMSAMP; 

and 

also 

begin 

writotal 'ending  saapla?'); 
raadtal STOP_SAMP ) ; 

for  I  :>  START  SAMP  to  STOP  SAMP  do 

SIGNAL!. I.)  :•  SIGNAL!. l")  ♦  stn<2*PI»FRE0*I/NUMSAMP) ; 


also 

begin 

writotal 'Vhat  is  tha  starting  f requency? ' ) ; 
raadtal ST ART_FREO); 

writotal ' wnot  is  tha  last  f requency? ' ) ; 
raadtal STOP  FRIO); 

for  FREQ  :•  START  FRIO  to  STOP  FREO  do 
for  I  ;■  1  to  NUMSAMP  do 

SIGNAL!. I.)  :■  SIGNAL!. I.)  ♦  S 1n( 2*PI *FREO* I /NUMSAMP ) ; 


4  ;  begin  (• 

writotal 'How  many  stops  in  tna  signal?'); 


x'k  A  wV  -  % 


^  K  A  V'  A 


•taps  and  impulses 
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readln(NUMSTEPS); 
for  J  :•  1  to  NUMSTEPS  do 
begin 

uriteln! 'Input  the  starting  senile  for  step'.j); 
reed) n< START_SAMP ) ; 
uriteln! 'And  the  test  sample?'); 
reed)n(STOP_SAMP); 

Mr 1 tel n( 'And  the  offset?'); 

reedln(STEP  OFFSET); 

for  I  :•  START  SAMP  to  STOP  SAMP  do 

SIGNAL!. I.)  :•  SIGNAL!  .  I  .  )  ♦  STEP  OFFSET; 


8  :  begin  (•  Gaussian  noise  *) 

uriteln! ‘What  shall  the  noise  sigma  squared  be  (gaussian  diet.)?'); 

Mr i tel n( 'Sigma  squared?'); 

reedln(S!GMA  SO): 

for  I  :■  1  to  NUMSAMP  do 


begin 

u  :•  rendom(O) ; 

R  oqrt(-2.0*  SIGMA_SO  •  ln(u)); 
u  :•  random(O); 

0  :•  2.0  •  PI  •  U: 

X  :•  R  •  coa(O); 

SIGNAK.I.)  :■  SIGNAK.I.)  ♦  X; 


«  :  begin  (•  Laplaeian! impulsive)  noise*) 

vr tteln( 'Whet  shell  the  probability  of  an  impulse  be?'); 
reedln(PROSlMP) ; 
if  PROG IMP  >  o  then 
begin 

uriteln! 'Do  you  uent  1)  all  positive  impulses  or  2)  both  negative', 

'  end  positive  impulses. '); 
reedln!SIQN); 
if  SIGN  -  2  then 
SIGN  :•  -1; 

uriteln! 'wnat  shall  the  impulse  stenoerd  deviation  be?'); 
roedln! IMPULSE  VALUE); 
for  U  :■  1  to  NUMSAMP  do 
begin 

U  :•  rendom(O) ; 

if  U  <■  PROGIMP  then 

begin 

U  :•  rtndon(O) : 

x  :•  •! IMPULSE  VALUE/sqr t ( 2 ) ) *  1 n! 1 -U ) : 

SIGNAL!. J)  :■  SIGNAL!. J)  ♦  X  •  SON! random(O)  * 

(0.9  *  SIGN)); 


?  :  begin  (•  simple  HR  lou  pass  filter*) 

uri teln! 'What  is  the  parameter  ALPHA?'); 
reedln(A) ; 

for  J  :■  NUMSAMP  dovnto  2  do 

SIGNAK.I.)  :•  A*SIQNAL(  .1*1.  )  *  SIGNAK  .I  ); 

end; 


•  : 

DONE  :■  true; 

OTHERWISE  (•  do  nothing  •) 

ur iteln! 'That  is  not  an  option'); 


of  ease  statement 


CMS! 'CLRSCRN' , I ) ; 
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until  DONE; 

•nd;  (•  Of  MAKEDATA  •) 


. . . . . . . . . 

(•  *) 

(*  procadura  WINDOW  *) 

(•  *) 

(*  This  procadura  windows  tho  input  data  to  tha  flltar  with  tha  *) 

("  Blackman  window.  This  window  is  glvan  by  *) 

(•  *) 

(•  w(n)  ■  0.42  -  0.5cos(2nPl/N-i )  ♦  0.08cos(4nPl/N-l )  •) 

(•  •) 

( •  whara  n-i  is  tha  signal  siza.  •) 

<•  *) 

. . . . . 

procadura  WINDOW  (  var  SIGNAL  :  SIGNAL  TYPE  ); 


const 

PI  ■  3.141 5926535898 ; 


var 


MODIFIER  :  REAL; 

I  :  INTEGER; 

bag  in 

for  I  :■  1  to  SIGNAL_SIZE  do 
bagin 

MODIFIER  :•  0. 42-0. 5*COs(2»PI 'I/SIGNAL^Sm )«-0 . 08*Cos(4»PI»I/ 

SIGNAL  SIZE); 

SIGNAL(.I.)  SIGNAL( . I . )  •  MODIFIER; 
and; 


and; 


. . . . . 

(•  •) 

(•  procadura  FREO  •) 

(•  *) 

(*  This  procadura  takas  tha  dlscrata  Fourlar  transform  (OFT)  of  tha  *) 
(*  signal  givan  to  it.  Tha  OFT  is  implamantad  by  tha  dacimation-in-*) 

(•  fraquancy  FFT  algorithm.  •) 

(•  •) 

. . . 

procadura  FREQ  (  SIGNAL  :  SIGNAL  TYPE; 

var  FREQX  :  SIGNAL  TYPE); 


const 

PI  ■  3.1415926535898; 
MAX_ARRAY_SIZE  ■  2048; 

var 

DUMMY. 

I. 

J. 

M. 

N. 

L. 

NV2 , 

NM1 , 

IP. 

LE. 

K  INTEGER; 
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LARGEST. 

TEST. 

LEI. 

UR. 

Ul. 

MR. 

MX. 

TR. 

TI. 

TMR, 

TMI 

MAG. 

PHASE. 

XR. 

XI 


REAL; 


array( . 1 . . MAX_ARRAY_S I ZE . )  of  REAL; 


(•« 

(• 

<• 

(• 

(• 

(• 

<• 


function  RAISE 


*) 

*) 

•) 

■) 


This  function  ralaas  A  to  the  1 ntagar  powar  N.  This  function  is  *) 
not  intrinsically  availabla  in  tna  Pascal  eonpilar.  •) 

•) 


(> 


function  RAISE  (  A.N  :  INTEGER  )  :  INTEGER; 
var 


-) 


TEMP. 

I  :  INTEGER; 


bag  in 


TEMP  :■  1; 
for  I  ;•  1  to  N  do 
TEMP  :•  TEMP  •  A; 
RAISE  :•  TEMP; 


and; 


bagin 


(• 


of  FREQ 


N  :•  SIGMAL_SIZE ; 
for  I  :•  1  to  N  do 
bagin 

XR( .1. )  :■  SIGNAL( . I. ); 
XI(  .  I .  )  :•  0; 


( *  calculata  M 


M  :■  0; 

TEST  :•  N; 

rapaat 

TEST  :•  TEST/2; 

M  :•  M  ♦  1; 
until  TEST  <•  1; 
if  TEST  <  1  than 

for  1  :■  H*  1  to  RAISE(2,M)  do 
XR(  .1.  )  :■  0.0; 

N  :•  RAISE ( 2 , M) ; 
for  L  :•  1  to  M  do 
jin 

LE  :■  RAISE(2.M*1-L); 

LEI  :■  LE  /  2.0; 

UR  :■  1; 

UI  : ■  O; 

WR  :■  cos( PI/LE 1 ) ; 

«i  :■  -Sin(PI/LE 1 )  : 
for  J  :■  1  to  round(LEi)  do 
jin 

for  K  :•  round ( ( J- 1 )/LE )  to  round! N/LE)  do 


(•  langth  /»2**M,  pad  witn  O  •) 
(•  calculata  FFT  •) 


tl 


a 


i 


a 


YJ 


a 


9 


a 
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begin 

I  :■  K*LE  ♦  d; 

If  I  <•  N  then 
begin 

IP  :•  I  ♦  round(LEI); 

TR  :■  XR( . I . )  ♦  XR( . IP . ) ; 

TI  :■  XI(.I.)  ♦  XI ( . IP . ) ; 

TMR  :»XR(  .1 .  )  -  XR( . IP. ); 

TMI  :■  XI(  .  I .  )  -  XI (  .  IP .  )  ; 
XR( . IP . )  :■  TMR*UR  -  TMt»UI; 
XI (.IP.)  :■  TMR*UI  ♦  TMI-UR; 
XR( . I . )  :•  TR; 

XI(  .  I  .  )  :■  TI; 
end: 


end; 

(•  NEXT 

K  •) 

TR 

■  UR*WR  • 

UI*WI : 

UI 

•  UR*WI  ♦ 

UI*WR; 

UR 

-  TR; 

end;  (•  NEXT  d  •) 
end:  (•  NEXT  L  •) 

NV2  :■  N  dlv  2; 

NM1  :■  N-1; 
d  :■  1; 

for  1  :•  1  to  NM1  do 
begin 

if  I  <  d  then 
begin 

TR  :■  XR(.I); 

XR( . I . )  :■  XR( . d. ) ; 
XR( . d . )  :•  TR; 

TI  XI ( . 1 .  )  ; 
XI(.I.)  :•  X I (  d.  ); 
XI(.d.)  :•  TI; 
end; 

K  ;•  NV2; 
while  K  <  d  do 
begin 

d  :«  d-K; 

K  K  dlv  2; 
end; 

d  d-Mt; 
end;  (•  NEXT  I  •) 


for  I  :•  1  to  N  do  (•  celc  meg  A  pnese  of  DFT  ■) 

***^MAG( .  I . )  :•  sqrt(eqr<XR< . I . ))  ♦  aqr ( X I (  .  I .  ) ) ) ; 

PHASE( . I. )  :•  0.0; 
if  XR(  .1 . )  <>  0.0  then 

PHASE(.I)  :■  erctan(xi( . I . )/xR( . I . ) )  •  ISO  /  PI; 
if  XR(  .1. )  <  0.0  then 
if  xi(  .1 . )  >  0  0  then 

PHASE ( . I . )  :•  ICO  ♦  PHASE ( . I . ) 

else 

PHASE ( . I . )  :•  -1*0  ♦  PHASE!. I.) 

end; 


for  I  -  1  to  SIGNAL  SIZE  do 
FRE0X( .1 . )  :■  MAQ( .! . ); 


Of  FREQ 


•) 


procedure  FILTER 


Purpose  :  Tnis  progree  wes  designed  to  iwpieeent  a  generalised 
digital  filter  eccoroing  to  tne  following  sci 


Once  a  window  size  nas  been  chosen  the  user  is  asked  to 
supply  tne  requested  weighting  values  or  choose  a 
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predefined  set  For  either  the  tin*  or  rank  weights. 


time  ordered  ---- — -- — - - 

input  — >  |  |  |  .  | 

t 

t 

t . 

t 

A(1)-0 

-0 

-0. 

. . A(n)-0 

T 

T 

t 

t 

|  SORT 

T 

t 

t . 

T 

B(1)-0 

-0 

-0. 

.  B(n)-0 

t 

t 

t 

t 

1 

♦ 

♦ 

The  procoduro  then  opens  the  external  data 

file  where  the  values  of  the  seep led 

signal  to  be  filtered  are  and  initializes  the  window 

Initializing  consists  of  filling  the  first  neif  of  the 

window  with  the  first  data  value  and  then  reading  in 

subsequent  values  until  the  window  is  filiod.  Most 

the  window  Is  weighted  by  the  W( )  set  of  coefficients 

The  window  is  then  rank  ordered.  After  this  the 

second  set  of  coefficients  are  applied  wnicn  trio  off 

the  snallest  and  largest  values  of  the  data 

After  this  last  triaoing  the  renaming  aanpiea  are 

BuQwad  and  gain  adjusted.  Finally,  the 

window  is  stopped  over  one  value  in  tiee  by  deleting 

the  oldest  value  and  inserting  the  now  value  m  the 

proper  rank  order.  The  filtering  and  stopping 

continue  until  the  end  of  the  data  file  is  reached 


pi 


ocedure  FILTER  (  var  SIGNAL  :  SIGNAL  TYPE  ); 


const 


MAX_W1ND0W_SIZE  ■  B4; 


type 


DATA  ELEMENT  POINTER  •  00ATA  ELEMENT; 
DATA_ELEMENT  •  record 

ACTUAL  VALUE 
HQOIFlfo  VAL 
NEXT  LARGEST 

next_in_time 

end; 


REAL; 

REAL  . 

DATA  ELEMENT  POINTER. 
OATA  ELEMENT  POINTER 


WEIGHT  ARRAY  TYPE 


array!  1  MAX_WIN00W_S1ZS  >  of  REAL 


var 


TEMP 

SMALLEST  VALUE. 
OLDEST  VALUE 

I. 

WINDOW  SIZE 
FILTER  input. 
filter"output 

TIME  WEIGHT  ARRAY. 
RANK~WE I GHT~ ARRAY 


array!  1  1034  >  of  REAL . 

DATA_ELEMENT_POINTER ; 
INTEGER; 

REAL; 

WEIGHT  ARRAY  TVPE. 


(• . . . 

(•  •  > 

(•  procedure  GST  PARAMETERS  • > 

(•  •  » 
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(*  This  procadura  raturna  tha  uaar  input  paramatars  that  will  ba  *) 


(*  uaaa  during  tha  prograa  run.  •) 

C  *) 

. . *) 


procadura  GET  PARAMETERS  (  var  WINDOW  SIZE  :  INTEGER; 

var  TIME  WEIGHT  ARRAY  :  WEIGHT  ARRAY  TYPE; 
var  RANK~WEIGHT~ARRAY  :  WEIGHT~ARRAY~TYPE ) ; 

var 


ANS. 

CHOICE. 

1 

CHECK. 
VALID  size 
GAIN.* 
DUMMY 

COfFF  FILE 


INTEGER; 

•OOLEAN; 

REAL; 

TEXT; 


prooadura  NANUAL_ INPUT 

Tn«a  pracadu ra  lata  tna  uaar  input  tha  vaighting 
caafttcianta  to  oa  uaaa. 


pracadyra  MANUAL  INPUT(  var  INPUT  ARRAY  :  WEIGHT  ARRAY  TYPE: 

WlNDOW.SIZE  :  INTEGER); 


*) 

*) 

*) 

•) 

*) 

•) 

•) 

*) 


1 . 

ANS  INTEGER - 


•r>ta»«(  Input  tna  aaairad  t i i tar  coaf f icianta . ' ) ; 
*ar  I  •  »  to  '■ounai  wlKDOW_SIZE/2 )  do 

wr • tal n<  At  .1  2 .  1  •  ?  '  )  ; 

HMInl  INPUT  ARRAY  (  I  I); 


•trie  (1- 


,  0-no )? ' ) ; 


•r»ta<nt  la  tna  ti’tar  at 

•*ROd<  n|  ANS  > 

:  •  rpunol  WlNOOP_SIZE/2  )♦  1  to  WINOOW.SIZE  0O 

t  ♦ 


INPUT 


VI  I  )  •IIPUT  ARRAY!  WINDOW  SIZE 


I  ♦  1.) 


•r  •  ta>«<  A <  I  2 .  >  •  ?  '  ) , 
-aaatn)  input  array!  i  n. 


. . ) 

'  *  •) 

!  •  procadura  00  alpha  trip  •) 

(•  ■  •) 

i*  t*mo  prpoadura  ca'cu'ataa  tna  nuadar  of  a tenant a  to  ba  *) 

t  •  t->anaa  of*  of  oner  ana  of  tna  rank  waignt  array  and  0'a  than.  •) 

(*  *) 

.  . . ) 


00  Al*HA  trip  (  var  RANK  WEIGHT  ARRAY  :  WEIGHT  ARRAY  TYPE; 
WINDOW  SIZE*  INTEGER); 


*a  *_n  *  »  ’A  _a  a  "  a  *  A  •  .  •  .  '  .  .‘Y  _\  .w  .-V 
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NUM  TRIMMED. 
TRIM. 

NUM  LEFT 


INTEGER; 


begin 


wr1te1n( ' Input  how  many  clamant*  to  trim  (muat  ba  avan  «))'); 
raadln(NUM  TRIMMED); 

TRIM  NUM  TRIMMED  dlv  2; 

NUM__LEFT  : «~WINDOW_SIZE  -  NUM_TRIMMED; 
for  I  :■  1  to  TRIM  do  ~ 

RANK  WEIGHT  ARRAY ( . I . )  :■  0; 
for  I  ;■  TRIM+T  to  TRIM+NUM  LEFT  do 
RANK  WEIGHT  ARRAY ( . I . }  1; 

for  I  :■  TRIM+NUM  LEFT+1  to  WINDOW  SIZE  do 
RANK  WEIGHT  ARRAY ( . I . )  0; 


and; 


(•- 

(• 

(• 

(* 

(• 

<• 

(* 

C 


procadura  CALCULATE 

This  procadura  calculates  tha  gain  factor  for  tha  filter 
coefficients  so  that  tha  output  is  proparly  scaled. 

procadura  CALCULATE  (  var  GAIN  ;  REAL; 

WINDOW  SIZE  :  INTEGER; 

TIME  WEIGHT  ARRAY  :  WEIGHT  ARRAY  TYPE; 
RANK_WEIGHT_ARRAY  :  WEIGHTJkRRAYjrYPE ) ; 

var 


•) 

•) 

•) 

•) 

•) 

*) 

*) 

•) 


I 

fs 


I 


$ 

s 

Aj 

I 


A. 

d. 

I  ;  INTEGER; 

ORDER  :  SIGNAL_TYPE; 

begin 

for  I  :■  1  to  WINDOW_SIZE  do 

0RDER( . I . )  :>  TIME  WEIGHT  ARRAY( . I . ) ; 

SORT (WINDOW  SIZE. ORDER); 

GAIN  ;■  0.0; 

for  1  1  to  WINDOW  SIZE  do 

GAIN  :•  GAIN  ♦  ORDER(.I.)  *  RANK_WEIGHT_ARRAY( . I . ) ; 

and; 
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wr iteln( ' 

end: 


A( ' , I: 2 . '  )  -'.TIME  HEIGHT  ARRAY (  I  )  :  1 2  » , 
B( ' . I : 2 . '  )  ■' .RANK_WElGHT_lRRAY( . I . ) : 12:1) ; 


. . 

(• 

(•  proeAdurt  MOOIFY 

(• 

(•  This  procMurt  Modifies  elthsr  the  whole  weight  srrsy  or 
(•  individual  coefficient  values  that  the  user  desires. 

<• 

(•--- . * . 

procedure  MOOIFY  (  var  MOOIFY_ARRAY  :  HEIGHT  ARRAY_TYPE; 

HINOOH  SIZE  :  INTEGER); 


var 


I  :  INTEGER; 

VAL  :  REAL; 

begin 

writeln( 'Hhich  element  (  0->  all  eleeents  )'); 
read1n( I ) ; 
if  I  «  O  then 

MANUAL_INPUT(MOOIFY_ARRAY ,HINOOH_SIZE ) 

else 

begin 

wrlte1n( 'What  Is  the  new  value'); 
readln( VAL) ; 

MODIFY_ARRAY( . I  .  )  :■  VAL; 
end; 
end; 


begin 


(•  of  GET_PARAMETERS  •) 


CMS('CLRSCRN'.I); 

wrlteln; 

write1n( 'This  program  implements  a  digital  filter.'); 
wrlteln; 

repeat 

wrl tel n( ' Input  the  desired  window  size(must  be  an  odd  integer)  '); 
wrlteln; 

read ( W I NOOW_S I ZE ) : 

VALID  SIZE  true; 

If  not(odd(WINOOW_SIZE))  then 
begin 

wrlte1n('The  window  size  must  be  an  odd  integer.’); 
VALID.SIZE  ;■  false; 
end; 

If  (WINOOW_SIZE  >  MAX_WINDOW_SIZE)  or  (HINOOH_SIZE  <  O)  then 
begin 

wrlte1n('The  size  must  be  positive  but  less  than 
MAX_WINOOH_SIZE); 

VALIO_SIZE  :■  false; 
end; 

until  VALID  SIZE; 

for  I  :•  1  to  WINDOW_SI ZE  do 

begin 

TIME  HEIGHT  ARRAY (  .  I  .  )  :■  0; 

RANKJ<EIGHT_ARRAY( . I . )  : •  0; 

end; 


wrlte1n( 'Would  you  like  to  initialize  a  predetermined  filter'); 
wrlte1n( ' ( 1 .e.  a  mean  or  median  filter)  or  manually  input  the'); 
writeln( 'filter  coefficients?'); 


wrl tel n( 
write1n( 
wr 1 te1n( 
wrlteln( 
wrlte1n( 


D 

2) 

3) 

4) 

5) 


mean  filter'); 
median  f 1 1  ter' ) ; 
alpha-tr Immed  mean  filter'); 
alpha- trimmed  linear  filter') 
other  coefficient  input'); 
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road!n( CHOICE  ) ; 

CMS( 'CLRSCRN' , I ) ; 

CAM  CHOICE  of 

1  :  for  I  :•  1  to  WlNDOW_SIZE  Oo 
0^0 1  n 

TIME  WEIGHT  ARRAY (  I)  :•  t/WINDOW  SIZE; 
RANK_WEIGHT"ARRAY( . I . )  :•  1; 
ond; 

2 ;  0*010 

for  I  :■  1  to  WINDOW_SI ZE  Oo 

TIME  WEIGHT  ARRAY (  .  I  .  )  :•  1; 

RANKJfEIGHT_ARRAY( . I . )  :•  0; 
onj; 

RANK_WEIGHT_ARRAY(  round(WINDOW_SIZE/2)  )  ;■  1; 
ond; 

2:  bogto 

for  I  ;■  «  to  WINDOW  SIZE  OO 
TIMf_WSIOHTJUMUvT.  I  .  )  :•  1; 

00_ALPMA_TRIN(RANK_WEIGHT_ARRAY . WINDOW.SIZE ) ; 

CALCULATE ( GAIN , WINOOW  SIZE. TIME  WEIOHT  ARRAY, RANK  WEIGHT  ARRAY); 
for  I  •  1  to  WINOOW  ?IZE  OO 

RANK_WEIGHT_ARRAvT. I . )  : ■  RANK_WE IGHTJLRRAY (  . I .  )/GAIN; 

OHO; 

4 :  bog 1 n 

CMS( ' CLRSCRM ' , I ) ; 
writoln; 

writotn(  AlpOo-tr loood  1 inoor  flltor'); 
wr 1  tain; 

wr 1 ta I n( ' Aro  too  dostrod  flltor  cooffldonto  aavad( f-yoa. 0-no)7 ' ) ; 
roodlo(ANS) ; 
if  AMS  ■  1  thon 
Dog  in 

rooot(COEFF  FILE. ' naoa-FIR . COEFF  .  ); 

for  I.  1  to  WINOOW  SIZE  Oo 

road! n( C0EFF_F I LE . T IME_WE IGHT— ARRAY ( .1. )); 

ond 

•  1M 

***^CMS(  'CLRSCRN'  .1); 

MANUAL_INPUT(TIME_WEIGHT_ARRAV , WIN00W_SIZE ) ; 
writainf'Do  you  want  to  aavo  thia  aat  ( i-yoa,0-no)?' ) ; 
roadlo(ANS); 
if  ANS  •  1  than 
bog  in 

rour ltd (COEFF  FILE. ' naoo-FIR . COEFF  •  ') ; 
for  I  1  to  WINOOW  SIZE  do 

wr1tOln(C0EFF_FlLE.TIME_WEIQHT_ARRAY( .1 . )); 

ond; 

ond; 

CMS( 'CLRSCRN' .1); 

DO_ALRHA_TRIM(RANK_WEIGHT_ARRAY ,WIN00W_SIZE ) ; 

CALCULATE (GAIN. WI ND0W_S I ZE . T I ME_WE I GHT_ARR A Y , RANK_WEI GHT_ARR A Y ) ; 
wrltoln(’Tho  gain  for”thia  aat  of  coofficianta  ia“.0AIN;7 :4, 

'  and  tho  rank ' ) ; 

wr ito!n( 'weight ing  coofficianta  nova  boon  dividad  by  It.'); 
for  I  1  to  WINOOW  SIZE  do 

RANK_WEIGHT_ARRAyT. I . )  :•  RANK_WEIGHTjkRRAY( . I . )/QAIN; 

and; 

8:  bag  1 n 

writatn('Ara  tho  daairad  filtar  coofficianta  aavadf i-yaa.O-no)? • ) ; 
road! n( ANS ) ; 
if  ANS  •  1  than 
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bagin 

raaat (COEFF  FILE. 'nama«FIR.COEFF.»' ); 
for  I  1  to  WINDOW  SIZE  do 

raadln(COEFF_FILE7TlME_WEIGHT_ARRAY(  .  I .  ) ) ; 

•nd 

till 

bagin 

CMS( 'CLRSCRN' , I ) ; 
writaln; 

wrlta1n('Now  Input  tha  '  ,  WINDOW_SIZE  :  2 .  '  weight  parametera  for  t  ho'): 
wr1teln( ' tin*  ordered  buffer');- 

MANU A  L_ I NPUT (TIME-_WEIGHT_ARRAY,WINDOW_SIZE): 
wrlteln('Do  you  want  to  aave  thia  aat  ( 1-yaa.O-no)?' ) ; 
raadl n( ANS); 

If  ANS  -  1  than 
bagin 

rawrl ta(COEFF  FILE . 'nama-FIR .COEFF . • ' ) ; 
for  2  i-l  to- WINDOW  SIZE  do 

wr 1 teln(COEFF_FILE , TIME_WEIGHT_ARRAY( .!.)); 

and; 

and; 

CMS( 'CLRSCRN' ,1); 
writaln; 

wrlte1n('Now  tha  aane  thing  for  tha  rank  ordarad  buffar'); 

MANUAL..  I NPUT ( RANK_WE I GHT_ARR A  Y , W I NDOW_S I Z  E ) ; 

CALCULATE(GAIN. WINDOW  SIZE. TIME  WEIGHT_ARRAY ,RANK_WEIGHT  ARRAY); 
for  I  :■  1  to  WINDOW  SIZE  do 

RANK_WEIGHT_ARRAY7.I. )  :•  RANK_WEIGHT_ARRAY( . I . )/GAIN; 

and; 


and;  (*  of  caaa  atatamant  •) 

writa1n( 'would  you  Ilka  to  aaa  tha  coafflclanta( 1-yaa.O-no)?' ) ; 
raadl n( I ) ; 

If  I  ■  1  than 

SHOW( TIME_WEIGHT_ARRAY . RANK_WEIGHT_ARRAY . WINDOW_SIZE ) ; 

wr Ita1n( 'Would  you  Ilka  to  changa  anyth 1ng( 1-yaa.O-no)? ') ; 
raadl n( ANS);  • 
whila  ANS  ■  1  do 
bagin 

wr1ta1n( 'Which  aat  of  coaff Iolanta  1->t1ma,  2->rank'); 
raadl n( CHOICE ) ; 
caaa  CHOICE  of 

1  :  MODIFY (TIME  WEIGHT  ARRAY. WINDOW  SIZE); 

2  :  M0DIFY(RANk7wEIGHT7aRRAY.WIND0W_SIZE); 

and:  ” 

SHOW  (TIME_WEIGHT_ARRAY.RANK_WEIGHT_ARRAY,WINDOW_SIZE) ; 
writaln! 'Changa  anything  alaa( 1-yaaTo-no)?' ) ; 
raadl n( ANS); 
and; 


{*  Of  GET  PARAMETERS  •) 


( . ) 

(•  *) 

(*  procadura  INITIALIZE  *) 

(•  •) 

(•  Thia  procadura  initial izaa  tha  buffar  window  by  raadlng  tha  *) 

(*  firat  data  valua  and  uaing  it  to  fill  tha  flrat  half  of  tha  *) 

(*  window  and  than  filling  tha  raat  of  tha  window  with  tha  *) 

(*  aubaaquant  data  valuaa  raad  out  of  tha  data  flla.  *) 

(•  •) 

( . - . - . . 

procadura  INITIALIZE  (  WINOOW.SIZE  :  INTEGER; 


var  OLDEST  VALUE  :  DATA_ELEMENT_PO INTER ) : 


var 
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ELEMENT. 

NEXT  ELEMENT  :  DATA  ELEMENT  POINTER; 

I . 

DUMMY, 

MIDOL£_OF_WINOOW  ;  INTEGER; 

begin  (*  of  INITIALIZE  •) 

new! ELEMENT) ;  (*  allocate  dynamic  variable  •) 

OLDEST  VALUE  :■  ELEMENT; 

MIODLE  OF  WINOOW  :■  round!WINDOW  SIZE/2);  (•  find  middle  of  window  •) 
ELEMENT*. ACTUAL_VALUE  SIGNAL(“l -MIDOLE_OF_WINOOW+ 1 . ) ; 

for  I  :•  2  to  WIND0W_SIZE  do  (■  read  in  data  *) 

begin  (•  for  last  half  of  window  •) 

n*w(NEXT  ELEMENT); 

NEXT  ELEMENT*. ACTUAL  VALUE  :■  SIGNAL! . I -MIDDLE  OF  WINOOW-M.); 
ELEMENT*. NEXT  IN  TImI  :•  NEXT  ELEMENT; 

ELEMENT  :•  NEXT_ELEMENT ; 
and; 


•nd;  (•  of  INITIALIZE  •) 

. . . 

<*  •> 

!*  procedure  WEIGHT  TIME  BUFFER  •) 

C  “  •> 

(*  This  procedure  weights  each  of  the  elements  in  the  time  odered  *) 

!*  buffer  by  the  user  Input  values  stored  In  the  TIME  WEIGHT  ARRAY*) 

(*  •) 

. . . 


procedure  WEIGHT  TIME  BUFFER  (  WINDOW  SIZE  :  INTEGER; 

OLDEST  VALUE  :  DATA  ELEMENT  POINTER; 

TIME  WEIGHTJkRRAY  :  WEIGHT  ARRAY  TYPE  ); 


var 


I  ;  INTEGER; 

ELEMENT  :  DATA_ELEM£NT_POINTER; 
begin 

ELEMENT  ;■  OLDEST  VALUE; 
for  I  :■  1  to  WINOOW_SIZE  do 
begin 

ELEMENT*. MODI F I EO  VAL  ELEMENT* . ACTUAL  VALUE  • 

TIME  WEIGHT  ARRAY! . I .  ) ; 
ELEMENT  :■  ELEMENT* . NEXT_IN_TIME ; 
end; 


end; 


. . . 

(•  •) 

!*  procedure  RANK  ORDER  *) 

(•  “  *) 
!*  This  procedure  takes  the  window  that  has  been  modified  by  the  *) 

!*  time  weights  &  uses  an  insertion  sort  to  rank  order  each  of  the*) 

!•  elements.  At  the  end  of  this  procedure  each  element  in  the  *) 

!*  window  is  ordered  by  size  of  the  value  as  well  as  in  time.  *) 

C  *) 

<• . . 

procedure  RANK  ORDER  (  WINDOW  SIZE  :  INTEGER; 

var  OLDEST~VALUE  ;  DATA  ELEMENT  POINTER; 


var  SMALLEST_VALUE  :  DATA_ELEMENT_POINTER) ; 


var 


ELEMENT, 

NEW  ELEMENT  ;  DATA  ELEMENT  POINTER; 
I  _  :  INTEGER; 
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|ln  (•  Of  RANK.OROER  •) 

SMALLtST  VALUt  :•  OLDEST  VALUE;  (•  start  with  first  value  as  •) 

SMALLEST  VALUE*. NEXT  LARGEST  :•  ml;  (•  smallest  and  inaart  each  •) 

NtW.ILEMlNT  :•  OLOEST.VALUEa.NEXT_IN.TIME;  (•  1st  naw  valu  to  mart*) 
for-!  :•  2  to  WINOOW.SIZE  do  (*  loop  through  sach  value  •) 

Degin  ~  (*  cheek  if  naw  vatu<saattaat* ) 

ELEMENT  :•  SMALLEST  VALUE;  (■  value  in  order  •) 

If  NEW.ELEMENTa.MOOIPIEO.VAL  <  SMALLEST_VALUEa.MOOiriEO.VAL 
then  Degin 

NEW  ELEMENT*. NEXT  LARGEST  :•  SMALLEST  VALUE; 
SMALLEST.VALUE  :*“nEW_ELEMSNT; 
and 

else  (*  else  find  proper  order  •) 

^wnile  (ELEMENT*. NEXT  LARGEST  <>  nil)  and 
(NEW  ELEMENT* .M001FIE0.VAL  > 

ELEMENT*. NEXT  LARGEST*. MOOT  PIED  VAL)  do 
ELEMENT  :■  ELEMENT*  NEXT.LAROEST ; 

(•  and  insert  in  place  *) 

NEW  ELEMENT*. NEXT  LARGEST  :•  ELEMENT* . NEXT  LARGEST; 
ELEMENT*. NEXT.LAROEST  :■  NEW.ELEMENT ; 
end; 

NEW.ELEMENT  :■  NEW.ELEMENT* .NEXT  IN.TIME ;  (•  then  go  on  to  •) 
end;  ~  ~  Tm  Tnsert  newt  slant  in  order*) 

9;  (•  Of  RANK  ORDER  •) 


procedure  WEIGHT  RANK.BUPPER  •] 

•) 

This  procedure  weights  each  of  the  elements  in  the  time  odered  •) 
buffer  by  the  user  input  values  stored  in  the  TIME.WEIGHT.ARRAV* ] 


procedure  WEIGHT  RANK  BUFFER  (  WINDOW  SIZE  .  INTEGER; 

SMALLEST  VALUE  .  DATA  ELEMENT  POINTER: 
RANK  WEIGHT  ARRAY  :  WEIGHT  ARRAY  TYPE  ); 


X  :  INTEGER; 

ELEMENT  :  OATA.ELEMENT.POINTER; 


begin 


ELEMENT  SMALLEST  VALUE; 
for  I  :•  1  to  WINOOW.SIZE  do 
begin 

ELEMENT*. MODI FI ED  VAL  :■  ELEMENT* . MODI  FI ED  VAL  • 

RANK  WEIGHT  ARRAyT-I.); 
ELEMENT  :•  ELEMENT* .NEXT.LARGEST; 
end; 


procedure  FILTER.M 

This  procedure  finds  the  middle  value  of  the  time  sequence, 
which  Is  the  Input  to  the  filter,  and  sums  the  values  of  the 
order  sequence,  tha  filter  output,  and  sends  the  two  values 
back  to  the  main  program  to  be  stored  in  a  graph  file. 


procedure  FILTER.M 


(  WINDOW  SIZE  : 
var  OLDEST  VALUE 


INTEGER; 

DATA  ELEMENT  POINTER; 


l 


V-V.vjy 


E 


R 


I 


n 
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var  SMALLEST  VALUE  :  DATA  ELEMENT  POINTER; 
var  FILTER  INPUT  :  REAL; 
var  F I LTER- OUTPUT  :  REAL); 


var 


I. 

MIDDLE  OF  WINDOW  :  INTEGER; 

INPUT  ELEMENT. 

OUTPUT.ELEMENT  :  DATA_ELEMENT_POINTER; 
bag in  (•  of  FILTER_M  •) 


MIDDLE  OF  WINDOW  round(WINOOW  SIZE/2); 

INPUT  ELEMENT  :•  OLDEST  VALUE; 

OUTPUT  ELEMENT  :•  SMALLEST  VALUE; 

FILTERJJUTPUT  :•  0; 

for  I  T»  1  to  WINOOW_SIZE  do  (*  atop  thru  window  to  calc  •) 

bag in  "  (•  output  •) 

if  I  «  MIDDLE  OF  WINDOW  than  (•  f 11  tar  input  la  mlddla  vat*) 

INPUT  ELEMENT": ■  INPUT  ELEMENT* . NEXT  IN  TIME; 

FILTER  OUTPUT  :•  FILTER  OUTPUT  *  OUTPUT~ELEMENT*. MODIFIED  VAL: 
OUTPUT'ELEMENT  :■  OUTPUT_ELEMENT* . NEXT_LARGEST ; 
and; 

FILTER_INPUT  :■  INPUT_ELEMENT* . ACTUAL_VALUE ; 


(• 


of  FILTER_M 


•) 


P 


procadura  STEP_ONE 

Stop  tna  window  o var  on*  data  vatua  by  dal at ing  out  tna  oidaat 
vatua  fro*  tha  t laa  sequence  and  than  creating 
a  new  dynamic  variabt*.  reading  in  ita  new  value,  and 
inserting  it  at  the  and  of  the  tie*  sequence. 


ocadure  STEP  ONE  (var  OLDEST  VALUE  :  DATA  ELEMENT  POINTER; 

NEXT  VAL  .  INTEGER;" 

WINDOW  SIZE  :  INTEGER); 


var 


J , 

DUMMY  :  INTEGER; 

ELEMENT. 

TARGET_ELEMENT  :  DATA_ELEMENT_PO INTER ; 
begin  (•  of  STEP_ONE  •) 

(*  *  *  *  delete  oldest  valus  fro*  window  •  •  •  •) 


TARGET  ELEMENT  :■  OLDEST  VALUE; 

OLOEST_VALUE  ;■  OLOEST_valuE# . NEXT_IN_TIME ;  (*  delete  fro*  time  seq* ) 
dispose(TARGET_ELEMENT ) ;  (•  de-al locate  dynamic  varbl  *) 

(•  •  •  •  put  new  element  into  window  •  •  •  •) 
new (ELEMENT); 

TAR«ET_ELEMENT  :■  OLDEST_valu£ ;  (•  start  at  beginning  to  find*) 

for  u  "■  1  to  (WINDOW  SIZE*2)  do  (*  the  most  recent  value  *) 

TARGET  ELEMENT  :■  TARGET  ELEMENT*  NEXT  IN  TIME; 

TAR0IT  ELEMENT*. NEXT_IN_JIM?  :•  ELEMENT;  "  "  (•  add  new  value  •) 

ELEMENT* , NEXT  IN  TIME  :■  nil;  (•  and  put  nil  value  on  end  *) 

ELEMENT* . ACTUAL  VALUE  :•  SIQNALt . NEXT  VAL.) 


(•  Of  STEP  ONE  •) 
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begin  (•  of  FILTER  -  MAIN  PROORAM  •) 

GET_PARAMETERS( WINDOW  SIZE. TIME  WEIGHT  ARRAY, RANK  WEIGHT  ARRAY); 
!NITlALIZE(W!NOOW  SIZE. OLDEST  VALUE); 
for  I  :•  1  to  SIGNAL.SIZE  do 

WEIGHT  TIME  BUFFER (WINDOW  SIZE, OLDEST  VALUE. TIME  WEIGHT  ARRAY); 
RANK  OROER(W!NOOW  SIZE, OLDEST  VALUE . SMALLEST  VALUE); 
WEIGHT.RANK.SUFFE5( WINDOW  SIZE , SMALLEST  VALUE, RANK  WEIGHT  ARRAY); 

F l LTER_M( WINDOW  SIZE. OLDEST  VALUE . SMALLEST  VALUE. FILTER  INPUT, 
FILTER  OUTPUT): 

TEMP(.I.)  :•  FILTER  OUTPUT; 

If  I  <>  SIGNAL  SIZE  than 

STEP_ONE(OLDEST_VALUE . ( I+round(WINDOW_SIZE/2 ) ) . WINDOW_SIZE ) ; 

and; 

for  I  ;•  1  to  SIGNAL  SIZE  do 
SIGNAL( . I . )  : •  TEMP(.I.); 


(*  Of  FILTER  -  MAIN  PROGRAM  *) 


procoduro  TRANS.AND.STORE 

This  procoduro  tokos  tho  transfor  function  of  tha  filtar  by 
dividing  tho  oagnltudo  of  tha  OPT  of  tho  output  of  tha  filtar 
by  tho  oagnltudo  of  tho  OPT  of  tha  input  to  tho  filtar.  This 
rosult  is  than  storod  for  tha  nunbar  of  trials  attooptod  so 
that  tno  individual  transfor  function  of  aach  trial  is 
raoaobarad  for  further  processing. 


procoduro  TRANS.AND  STORE  (  PREQX ,  FREQY  ;  SIGNAL  TYPE; 

vor  Tyx  SUM  :  MEDIAN  TYPE; 

TRIAL  :  INTEGER  T; 


vor 

I  :  INTEGER; 

FREOT  :  REAL; 


bog  In 

for  I  :>  1  to  SIGNAL.SIZE  dtv  2  do 
FREOT  :•  FREOY( . I . )/FREOX( . I .  ) ; 

Tyx_SUM< . I. TRIAL. )  :■  Tyx_SUM( . I . TRIAL .  )  ♦  FREOT; 
and; 


and; 


••) 

•) 

•) 

•) 

•) 

•) 

•) 

•) 

•) 

'•) 


procoduro  MEDIAN.ANO.PRINT 

This  procoduro  sorts  tho  data  of  tho  previously  calculated 
transfor  functions  on  a  point  by  point  basis.  Tho  median  value 
of  each  point  is  selactad  as  tho  output  of  the  filter  and  stored 
in  a  fl Jo  for  further  graphing. 


procoduro  MEDIAN. ANO  PRINT  (  OROER  ARRAY  :  MEDIAN  TYPE; 

TRIALS  :  INTEGER  );“ 


PILE  :  MGA2 


PULL 
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Tyx  .  SIGNAL  TYPI ; 

AN*. 

NUN. 

1NOCX . 

iNoexa. 

INC. 

T^INOEX. 

d. 

!  :  INTCGIR; 

TUVALU! . 

LARGEST  T. 

NORM  T  :  REAL; 

PILE?  :  TEXT; 


function  LOG 

Thio  function  calculates  tha  base  10  logaritha  of  X.  The  Pascal 
eoMpitar  only  provides  the  natural  logaritha  as  an  Intrinsic 
function.  Tha  base  10  logaritha  is  calculated  by 

X  ■  in  X 


in  10 


function  LOG  (X : REAL ): REAL; 


*) 

•) 

*) 

•) 

*) 

•) 

•) 

•) 

•) 

•) 

•) 

•) 

•) 


begin 

if  x  •  00  then 
X  :■  000000001; 

LOG  :•  ln(X)/1n( 10.0) ; 
end; 


In 


(•  Of  MED I AN_AND_PR I NT 


•) 


INDEX  :•  SIGNAL  SI2E  dlv  3; 

IN0EX3  :•  SIGNAL.SIZE; 

for  I  :•  1  to  INDEX  do  (•  sort  each  point  of  Tyx  •) 

begin 

for  J  :■  1  to  TRIALS  do 

Tyx(.d)  :■  ORDER_ARRAY( . I , d . ) ; 

SORT (TRIALS. Tyx); 

for  d  :•  1  to  trials  do 

OROER_ARRAY( .I.d. )  :■  Tyx(.d); 

end; 

for  I  :■  1  to  INDEX  do  (*  choose  median  as  response*) 

Tyx( . I .  )  :«  OROER_ARRAY( .1 . round( TRIALS/3 ) . ); 
readln(ANS) ; 
if  ANS  ■  1  then 
begin 

LARGEST_T  :•  Tyx( . 1 . ); 

T  VALUE  ;•  LARGEST  T; 

ON0EX  :•  1; 

for  I  :•  1  to  INOEX  do  (•  find  largest  value  of  Tyx  *) 

begin  (*  for  normalization  *) 

LARGEST  T  :•  max (LARGEST  T.Tyx(.I.)); 

If  LARGE ST—T  <>  T— VALUE  then 
begin 

T  INOEX  :•  I; 

T_VALUE  :■  LARGE ST_T ; 
end; 
end; 

N0RM_T  :■  0; 

INC  ;■  rounddNOEX  *0.02);  (•  choose  interval  around  *) 

NUM  ;■  0;  (*  largest  value  for  normal iz*) 

for  I  :•  max(1.T_INDEX-INC)  to  m1n(T_INDEX*INC. INDEX)  do 
begin  ”  {•  calculate  norm  factor  •) 

N0RM_T  : •  N0RM_T  ♦  Tyx( . I . ) ; 


PILE:  MEGA2 


PULL 
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NUM  :•  NUM  ♦  1; 
and; 

NORM_T  :•  NORN_T/NUM; 

for  T  :•  1  to  TnOEX  do  (•  normal izo  Tyx  *) 

Tyx( . I . )  :•  20* LOG ( Tyx ( . I . )/NORM_T ) : 

and; 


reurltelFILET. ' name-MEDTRANS . POATA . • ' ) ; 
for  I  -  1  to  INDEX  do 

«M te1n(PILET, ( I/INOEX2) : 10: A. '  ' .Tyx( .1  .  )); 


and; 


procadura  AVERAGE^ANO_PRINT 

This  procadura  averages  tha  data  of  tna  previously  calculated 
tranafar  functions  on  a  point  by  point  basis.  Tha  average  value 
at  aacri  point  is  calculated  tha  output  of  tha  ftltar  and  storad 
in  a  file  for  further  graphing. 


procadura  AVERAGE  ANO  PRINT  (  ARR  :  MEDIAN  TYPE; 

TRIALS  :  INTEGER  ); 


var 


Tyx  S  :  SIGNAL  TYPE: 

ansT 

NUM. 

INDEX. 

IN0EX2. 

INC. 

T_INDEX, 

U. 

I  ;  INTEGER; 

T  VALUE. 

LARGEST  T. 

NORM  T  "  :  REAL; 

FILET  :  TEXT; 


. . . 

<•  •) 

(•  function  LOG  •) 

(•  •> 

(*  This  function  calculates  tha  base  10  logarithm  of  X.  Tha  Pascal  *) 
(•  compiler  only  provides  tha  natural  logarithm  as  an  intrinsic  *) 

(*  function.  Tha  base  10  logarithm  is  calculated  by  *) 

(•  •) 

(•  log  X  »  in  x  •) 

(•  .  •) 

(•  in  10  •) 

(•  •> 

. . . 

function  LOG  (X:REAL) .REAL; 


begin 

If  X  •  0.0  than 
X  :■  0.00000001; 

LOG  :•  ln(X)/ln( 10.0) ; 

and; 


begin  (*  Of  A VERAGE_AND_PR I NT  " ) 

INOEX  :•  SIGNAL  SIZE  dlv  2; 

IN0EX2  ;>  SIGNAL  SIZE; 
for  I  ;■  1  to  INOEX  do 


FILE:  ME GAS 


FULL 
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begin  (*  average  each  point  to  •) 

Tyx_5(.I.)  :•  0.0:  (•  calculate  Tyx  •) 

fo r"j  :■  1  to  TRIALS  do 

Tyx_S(.I.)  :•  Tyx_S( - 1 . )  ♦  ARR(.1.U.); 

end: 

readln(ANS) ; 
if  ANS  -  1  then 

begin  (*  ehooae  largeat  value  of  •) 

LARGEST  T  :■  Tyx  S(.1.):  (•  Tyx  to  normalize  •) 

T  VALUE" :■  LARGEST  T; 

T  INOEX  :•  1; 

for  I  1  to  INOEX  do 

LARGEST  T  max ( LARGEST  T , Tyx  HI.)); 

If  LARGEST_T  <>  T_VALUE  then 
begin 

T  INOEX  :•  I; 

T"VALUE  :•  LARGE ST_T ; 
end;" 
end: 

NORM  T  :■  O; 

INC  :■  rounddNDEX  •  0.02):  (*  calculate  normal  izat ion  •) 

NUN  :■  0:  (*  factor  over  range  near  max*) 

for  I  max( 1 ,T_IN0EX-INC)  to  m1n(T_INDEX* INC. INDEX)  do 
1  n 

NORM  T  :•  NORM  T  *  Tyx  S(.l.): 

NUN  ?-  NUN  ♦  l7 

end: 

N0RM_ T  N0RM_T/NUM; 

for  T  :■  1  to  INOEX  do  (*  normalize  Tyx  *) 

Tyx_S( . I . )  :•  2O*L0G(Tyx_S(.I.)/N0RM_T); 
end  " 

elae 

for  I  :•  1  to  INOEX  do 

Tyx_S(.I.)  :•  Tyx_S( . I . )/TRIAL5: 

rewr 1 te( F I LET . ' name- SUMTR ANS . POATA . * ' ) ; 
for  I  1  to  INOEX  do 

wrlteln(FILET, ( I/IN0EX2) : 10:8, '  * ,Tyx_S( . I . ) ) ; 


end; 


begin 


(  •  Of  ME GAMED I AN 


•) 


termout ( OUTPUT ) ; 

reaetC INPUT, 'name* TRIALS. MEGA .8' ) : 
read ln( TRIALS) : 
datet 1me(0ATE .TIME ) ; 

readatr(str(TIME ) , 1 : 2.CH. SEED1 : 2.CH.SEED2) ; 

I  :■  I *SEED1*SEED2*I*SEED1*SEED2 ; 

DUMMY  random(I); 

SIGNAL_SI2E  :■  256; 
for  I  :•  1  to  TRIALS  do 
begin 

reaet( INPUT, 'name* INOATA. MEGA. B' ); 

MAKEDATA( SIGNAL): 

WINDOW!  SIGNAL); 

FREQ! SIGNAL . FFTX_MAG) ; 

FILTER! SIGNAL); 

FRE0(SIGNAL.FFTY_MAG); 

TRANS_AN0_ST0RE! FFTX_MAG, FFTY_MAG,Tyx_SUM, I ) ; 
end; 

MEDIAN  ANO  PR I NT (Tyx  SUM, TRIALS); 

AVERAGE  ANO  PRINT! Tyx  SUM, TRIALS) ; 


Of  MEGAMEDIAN 


*) 


END 

J)A'Te. 
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